Identification of potent BRD4-BD1 inhibitors using classical and steered molecular dynamics based free energy analysis.

In the present work a combination of traditional and steered molecular dynamics based techniques were employed to identify potential inhibitors against the human BRD4 protein (BRD4- BD1); an established drug target for multiple illnesses including various malignancies. Quinoline derivatives that were synthesized in-house were tested for their potential as new BRD4-BD1 inhibitors. Initially molecular docking experiments were performed to determine the binding poses of BRD4-BD1 inhibitors. To learn more about the thermodynamics of inhibitor binding to the BRD4-BD1 active site, the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) free energy calculations were conducted afterwards. The findings of the MM-PBSA analysis were further reinforced by performing steered umbrella sampling simulations which revealed crucial details about the binding/unbinding process of the most potent quinoline derivatives at the BRD4-BD1 active site. We report a novel quinoline derivative which can be developed into a fully functional BRD4-BD1 inhibitor after experimental validation. The identified compound (4 g) shows better properties than the standard BRD4-BD1 inhibitors considered in the study. The study also highlights the crucial role of Gln78, Phe79, Trp81, Pro82, Phe83, Gln84, Gln85, Val87, Leu92, Leu94, Tyr97, Met105, Cys136, Asn140, Ile146 and Met149 in inhibitor binding. The study provides a possible lead candidate and key amino acids involved in inhibitor recognition and binding at the active site of BRD4-BD1 protein. The findings might be of significance to medicinal chemists involved in the development of potent BRD4-BD1 inhibitors.

Thermal isomerization rates in retinal analogues using Ab-Initio molecular dynamics.

For a detailed understanding of chemical processes in nature and industry, we need accurate models of chemical reactions in complex environments. While Eyring transition state theory is commonly used for modeling chemical reactions, it is most accurate for small molecules in the gas phase. A wide range of alternative rate theories exist that can better capture reactions involving complex molecules and environmental effects. However, they require that the chemical reaction is sampled by molecular dynamics simulations. This is a formidable challenge since the accessible simulation timescales are many orders of magnitude smaller than typical timescales of chemical reactions. To overcome these limitations, rare event methods involving enhanced molecular dynamics sampling are employed. In this work, thermal isomerization of retinal is studied using tight-binding density functional theory. Results from transition state theory are compared to those obtained from enhanced sampling. Rates obtained from dynamical reweighting using infrequent metadynamics simulations were in close agreement with those from transition state theory. Meanwhile, rates obtained from application of Kramers' rate equation to a sampled free energy profile along a torsional dihedral reaction coordinate were found to be up to three orders of magnitude higher. This discrepancy raises concerns about applying rate methods to one-dimensional reaction coordinates in chemical reactions.

Artificial Intelligence-Powered Molecular Docking and Steered Molecular Dynamics for Accurate scFv Selection of Anti-CD30 Chimeric Antigen Receptors.

Chimeric antigen receptor (CAR) T cells represent a revolutionary immunotherapy that allows specific tumor recognition by a unique single-chain fragment variable (scFv) derived from monoclonal antibodies (mAbs). scFv selection is consequently a fundamental step for CAR construction, to ensure accurate and effective CAR signaling toward tumor antigen binding. However, conventional in vitro and in vivo biological approaches to compare different scFv-derived CARs are expensive and labor-intensive. With the aim to predict the finest scFv binding before CAR-T cell engineering, we performed artificial intelligence (AI)-guided molecular docking and steered molecular dynamics analysis of different anti-CD30 mAb clones. Virtual computational scFv screening showed comparable results to surface plasmon resonance (SPR) and functional CAR-T cell in vitro and in vivo assays, respectively, in terms of binding capacity and anti-tumor efficacy. The proposed fast and low-cost in silico analysis has the potential to advance the development of novel CAR constructs, with a substantial impact on reducing time, costs, and the need for laboratory animal use.

Molecular dynamics study on micelle-small molecule interactions: developing a strategy for an extensive comparison.

Theoretical predictions of the solubilizing capacity of micelles and vesicles present in intestinal fluid are important for the development of new delivery techniques and bioavailability improvement. A balance between accuracy and computational cost is a key factor for an extensive study of numerous compounds in diverse environments. In this study, we aimed to determine an optimal molecular dynamics (MD) protocol to evaluate small-molecule interactions with micelles composed of bile salts and phospholipids. MD simulations were used to produce free energy profiles for three drug molecules (danazol, probucol, and prednisolone) and one surfactant molecule (sodium caprate) as a function of the distance from the colloid center of mass. To address the challenges associated with such tasks, we compared different simulation setups, including freely assembled colloids versus pre-organized spherical micelles, full free energy profiles versus only a few points of interest, and a coarse-grained model versus an all-atom model. Our findings demonstrate that combining these techniques is advantageous for achieving optimal performance and accuracy when evaluating the solubilization capacity of micelles. All-atom (AA) and coarse-grained (CG) umbrella sampling (US) simulations and point-wise free energy (FE) calculations were compared to their efficiency to computationally analyze the solubilization of active pharmaceutical ingredients in intestinal fluid colloids.

Studying specificity in protein-glycosaminoglycan recognition with umbrella sampling.

In the past few decades, glycosaminoglycan (GAG) research has been crucial for gaining insights into various physiological, pathological, and therapeutic aspects mediated by the direct interactions between the GAG molecules and diverse proteins. The structural and functional heterogeneities of GAGs as well as their ability to bind specific proteins are determined by the sugar composition of the GAG, the size of the GAG chains, and the degree and pattern of sulfation. A deep understanding of the interactions in protein-GAG complexes is essential to explain their biological functions. In this study, the umbrella sampling (US) approach is used to pull away a GAG ligand from the binding site and then pull it back in. We analyze the binding interactions between GAGs of three types (heparin, desulfated heparan sulfate, and chondroitin sulfate) with three different proteins (basic fibroblast growth factor, acidic fibroblast growth factor, and cathepsin K). The main focus of our study was to evaluate whether the US approach is able to reproduce experimentally obtained structures, and how useful it can be for getting a deeper understanding of GAG properties, especially protein recognition specificity and multipose binding. We found that the binding free energy landscape in the proximity of the GAG native binding pose is complex and implies the co-existence of several binding poses. The sliding of a GAG chain along a protein surface could be a potential mechanism of GAG particular sequence recognition by proteins.

Computational approaches to investigate fluoride binding, selectivity and transport across the membrane.

The use of molecular dynamics (MD) simulations to study biomolecular systems has proven reliable in elucidating atomic-level details of structure and function. In this chapter, MD simulations were used to uncover new insights into two phylogenetically unrelated bacterial fluoride (F-) exporters: the CLCF F-/H+ antiporter and the Fluc F- channel. The CLCF antiporter, a member of the broader CLC family, has previously revealed unique stoichiometry, anion-coordinating residues, and the absence of an internal glutamate crucial for proton import in the CLCs. Through MD simulations enhanced with umbrella sampling, we provide insights into the energetics and mechanism of the CLCF transport process, including its selectivity for F- over HF. In contrast, the Fluc F- channel presents a novel architecture as a dual topology dimer, featuring two pores for F- export and a central non-transported sodium ion. Using computational electrophysiology, we simulate the electrochemical gradient necessary for F- export in Fluc and reveal details about the coordination and hydration of both F- and the central sodium ion. The procedures described here delineate the specifics of these advanced techniques and can also be adapted to investigate other membrane protein systems.

A molecular arm: the molecular bending-unbending mechanism of integrin.

The balance of integrin activation and deactivation regulates its function and mediates cell behaviors. Mechanical force triggers the unbending and activation of integrin. However, how an activated and extended integrin spontaneously bends back is unclear. I performed all-atom molecular dynamics simulations on an integrin or its subunits to reveal the bending-unbending mechanism of integrin. According to the simulations, the integrin structure works like a human arm. The integrin α subunit serves as the bones, while the ß leg serves as the bicep. The integrin extension results in the stretching of the ß leg, and the extended integrin spontaneously bends as a consequence of the contraction of the ß leg. This study provides new insights into the mechanism of how the integrin secures in the bent inactivated state and sheds light on how the integrin could achieve a stable extended state.

Graphene Exfoliation in Binary NMP/Water Mixtures by Molecular Dynamics Simulations.

We investigate the molecular mechanism underlying the liquid-phase exfoliation of graphene in aqueous/N-methyl-2-pyrrolidone (NMP) solvent mixtures and calculate the associated free energies, considering different NMP concentrations and exfoliation temperatures. We employ steered molecular dynamics to establish a path for the exfoliation of a graphene sheet from graphite within each solvent environment. Then, we conduct umbrella sampling simulations throughout the created paths to compute the potential of mean force (PMF) of the graphene sheet. As the exfoliated nanosheet disperses into the liquid, it becomes fully covered by an adsorbed solvent monolayer. We analyze the composition of the monolayer by measuring the direct contacts of either NMP or water molecules with the carbon surface. The carbon surface exhibits a preference for adsorbing NMP over water. The NMP molecules form a hydrophobic compact monolayer structure, effectively protecting the carbon interface from unfavorable interactions with water. The creation of the hydrophobic monolayer is a key factor in the exfoliation process, as it effectively inhibits the restacking of exfoliated nanosheets. An adequate level of graphene solubility is achieved through the addition of 20 % to 30 % water by weight to the NMP solvent. This finding holds significant importance for improving production efficiency and reducing dependence on organic solvents in the industrial manufacturing of graphene.

Structure-based virtual screening of novel USP5 inhibitors targeting the zinc finger ubiquitin-binding domain.

The equilibrium of cellular protein levels is pivotal for maintaining normal physiological functions. USP5 belongs to the deubiquitination enzyme (DUBs) family, controlling protein degradation and preserving cellular protein homeostasis. Aberrant expression of USP5 is implicated in a variety of diseases, including cancer, neurodegenerative diseases, and inflammatory diseases. In this paper, a multi-level virtual screening (VS) approach was employed to target the zinc finger ubiquitin-binding domain (ZnF-UBD) of USP5, leading to the identification of a highly promising candidate compound 0456-0049. Molecular dynamics (MD) simulations were then employed to assess the stability of complex binding and predict hotspot residues in interactions. The results indicated that the candidate stably binds to the ZnF-UBD of USP5 through crucial interactions with residues ARG221, TRP209, GLY220, ASN207, TYR261, TYR259, and MET266. Binding free energy calculations, along with umbrella sampling (US) simulations, underscored a superior binding affinity of the candidate relative to known inhibitors. Moreover, US simulations revealed conformational changes of USP5 during ligand dissociation. These insights provide a valuable foundation for the development of novel inhibitors targeting USP5.

Dissecting the role of ALK double mutations in drug resistance to lorlatinib with in-depth theoretical modeling and analysis.

Anaplastic lymphoma kinase (ALK) is implicated in the genesis of multiple malignant tumors. Lorlatinib stands out as the most advanced and effective inhibitor currently used in the clinic for the treatment of ALK-positive non-small cell lung cancer. However, resistance to lorlatinib has inevitably manifested over time, with double/triple mutations of G1202, L1196, L1198, C1156 and I1171 frequently observed in clinical practice, and tumors regrow within a short time after treatment with lorlatinib. Therefore, elucidating the mechanism of resistance to lorlatinib is paramount in paving the way for innovative therapeutic strategies and the development of next-generation drugs. In this study, we leveraged multiple computational methodologies to delve into the resistance mechanisms of three specific double mutations of ALKG1202R/L1196M, ALKG1202R/L1198F and ALKI1171N/L1198F to lorlatinib. We analyzed these mechanisms through qualitative (PCA, DCCM) and quantitative (MM/GBSA, US) kinetic analyses. The qualitative analysis shows that these mutations exert minimal perturbations on the conformational dynamics of the structural domains of ALK. The energetic and structural assessments show that the van der Waals interactions, formed by the conserved residue Leu1256 within the ATP-binding site and the residues Glu1197 and Met1199 in the hinge domain with lorlatinib, play integral roles in the occurrence of drug resistance. Furthermore, the US simulation results elucidate that the pathways through which lorlatinib dissociates vary across mutant systems, and the distinct environments during the dissociation process culminate in diverse resistance mechanisms. Collectively, these insights provide important clues for the design of novel inhibitors to combat resistance.

Mechanistic insights into the conversion of flavin adenine dinucleotide (FAD) to 8-formyl FAD in formate oxidase: a combined experimental and in-silico study.

Formate oxidase (FOx), which contains 8-formyl flavin adenine dinucleotide (FAD), exhibits a distinct advantage in utilizing ambient oxygen molecules for the oxidation of formic acid compared to other glucose-methanol-choline (GMC) oxidoreductase enzymes that contain only the standard FAD cofactor. The FOx-mediated conversion of FAD to 8-formyl FAD results in an approximate 10-fold increase in formate oxidase activity. However, the mechanistic details underlying the autocatalytic formation of 8-formyl FAD are still not well understood, which impedes further utilization of FOx. In this study, we employ molecular dynamics simulation, QM/MM umbrella sampling simulation, enzyme activity assay, site-directed mutagenesis, and spectroscopic analysis to elucidate the oxidation mechanism of FAD to 8-formyl FAD. Our results reveal that a catalytic water molecule, rather than any catalytic amino acids, serves as a general base to deprotonate the C8 methyl group on FAD, thus facilitating the formation of a quinone-methide tautomer intermediate. An oxygen molecule subsequently oxidizes this intermediate, resulting in a C8 methyl hydroperoxide anion that is protonated and dissociated to form OHC-RP and OH-. During the oxidation of FAD to 8-formyl FAD, the energy barrier for the rate-limiting step is calculated to be 22.8 kcal/mol, which corresponds to the required 14-hour transformation time observed experimentally. Further, the elucidated oxidation mechanism reveals that the autocatalytic formation of 8-formyl FAD depends on the proximal arginine and serine residues, R87 and S94, respectively. Enzymatic activity assay validates that the mutation of R87 to lysine reduces the kcat value to 75% of the wild-type, while the mutation to histidine results in a complete loss of activity. Similarly, the mutant S94I also leads to the deactivation of enzyme. This dependency arises because the nucleophilic OH- group and the quinone-methide tautomer intermediate are stabilized through the noncovalent interaction provided by R87 and S94. These findings not only explain the mechanistic details of each reaction step but also clarify the functional role of R87 and S94 during the oxidative maturation of 8-formyl FAD, thereby providing crucial theoretical support for the development of novel flavoenzymes with enhanced redox properties.

Driving forces and large scale affinity calculations for piperine/γ-cyclodxetrin complexes: Mechanistic insights from umbrella sampling simulation and DFT calculations.

Piperine (PiP), a bioactive molecule, exhibits numerous health benefits and is frequently employed as a co-delivery agent with various phytomedicines (e.g., curcumin) to enhance their bioavailability. This is attributed to PiP's inhibitory activity against drug-metabolizing proteins, notably CYP3A4. Nevertheless, PiP encounters solubility challenges addressed in this study using cyclodextrins (CDs). Specifically, γ-CD and its derivatives, Hydroxypropyl-γ-CD (HP-γ-CD), and Octakis (6-O-sulfo)-γ-CD (Octakis-S-γ-CD), were employed to form supramolecular complexes with PiP. The conformational space of the complexes was assessed through 1 µs molecular dynamics simulations and umbrella sampling. Additionally, quantum mechanical calculations using wB97X-D dispersion-corrected DFT functional and 6-311 + G(d,p) basis set were conducted on the complexes to examine the thermodynamics and kinetic stability. Results indicated that Octakis-S-γ-CD exhibits superior host capabilities for PiP, with the most favorable complexation energy (-457.05 kJ/mol), followed by HP-γ-CD (-249.16 kJ/mol). Furthermore, two conformations of the Octakis-S-γ-CD/PiP complex were explored to elucidate the optimal binding orientation of PiP within the binding pocket of Octakis-S-γ-CD. Supramolecular chemistry relies significantly on non-covalent interactions. Therefore, our investigation extensively explores the critical atoms involved in these interactions, elucidating the influence of substituted groups on the stability of inclusion complexes. This comprehensive analysis contributes to emphasizing the γ-CD derivatives with improved host capacity.

The mechanism and energetics of the dynein priming stroke.

Dyneins are an AAA+ motor responsible for motility and force generation toward the minus end of microtubules. Dynein motility is powered by nucleotide-dependent transitions of its linker domain, which transitions between straight (post-powerstroke) and bent (pre-powerstroke) conformations. To understand the dynamics and energetics of the linker, we performed all-atom molecular dynamics simulations of human dynein-2 primed for its power stroke. Simulations revealed that the linker can adopt either a bent conformation or a semi-bent conformation, separated by a 5.7 kT energy barrier. The linker cannot switch back to its straight conformation in the pre-powerstroke state due to a steric clash with the AAA+ ring. Simulations also showed that an isolated linker has a free energy minimum near the semi-bent conformation in the absence of the AAA+ ring, indicating that the linker stores energy as it bends and releases this energy during the powerstroke.

Imidazole[1,5-a]pyridine derivatives as EGFR tyrosine kinase inhibitors unraveled by umbrella sampling and steered molecular dynamics simulations.

Although the use of the tyrosine kinase inhibitors (TKIs) has been proved that it can save live in a cancer treatment, the currently used drugs bring in many undesirable side-effects. Therefore, the search for new drugs and an evaluation of their efficiency are intensively carried out. Recently, a series of eighteen imidazole[1,5-a]pyridine derivatives were synthetized by us, and preliminary analyses pointed out their potential to be an important platform for pharmaceutical development owing to their promising actions as anticancer agents and enzyme (kinase, HIV-protease,…) inhibitors. In the present theoretical study, we further analyzed their efficiency in using a realistic scenario of computational drug design. Our protocol has been developed to not only observe the atomistic interaction between the EGFR protein and our 18 novel compounds using both umbrella sampling and steered molecular dynamics simulations, but also determine their absolute binding free energies. Calculated properties of the 18 novel compounds were in detail compared with those of two known drugs, erlotinib and osimertinib, currently used in cancer treatment. Inspiringly the simulation results promote three imidazole[1,5-a]pyridine derivatives as promising inhibitors into a further step of clinical trials.