Rings, Hexagons, Hetals, and Dipolar Moment Sink-Sources: The Fanciful Behavior of Water around Cyclodextrin Complexes

Cyclodextrins: Establishing building blocks for AI-driven drug design by determining affinity constants in silico

Cyclodextrins (CDs) are cyclic carbohydrate polymers that hold significant promise for drug delivery and industrial applications. Their effectiveness depends on their ability to encapsulate target molecules with strong affinity and specificity, but quantifying affinities in these systems accurately is challenging for a variety of reasons. Computational methods represent an exceptional complement to in vitro assays because they can be employed for existing and hypothetical molecules, providing high resolution structures in addition to a mechanistic, dynamic, kinetic, and thermodynamic characterization. Here, we employ potential of mean force (PMF) calculations obtained from guided metadynamics simulations to characterize the 1:1 inclusion complexes between four different modified βCDs, with different type, number, and location of substitutions, and two sterol molecules (cholesterol and 7-ketocholesterol). Our methods, validated for reproducibility through four independent repeated simulations per system and different post processing techniques, offer new insights into the formation and stability of CD-sterol inclusion complexes. A systematic distinct orientation preference where the sterol tail projects from the CD's larger face and significant impacts of CD substitutions on binding are observed. Notably, sampling only the CD cavity's wide face during simulations yielded comparable binding energies to full-cavity sampling, but in less time and with reduced statistical uncertainty, suggesting a more efficient approach. Bridging computational methods with complex molecular interactions, our research enables predictive CD designs for diverse applications. Moreover, the high reproducibility, sensitivity, and cost-effectiveness of the studied methods pave the way for extensive studies of massive CD-ligand combinations, enabling AI algorithm training and automated molecular design.

Simulating Bacterial Membrane Models at the Atomistic Level: A Force Field Comparison

Molecular dynamics (MD) simulations are currently an indispensable tool to understand both the dynamic and nanoscale organization of cell membrane models. A large number of quantitative parameters can be extracted from these simulations, but their reliability is determined by the quality of the employed force field and the simulation parameters. Much of the work on parametrizing and optimizing force fields for biomembrane modeling has been focused on homogeneous bilayers with a single phospholipid type. However, these may not perform effectively or could even be unsuitable for lipid mixtures commonly employed in membrane models. This work aims to fill this gap by comparing MD simulation results of several bacterial membrane models using different force fields and simulation parameters, namely, CHARMM36, Slipids, and GROMOS-CKP. Furthermore, the hydrogen isotope exchange (HIE) method, combined with GROMOS-CKP (GROMOS-H2Q), was also tested to check for the impact of this acceleration strategy on the performance of the force field. A common set of simulation parameters was employed for all of the force fields in addition to those corresponding to the original parametrization of each of them. Furthermore, new experimental order parameter values determined from NMR of several lipid mixtures are also reported to compare them with those determined from MD simulations. Our results reveal that most of the calculated physical properties of bacterial membrane models from MD simulations are substantially force field and lipid composition dependent. Some lipid mixtures exhibit nearly ideal behaviors, while the interaction of different lipid types in other mixtures is highly synergistic. None of the employed force fields seem to be clearly superior to the other three, each having its own strengths and weaknesses. Slipids are notably effective at replicating the order parameters for all acyl chains, including those in lipid mixtures, but they offer the least accurate results for headgroup parameters. Conversely, CHARMM provides almost perfect estimates for the order parameters of the headgroups but tends to overestimate those of the lipid tails. The GROMOS parametrizations deliver reasonable order parameters for entire lipid molecules, including multicomponent bilayers, although they do not reach the accuracy of Slipids for tails or CHARMM for headgroups. Importantly, GROMOS-H2Q stands out for its computational efficiency, being at least 3 times faster than GROMOS, which is already faster than both CHARMM and Slipids. In turn, GROMOS-H2Q yields much higher compressibilities compared to all other parametrizations.

Unraveling the molecular dynamics of sugammadex-rocuronium complexation: A blueprint for cyclodextrin drug design

Sugammadex, marketed as Bridion™, is an approved cyclodextrin (CD) based drug for the reversal of neuromuscular blockade in adults undergoing surgery. Sugammadex forms an inclusion complex with the neuromuscular blocking agent (NMBA) rocuronium, allowing rapid reversal of muscle paralysis. In silico methods have been developed for studying CD inclusion complexes, aimed at accurately predicting their structural, energetic, dynamic, and kinetic properties, as well as binding constants. Here, a computational study aimed at characterizing the sugammadex-rocuronium system from the perspective of docking calculations, free molecular dynamics (MD) simulations, and biased metadynamics simulations with potential of mean force (PMF) calculations is presented. The aim is to provide detailed information about this system, as well as to use it as a model system for validation of the methods. This method predicts results in line with experimental evidence for both the optimal structure and the quantitative value for the binding constant. Interestingly, there is a less profound preference for the orientation than might be assumed based on electrostatic interactions, suggesting that both orientations may exist in solution. These results show that this technology can efficiently analyze CD inclusion complexes and could be used to facilitate the development and optimization of novel applications for CDs.

Effect of temperature on the structure and drug-release behaviour of inclusion complex of β-cyclodextrin with cyclophosphamide: a molecular dynamics study

Cyclodextrins (CDs) are suitable drug carriers because of their doughnut-shaped cavities with hydrophilic outer and hydrophobic inner surfaces. Temperature-responsive CD-based drug carriers are expected to be one of the most promising candidates for drug delivery systems. In this study, we performed molecular dynamics simulations of the inclusion complex of β-CD with cyclophosphamide (CP) at temperatures from 300 K to 400 K to investigate the temperature dependency of the release behaviour of CP and structural changes of β-CD in an aqueous solution. We analysed the distance between the centres of mass of β-CD and CP and the radius of gyration of β-CD. The CP molecule was released from the β-CD cavity at 400 K, whereas two different inclusion complexes, partially and completely, were observed at T < 400 K. β-CD encapsulating a CP molecule had a more spherical shape and rigidity than β-CD without a CP, and the rigidity of their inclusion complex decreased with increasing temperature. Our findings provide fundamental insights into the behaviours of the β-CD/CP complex and drug release at the molecular level and can facilitate the development of new temperature-responsive drug delivery systems with CD nanocarriers triggered by localised temperature increases using focused ultrasound.

Cyclodextrin Inclusion Complexes and Their Application in Food Safety Analysis: Recent Developments and Future Prospects

Food safety issues are a major threat to public health and have attracted much attention. Therefore, exploring accurate, efficient, sensitive, and economical detection methods is necessary to ensure consumers' health. In this regard, cyclodextrins (CDs) are promising candidates because they are nontoxic and noncaloric. The main body of CDs is a ring structure with hydrophobic cavity and hydrophilic exterior wall. Due to the above characteristics, CDs can encapsulate small guest molecules into their cavities, enhance their stability, avoid agglomeration and oxidation, and, at the same time, interact through hydrogen bonding and electrostatic interactions. Additionally, they can selectively capture the target molecules to be detected and improve the sensitivity of food detection. This review highlights recent advances in CD inclusion technology in food safety analysis, covering various applications from small molecule and heavy metal sensing to amino acid and microbial sensing. Finally, challenges and prospects for CDs and their derivatives are presented. The current review can provide a reference and guidance for current research on CDs in the food industry and may inspire breakthroughs in this field.

Comparative bioavailability study following a single dose intravenous and buccal administration of remdesivir in rabbits

As remdesivir, the first FDA-approved drug for SARS-CoV-2 infection, can be used only for hospitalized patients due to intravenous administration, there is an urgent need of effective oral antiviral formulations to be used at early stage of infection in an outpatient setting. The present paper reports on the comparative pharmacokinetics of the electrospun nanofiber remdesivir/sulfobutyl ether beta-cyclodextrin formulation after intravenous and buccal administration. It was postulated that oral transmucosal administration avoids remdesivir from metabolic transformation and intact remdesivir can be detected in plasma, but only the active metabolite GS-441524 could be experimentally detected at a significantly lower plasma level, than that provided by the intravenous route. In buccally treated animals, the metabolite GS-441524 appeared only at 1 h after treatment, while in intravenously treated animals, GS-441524 was possible to quantify even at the first time-point of blood collection. Further optimization of formulation is required to improve pharmacokinetics of remdesivir-sulfobutyl ether beta-cyclodextrin formulation upon buccal administration.

Molecular Dynamics Simulations of Transmembrane Cyclic Peptide Nanotubes Using Classical Force Fields, Hydrogen Mass Repartitioning, and Hydrogen Isotope Exchange Methods: A Critical Comparison

Self-assembled cyclic peptide nanotubes with alternating D- and L-amino acid residues in the sequence of each subunit have attracted a great deal of attention due to their potential for new nanotechnology and biomedical applications, mainly in the field of antimicrobial peptides. Molecular dynamics simulations can be used to characterize these systems with atomic resolution at different time scales, providing information that is difficult to obtain via wet lab experiments. However, the performance of classical force fields typically employed in the simulation of biomolecules has not yet been extensively tested with this kind of highly constrained peptide. Four different classical force fields (AMBER, CHARMM, OPLS, and GROMOS), using a nanotube formed by eight D,L-α-cyclic peptides inserted into a lipid bilayer as a model system, were employed here to fill this gap. Significant differences in the pseudo-cylindrical cavities formed by the nanotubes were observed, the most important being the diameter of the nanopores, the number and location of confined water molecules, and the density distribution of the solvent molecules. Furthermore, several modifications were performed on GROMOS54a7, aiming to explore acceleration strategies of the MD simulations. The hydrogen mass repartitioning (HMR) and hydrogen isotope exchange (HIE) methods were tested to slow down the fastest degrees of freedom. These approaches allowed a significant increase in the time step employed in the equation of the motion integration algorithm, from 2 fs up to 5-7 fs, with no serious changes in the structural and dynamical properties of the nanopores. Subtle differences with respect to the simulations with the unmodified force fields were observed in the concerted movements of the cyclic peptides, as well as in the lifetime of several H-bonds. All together, these results are expected to contribute to better understanding of the behavior of self-assembled cyclic peptide nanotubes, as well as to support the methods tested to speed up general MD simulations; additionally, they do provide a number of quantitative descriptors that are expected to be used as a reference to design new experiments intended to validate and complement computational studies of antimicrobial cyclic peptides.

Aggregation versus inclusion complexes to solubilize drugs with cyclodextrins. A case study using sulphobutylether-β-cyclodextrins and remdesivir

The formation of small hybrid aggregates between excipient and drug molecules is one of the mechanisms that contributes to the solubilization of active principles in pharmaceutical formulations. The characterization of the formation, governing interactions and structure of such entities is not trivial since they are highly flexible and dynamic, quickly exchanging molecules from one to another. In the case of cyclodextrins, this mechanism and the formation of inclusion complexes synergistically cooperate to favour the bioavailability of drugs. In a previous study we reported a detailed characterization of the possible formation of inclusion complexes with 1:1 stoichiometry between remdesivir, the only antiviral medication currently approved by the United States Food and Drug Administration for treating COVID-19, and sulphobutylether-β-cyclodextrins. Here we extend our study to assess the role of the spontaneous aggregation in the solubilization of the same drug, by molecular dynamics simulations at different relative concentrations of both compounds. The number of sulphobutylether substitutions in the cyclodextrin structure and two different protonation states of the remdesivir molecule are considered. We aim to shed light in the solubilization mechanism of sulphobutylether-β-cyclodextrins, broadly used as an excipient in many pharmaceutical formulations, in particular in the case of remdesivir as an active compound.

Cyclodextrin dimers: A versatile approach to optimizing encapsulation and their application to therapeutic extraction of toxic oxysterols

We have developed a novel class of specifically engineered, dimerized cyclodextrin (CD) nanostructures for the encapsulation of toxic biomolecules such as 7-ketocholesterol (7KC). 7KC accumulates over time and causes dysfunction in many cell types, linking it to several age-related diseases including atherosclerosis and age-related macular degeneration (AMD). Presently, treatments for these diseases are invasive, expensive, and show limited benefits. CDs are cyclic glucose oligomers utilized to capture small, hydrophobic molecules. Here, a combination of in silico, in vitro, and ex vivo methods is used to implement a synergistic rational drug design strategy for developing CDs to remove atherogenic 7KC from cells and tissues. Mechanisms by which CDs encapsulate sterols are discussed, and we conclude that covalently linked head-to-head dimers of βCDs have substantially improved affinity for 7KC compared to monomers. We find that inclusion complexes can be stabilized or destabilized in ways that allow the design of CD dimers with increased 7KC selectivity while maintaining an excellent safety profile. These CD dimers are being developed as therapeutics to treat atherosclerosis and other debilitating diseases of aging.