ZEPPI: Proteome-scale sequence-based evaluation of protein-protein interaction models.

We introduce ZEPPI (Z-score Evaluation of Protein-Protein Interfaces), a framework to evaluate structural models of a complex based on sequence coevolution and conservation involving residues in protein-protein interfaces. The ZEPPI score is calculated by comparing metrics for an interface to those obtained from randomly chosen residues. Since contacting residues are defined by the structural model, this obviates the need to account for indirect interactions. Further, although ZEPPI relies on species-paired multiple sequence alignments, its focus on interfacial residues allows it to leverage quite shallow alignments. ZEPPI can be implemented on a proteome-wide scale and is applied here to millions of structural models of dimeric complexes in the Escherichia coli and human interactomes found in the PrePPI database. PrePPI's scoring function is based primarily on the evaluation of protein-protein interfaces, and ZEPPI adds a new feature to this analysis through the incorporation of evolutionary information. ZEPPI performance is evaluated through applications to experimentally determined complexes and to decoys from the CASP-CAPRI experiment. As we discuss, the standard CAPRI scores used to evaluate docking models are based on model quality and not on the ability to give yes/no answers as to whether two proteins interact. ZEPPI is able to detect weak signals from PPI models that the CAPRI scores define as incorrect and, similarly, to identify potential PPIs defined as low confidence by the current PrePPI scoring function. A number of examples that illustrate how the combination of PrePPI and ZEPPI can yield functional hypotheses are provided.

Automated optimization of the solubility of a hyper-stable α-amylase.

Most successes in computational protein engineering to date have focused on enhancing one biophysical trait, while multi-trait optimization remains a challenge. Different biophysical properties are often conflicting, as mutations that improve one tend to worsen the others. In this study, we explored the potential of an automated computational design strategy, called CamSol Combination, to optimize solubility and stability of enzymes without affecting their activity. Specifically, we focus on Bacillus licheniformis α-amylase (BLA), a hyper-stable enzyme that finds diverse application in industry and biotechnology. We validate the computational predictions by producing 10 BLA variants, including the wild-type (WT) and three designed models harbouring between 6 and 8 mutations each. Our results show that all three models have substantially improved relative solubility over the WT, unaffected catalytic rate and retained hyper-stability, supporting the algorithm's capacity to optimize enzymes. High stability and solubility embody enzymes with superior resilience to chemical and physical stresses, enhance manufacturability and allow for high-concentration formulations characterized by extended shelf lives. This ability to readily optimize solubility and stability of enzymes will enable the rapid and reliable generation of highly robust and versatile reagents, poised to contribute to advancements in diverse scientific and industrial domains.

Synthesis, Structure and Antibacterial Activity of Two Novel Coordination Polymers Based on N,N'-bis(4-carbozvlbenzvl)-4-aminotoluene and Heterocyclic Ligand against S. aureus.

As the use of antibiotics increases, the increasing resistance of bacteria is the main reason for the reduced efficiency of antibacterial drugs, making the research of new antibacterial materials become new hot spot. In this article, two novel coordination polymers (CPs), namely, [Cd2(L)2(bibp)2]n (1) and [Ni(L)(bib)]n (2), where H2L = N,N'-bis(4-carbozvlbenzvl)-4-aminotoluene, bibp = 4,4'-bis(imidazol-1-yl)biphenyl, and bib = 1,3-bis(1-imidazoly)benzene, have been synthesized under solvothermal and hydrothermal condition. Structural clarification was performed through infrared spectrum and single-crystal X-ray diffraction analysis, while thermal analysis and XRD technology were used for the performance assessment of compounds 1 and 2. In addition, antibacterial performance experiments showed that compounds 1 and 2 have certain selectivity in their antibacterial properties and have good antibacterial properties against S. aureus. As the concentration of the compound increases, the inhibitory effect gradually strengthens, and when the concentration of the compound reaches 500 µg/mL and 400 µg/mL, the concentration of the S. aureus solution no longer increases and has been completely inhibited.

Carbonic Anhydrases: Different Active Sites, Same Metal Selectivity Rules.

Carbonic anhydrases are mononuclear metalloenzymes catalyzing the reversible hydration of carbon dioxide in organisms belonging to all three domains of life. Although the mechanism of the catalytic reaction is similar, different families of carbonic anhydrases do not have a common ancestor nor do they exhibit significant resemblance in the amino acid sequence or the structure and composition of the metal-binding sites. Little is known about the physical principles determining the metal affinity and selectivity of the catalytic centers, and how well the native metal is protected from being dislodged by other metal species from the local environment. Here, we endeavor to shed light on these issues by studying (via a combination of density functional theory calculations and polarizable continuum model computations) the thermodynamic outcome of the competition between the native metal cation and its noncognate competitor in various metal-binding sites. Typical representatives of the competing cations from the cellular environments of the respective classes of carbonic anhydrases are considered. The calculations reveal how the Gibbs energy of the metal competition changes when varying the metal type, structure, composition, and solvent exposure of the active center. Physical principles governing metal competition in different carbonic anhydrase metal-binding sites are delineated.

Adenosine Encapsulation and Characterization through Layer-by-Layer Assembly of Hydroxypropyl-ß-Cyclodextrin and Whey Protein Isolate as Wall Materials.

Adenosine, as a water-soluble active substance, has various pharmacological effects. This study proposes a layer-by-layer assembly method of composite wall materials, using hydroxypropyl-ß-cyclodextrin as the inner wall and whey protein isolate as the outer wall, to encapsulate adenosine within the core material, aiming to enhance adenosine microcapsules' stability through intermolecular interactions. By combining isothermal titration calorimetry with molecular modeling analysis, it was determined that the core material and the inner wall and the inner wall and the outer wall interact through intermolecular forces. Adenosine and hydroxypropyl-ß-cyclodextrin form an optimal 1:1 complex through hydrophobic interactions, while hydroxypropyl-ß-cyclodextrin and whey protein isolate interact through hydrogen bonds. The embedding rate of AD/Hp-ß-CD/WPI microcapsules was 36.80%, and the 24 h retention rate under the release behavior test was 76.09%. The method of preparing adenosine microcapsules using composite wall materials is environmentally friendly and shows broad application prospects in storage and delivery systems with sustained release properties.

The Structural and Molecular Mechanisms of Mycobacterium tuberculosis Translational Elongation Factor Proteins.

Targeting translation factor proteins holds promise for developing innovative anti-tuberculosis drugs. During protein translation, many factors cause ribosomes to stall at messenger RNA (mRNA). To maintain protein homeostasis, bacteria have evolved various ribosome rescue mechanisms, including the predominant trans-translation process, to release stalled ribosomes and remove aberrant mRNAs. The rescue systems require the participation of translation elongation factor proteins (EFs) and are essential for bacterial physiology and reproduction. However, they disappear during eukaryotic evolution, which makes the essential proteins and translation elongation factors promising antimicrobial drug targets. Here, we review the structural and molecular mechanisms of the translation elongation factors EF-Tu, EF-Ts, and EF-G, which play essential roles in the normal translation and ribosome rescue mechanisms of Mycobacterium tuberculosis (Mtb). We also briefly describe the structure-based, computer-assisted study of anti-tuberculosis drugs.

Theoretical Study of Cyanidin-Resveratrol Copigmentation by the Functional Density Theory.

Anthocyanins are colored water-soluble plant pigments. Upon consumption, anthocyanins are quickly absorbed and can penetrate the blood-brain barrier (BBB). Research based on population studies suggests that including anthocyanin-rich sources in the diet lowers the risk of neurodegenerative diseases. The copigmentation caused by copigments is considered an effective way to stabilize anthocyanins against adverse environmental conditions. This is attributed to the covalent and noncovalent interactions between colored forms of anthocyanins (flavylium ions and quinoidal bases) and colorless or pale-yellow organic molecules (copigments). The present work carried out a theoretical study of the copigmentation process between cyanidin and resveratrol (CINRES). We used three levels of density functional theory: M06-2x/6-31g+(d,p) (d3bj); ωB97X-D/6-31+(d,p); APFD/6-31+(d,p), implemented in the Gaussian16W package. In a vacuum, the CINRES was found at a copigmentation distance of 3.54 Å between cyanidin and resveratrol. In water, a binding free energy ∆G was calculated, rendering -3.31, -1.68, and -6.91 kcal/mol, at M06-2x/6-31g+(d,p) (d3bj), ωB97X-D/6-31+(d,p), and APFD/6-31+(d,p) levels of theory, respectively. A time-dependent density functional theory (TD-DFT) was used to calculate the UV spectra of the complexes and then compared to its parent molecules, resulting in a lower energy gap at forming complexes. Excited states' properties were analyzed with the ωB97X-D functional. Finally, Shannon aromaticity indices were calculated and isosurfaces of non-covalent interactions were evaluated.

Human Plasma Butyrylcholinesterase Hydrolyzes Atropine: Kinetic and Molecular Modeling Studies.

The participation of butyrylcholinesterase (BChE) in the degradation of atropine has been recurrently addressed for more than 70 years. However, no conclusive answer has been provided for the human enzyme so far. In the present work, a steady-state kinetic analysis performed by spectrophotometry showed that highly purified human plasma BChE tetramer slowly hydrolyzes atropine at pH 7.0 and 25 °C. The affinity of atropine for the enzyme is weak, and the observed kinetic rates versus the atropine concentration was of the first order: the maximum atropine concentration in essays was much less than Km. Thus, the bimolecular rate constant was found to be kcat/Km = 7.7 × 104 M-1 min-1. Rough estimates of catalytic parameters provided slow kcat < 40 min-1 and high Km = 0.3-3.3 mM. Then, using a specific organophosphoryl agent, echothiophate, the time-dependent irreversible inhibition profiles of BChE for hydrolysis of atropine and the standard substrate butyrylthiocholine (BTC) were investigated. This established that both substrates are hydrolyzed at the same site, i.e., S198, as for all substrates of this enzyme. Lastly, molecular docking provided evidence that both atropine isomers bind to the active center of BChE. However, free energy perturbations yielded by the Bennett Acceptance Ratio method suggest that the L-atropine isomer is the most reactive enantiomer. In conclusion, the results provided evidence that plasma BChE slowly hydrolyzes atropine but should have no significant role in its metabolism under current conditions of medical use and even under administration of the highest possible doses of this antimuscarinic drug.

Structural basis for the regulation of plant transcription factor WRKY33 by the VQ protein SIB1.

The WRKY transcription factors play essential roles in a variety of plant signaling pathways associated with biotic and abiotic stress response. The transcriptional activity of many WRKY members are regulated by a class of intrinsically disordered VQ proteins. While it is known that VQ proteins interact with the WRKY DNA-binding domains (DBDs), also termed as the WRKY domains, structural information regarding VQ-WRKY interaction is lacking and the regulation mechanism remains unknown. Herein we report a solution NMR study of the interaction between Arabidopsis WRKY33 and its regulatory VQ protein partner SIB1. We uncover a SIB1 minimal sequence neccessary for forming a stable complex with WRKY33 DBD, which comprises not only the consensus "FxxhVQxhTG" VQ motif but also its preceding region. We demonstrate that the ßN-strand and the extended ßN-ß1 loop of WRKY33 DBD form the SIB1 docking site, and build a structural model of the complex based on the NMR paramagnetic relaxation enhancement and mutagenesis data. Based on this model, we further identify a cluster of positively-charged residues in the N-terminal region of SIB1 to be essential for the formation of a SIB1-WRKY33-DNA ternary complex. These results provide a framework for the mechanism of SIB1-enhanced WRKY33 transcriptional activity.

Structural basis and synergism of ATP and Na+ activation in bacterial K+ uptake system KtrAB.

The K+ uptake system KtrAB is essential for bacterial survival in low K+ environments. The activity of KtrAB is regulated by nucleotides and Na+. Previous studies proposed a putative gating mechanism of KtrB regulated by KtrA upon binding to ATP or ADP. However, how Na+ activates KtrAB and the Na+ binding site remain unknown. Here we present the cryo-EM structures of ATP- and ADP-bound KtrAB from Bacillus subtilis (BsKtrAB) both solved at 2.8 Å. A cryo-EM density at the intra-dimer interface of ATP-KtrA was identified as Na+, as supported by X-ray crystallography and ICP-MS. Thermostability assays and functional studies demonstrated that Na+ binding stabilizes the ATP-bound BsKtrAB complex and enhances its K+ flux activity. Comparing ATP- and ADP-BsKtrAB structures suggests that BsKtrB Arg417 and Phe91 serve as a channel gate. The synergism of ATP and Na+ in activating BsKtrAB is likely applicable to Na+-activated K+ channels in central nervous system.

Identification of domains in Plasmodium falciparum proteins of unknown function using DALI search on AlphaFold predictions.

Plasmodium falciparum, the causative agent of malaria, poses a significant global health challenge, yet much of its biology remains elusive. A third of the genes in the P. falciparum genome lack annotations regarding their function, impeding our understanding of the parasite's biology. In this study, we employ structure predictions and the DALI search algorithm to analyse proteins encoded by uncharacterized genes in the reference strain 3D7 of P. falciparum. By comparing AlphaFold predictions to experimentally determined protein structures in the Protein Data Bank, we found similarities to known domains in 353 proteins of unknown function, shedding light on their potential functions. The lowest-scoring 5% of similarities were additionally validated using the size-independent TM-align algorithm, confirming the detected similarities in 88% of the cases. Notably, in over 70 P. falciparum proteins the presence of domains resembling heptatricopeptide repeats, which are typically involvement in RNA binding and processing, was detected. This suggests this family, which is important in transcription in mitochondria and apicoplasts, is much larger in Plasmodium parasites than previously thought. The results of this domain search provide a resource to the malaria research community that is expected to inform and enable experimental studies.

Cyanobacterial α-carboxysome carbonic anhydrase is allosterically regulated by the Rubisco substrate RuBP.

Cyanobacterial CO2 concentrating mechanisms (CCMs) sequester a globally consequential proportion of carbon into the biosphere. Proteinaceous microcompartments, called carboxysomes, play a critical role in CCM function, housing two enzymes to enhance CO2 fixation: carbonic anhydrase (CA) and Rubisco. Despite its importance, our current understanding of the carboxysomal CAs found in α-cyanobacteria, CsoSCA, remains limited, particularly regarding the regulation of its activity. Here, we present a structural and biochemical study of CsoSCA from the cyanobacterium Cyanobium sp. PCC7001. Our results show that the Cyanobium CsoSCA is allosterically activated by the Rubisco substrate ribulose-1,5-bisphosphate and forms a hexameric trimer of dimers. Comprehensive phylogenetic and mutational analyses are consistent with this regulation appearing exclusively in cyanobacterial α-carboxysome CAs. These findings clarify the biologically relevant oligomeric state of α-carboxysomal CAs and advance our understanding of the regulation of photosynthesis in this globally dominant lineage.

Structure-property modeling of coumarins and coumarin-related compounds in pharmacotherapy of cancer by employing graphical topological indices.

Coumarins, a subgroup of colorless and crystalline oxygenated heterocyclic compounds originally discovered in the plant Dipteryx odorata, were the subject of a recent study investigating their quantitative structure-activity relationship (QSAR) in cancer pharmacotherapy. This study utilized graph theoretical molecular descriptors, also known as topological indices, as a numerical representation method for the chemical structures embedded in molecular graphs. These descriptors, derived from molecular graphs, play a pivotal role in quantitative structure-property relationship (QSPR) analysis. In this paper, intercorrelation between the Balban index, connective eccentric index, eccentricity connectivity index, harmonic index, hyper Zagreb index, first path Zagreb index, second path Zagreb index, Randic index, sum connectivity index, graph energy and Laplacian energy is studied on the set of molecular graphs of coumarins. It is found that the pairs of degree-based indices are highly intercorrelated. The use of these molecular descriptors in structure-boiling point modeling was analyzed. Finally, the curve-linear regression between considered molecular descriptors with physicochemical properties of coumarins and coumarin-related compounds is obtained.

Computation-Driven Rational Design of Self-Assembled Short Peptides for Catalytic Hydrogen Production.

Self-assembling peptides represent a captivating area of study in nanotechnology and biomaterials. This interest is largely driven by their unique properties and the vast application potential across various fields such as catalytic functions. However, design complexities, including high-dimensional sequence space and structural diversity, pose significant challenges in the study of such systems. In this work, we explored the possibility of self-assembled peptides to catalyze the hydrolysis of hydrosilane for hydrogen production using ab initio calculations and carried out wet-lab experiments to confirm the feasibility of these catalytic reactions under ambient conditions. Further, we delved into the nuanced interplay between sequence, structural conformation, and catalytic activity by combining modeling with experimental techniques such as transmission electron microscopy and nuclear magnetic resonance and proposed a dual mode of the microstructure of the catalytic center. Our results reveal that although research in this area is still at an early stage, the development of self-assembled peptide catalysts for hydrogen production has the potential to provide a more sustainable and efficient alternative to conventional hydrogen production methods. In addition, this work also demonstrates that a computation-driven rational design supplemented by experimental validation is an effective protocol for conducting research on functional self-assembled peptides.

Advancing molecular modeling and reverse vaccinology in broad-spectrum yellow fever virus vaccine development.

Yellow fever outbreaks are prevalent, particularly in endemic regions. Given the lack of an established treatment for this disease, significant attention has been directed toward managing this arbovirus. In response, we developed a multiepitope vaccine designed to elicit an immune response, utilizing advanced immunoinformatic and molecular modeling techniques. To achieve this, we predicted B- and T-cell epitopes using the sequences from all structural (E, prM, and C) and nonstructural proteins of 196 YFV strains. Through comprehensive analysis, we identified 10 cytotoxic T-lymphocyte (CTL) and 5T-helper (Th) epitopes that exhibited overlap with B-lymphocyte epitopes. These epitopes were further evaluated for their affinity to a wide range of human leukocyte antigen system alleles and were rigorously tested for antigenicity, immunogenicity, allergenicity, toxicity, and conservation. These epitopes were linked to an adjuvant ( ß -defensin) and to each other using ligands, resulting in a vaccine sequence with appropriate physicochemical properties. The 3D structure of this sequence was created, improved, and quality checked; then it was anchored to the Toll-like receptor. Molecular Dynamics and Quantum Mechanics/Molecular Mechanics simulations were employed to enhance the accuracy of docking calculations, with the QM portion of the simulations carried out utilizing the density functional theory formalism. Moreover, the inoculation model was able to provide an optimal codon sequence that was inserted into the pET-28a( +) vector for in silico cloning and could even stimulate highly relevant humoral and cellular immunological responses. Overall, these results suggest that the designed multi-epitope vaccine can serve as prophylaxis against the yellow fever virus.

Intermolecular hydrogen bonding behavior of amino acid radical cations.

Amino acid and peptide radicals are of broad interest due to their roles in biochemical oxidative damage, pathogenesis and protein radical catalysis, among others. Using density functional theory (DFT) calculations at the ωB97X-D/def2-QZVPPD//ωB97X-D/def2-TZVPP level of theory, we systematically investigated the hydrogen bonding between water and fourteen α-amino acids (Ala, Asn, Cys, Gln, Gly, His, Met, Phe, Pro, Sel, Ser, Thr, Trp, and Tyr) in both neutral and radical cation forms. For all amino acids surveyed, stronger hydrogen-bonding interactions with water were observed upon single-electron oxidation, with the greatest increases in hydrogen-bonding strength occurring in Gly, Ala and His. We demonstrate that the side chain has a significant impact on the most favorable hydrogen-bonding modes experienced by amino acid radical cations. Our computations also explored the fragmentation of amino acid radical cations through the loss of a COOH radical facilitated by hydrogen bonding. The most favorable pathways provided stabilization of the resulting cationic fragments through hydrogen bonding, resulting in more favorable thermodynamics for the fragmentation process. These results indicate that non-covalent interactions with the environment have a profound impact on the structure and chemical fate of oxidized amino acids.

Affinity of Nicotinoids to a Model Nicotinic Acetylcholine Receptor (nAChR) Binding Pocket in the Human Brain.

The binding affinity of nicotinoids to the binding residues of the α4ß2 variant of the nicotinic acetylcholine receptor (nAChR) was identified as a strong predictor of the nicotinoid's addictive character. Using ab initio calculations for model binding pockets of increasing size composed of 3, 6, and 14 amino acids (3AA, 6AA, and 14AA) that are derived from the crystal structure, the differences in binding affinity of 6 nicotinoids, namely, nicotine (NIC), nornicotine (NOR), anabasine (ANB), anatabine (ANT), myosmine (MYO), and cotinine (COT) were correlated to their previously reported doses required for increases in intracranial self-stimulation (ICSS) thresholds, a metric for their addictive function. By employing the many-body decomposition, the differences in the binding affinities of the various nicotinoids could be attributed mainly to the proton exchange energy between the pyridine and non-pyridine rings of the nicotinoids and the interactions between them and a handful of proximal amino acids, namely Trp156, Trpß57, Tyr100, and Tyr204. Interactions between the guest nicotinoid and the amino acids of the binding pocket were found to be mainly classical in nature, except for those between the nicotinoid and Trp156. The larger pockets were found to model binding structures more accurately and predicted the addictive character of all nicotinoids, while smaller models, which are more computationally feasible, would only predict the addictive character of nicotinoids that are similar to nicotine. The present study identifies the binding affinity of the guest nicotinoid to the host binding pocket as a strong descriptor of the nicotinoid's addiction potential, and as such it can be employed as a fast-screening technique for the potential addiction of nicotine analogs.

Structural insights in to the atypical type-I ABC Glucose-6-phosphate importer VCA0625-27 of Vibrio cholerae.

Sugar phosphates are potential sources of carbon and phosphate for bacteria. Despite that the process of internalization of Glucose-6-Phosphate (G6P) through plasma membrane remained elusive in several bacteria. VCA0625-27, made of periplasmic ligand binding protein (PLBP) VCA0625, an atypical monomeric permease VCA0626, and a cytosolic ATPase VCA0627, recently emerged as hexose-6-phosphate uptake system of Vibrio cholerae. Here we report high resolution crystal structure of VCA0625 in G6P bound state that largely resembles AfuA of Actinobacillus pleuropneumoniae. MD simulations on VCA0625 in apo and G6P bound states unraveled an 'open to close' and swinging bi-lobal motions, which are diminished upon G6P binding. Mutagenesis followed by biochemical assays on VCA0625 underscored that R34 works as gateway to bind G6P. Although VCA0627 binds ATP, it is ATPase deficient in the absence of VCA0625 and VCA0626, which is a signature phenomenon of type-I ABC importer. Further, modeling, docking and systematic sequence analysis allowed us to envisage the existence of similar atypical type-I G6P importer with fused monomeric permease in 27 other gram-negative bacteria.

Conformational analysis and spectroscopic properties of antichagasic nifurtimox.

Considering the health relevance of Chagas' disease, recent research efforts have focused on developing more efficient drug delivery systems containing nifurtimox (NFX). This paper comprehensively investigates NFX through conformational analysis and spectroscopic characterization. Using a conformer-rotamer ensemble sampling tool (CREST-xtb), five distinct conformers of NFX were sampled within a 3.0 kcal mol-1 relative energy window. Subsequently, such structures were used as inputs for geometry optimization by density functional theory (DFT) at B3LYP-def2-TZVP level of theory. Notably, harmonic vibrational frequencies were calculated to establish an in-depth comparison with experimental results and existing literature for the NFX or similar molecules and functional groups, thereby achieving a widely reasoned assignment of the mid-infrared band absorptions for the first time. Moreover, UV-VIS spectra of NFX were obtained in several solvents, enabling the determination of the molar absorptivity coefficient for the two electronic transitions observed for NFX. Among the aprotic solvents, a bathochromic effect was observed in the function of the dielectric constants. Furthermore, a hypochromic effect was observed when the drug was dissolved in protic solvents. These findings offer crucial support for new drug delivery systems containing NFX while demonstrating the potential of spectrophotometric studies in establishing quality control assays for NFX drug products.

New tetrahydroisoquinoline-4-carbonitrile derivatives as potent agents against cyclin-dependent kinases, crystal structures, and computational studies.

The synthesis of two new hexahydroisoquinoline-4-carbonitrile derivatives (3a and 3b) is reported along with spectroscopic data and their crystal structures. In compound 3a, the intramolecular O-H···O hydrogen bond constraints the acetyl and hydroxyl groups to be syn. In the crystal, inversion dimers are generated by C-H···O hydrogen bonds and are connected into layers parallel to (10-1) by additional C-H···O hydrogen bonds. The layers are stacked with Cl···S contacts 0.17 Å less than the sum of the respective van der Waals radii. The conformation of the compound 3b is partially determined by the intramolecular O-H···O hydrogen bond. A puckering analysis of the tetrahydroisoquinoline unit was performed. In the crystal, O-H···O and C-H···O hydrogen bonds together with C-H···π(ring) interactions form layers parallel to (01-1) which pack with normal van der Waals interactions. To understand the binding efficiency and stability of the title molecules, molecular docking, and 100 ns dynamic simulation analyses were performed with CDK5A1. To rationalize their structure-activity relationship(s), a DFT study at the B3LYP/6-311++G** theoretical level was also done. The 3D Hirshfled surfaces were also taken to investigate the crystal packings of both compounds. In addition, their ADMET properties were explored.Communicated by Ramaswamy H. Sarma.