Calculating potentials of mean force from steered molecular dynamics simulations

Specific interactions of peripheral membrane proteins with lipids: what can molecular simulations show us?

Peripheral membrane proteins (PMPs) can reversibly and specifically bind to biological membranes to carry out functions such as cell signalling, enzymatic activity, or membrane remodelling. Structures of these proteins and of their lipid-binding domains are typically solved in a soluble form, sometimes with a lipid or lipid headgroup at the binding site. To provide a detailed molecular view of PMP interactions with the membrane, computational methods such as molecular dynamics (MD) simulations can be applied. Here, we outline recent attempts to characterise these binding interactions, focusing on both intracellular proteins, such as phosphatidylinositol phosphate (PIP)-binding domains, and extracellular proteins such as glycolipid-binding bacterial exotoxins. We compare methods used to identify and analyse lipid-binding sites from simulation data and highlight recent work characterising the energetics of these interactions using free energy calculations. We describe how improvements in methodologies and computing power will help MD simulations to continue to contribute to this field in the future.

Ultrahigh strength, modulus, and conductivity of graphitic fibers by macromolecular coalescence

Theoretical considerations suggest that the strength of carbon nanotube (CNT) fibers be exceptional; however, their mechanical performance values are much lower than the theoretical values. To achieve macroscopic fibers with ultrahigh performance, we developed a method to form multidimensional nanostructures by coalescence of individual nanotubes. The highly aligned wet-spun fibers of single- or double-walled nanotube bundles were graphitized to induce nanotube collapse and multi-inner walled structures. These advanced nanostructures formed a network of interconnected, close-packed graphitic domains. Their near-perfect alignment and high longitudinal crystallinity that increased the shear strength between CNTs while retaining notable flexibility. The resulting fibers have an exceptional combination of high tensile strength (6.57 GPa), modulus (629 GPa), thermal conductivity (482 W/m·K), and electrical conductivity (2.2 MS/m), thereby overcoming the limits associated with conventional synthetic fibers.

Picking on Carbonate: Kinetic Selectivity in the Encapsulation of Anions

Supramolecular hosts bind to inorganic anions at a fast rate and select them in proportion with thermodynamic stability of the corresponding [anion⊂host] complexes, forming in a reversible manner. In this study, we describe the action of hexapodal capsule 1 and its remarkable ability to select anions based on a large span of rates by which they enter this host. The thermodynamic affinity of 1 toward eighteen anions extends over eight orders of magnitude (0<K_a <10⁸  M^-1 ; ¹ H NMR spectroscopy). The capsule would retain CO₃ ^2- (K_a =10⁷  M^-1 ) for hours in the presence of eleven competing anions, including stronger binding SO₄ ^2- , HAsO₄ ^2- and HPO₄ ^2- (K_a =10⁷ -10⁸  M^-1 ). The observed selection resulted from 1 possessing narrow apertures (ca. 3×6 Å) comparable in size to anions (d=3.5-7.1 Å) slowing down the encapsulation to last from seconds to days. The unorthodox mode of action of 1 sets the stage for creating hosts that pick anions by their ability to access the host.

Reconstruction of the Free Energy Profile for SUMO1 from Nonequilibrium Single-Molecule Pulling Experiments

Free energy profiles form the cornerstone in the study of protein folding and function. In this study, the free energy profile of SUMO1 protein is directly reconstructed using an extension of the Jarzynski equality from atomic force microscope (AFM) based single-molecule force spectroscopy (SMFS) experiments. SUMO1 is a ubiquitin-like posttranslational modifier protein having a β clamp motif in its structure, imparting it with mechanical stability. We use the Jarzynski equality to obtain the equilibrium free energy profile from repeated nonequilibrium single-molecule pulling experiments. Indeed, the free energy values determined by the Jarzynski equality are lesser than the normal work average at all extensions. The free energy profiles constructed for the two velocities (100 and 400 nm/s) overlap with each other. The unfolding free energy barrier is estimated to be ∼7.5 kcal/mol. We anticipate that the Jarzynski equality can be applied in a similar manner to other ubiquitin-like proteins to extract their differences in the free energy profile, and hence, the effect of sequence diversity of structurally homologous proteins on the free energy landscape can be studied.

Mechanical force can enhance c-Src kinase activity by impairing autoinhibition

Cellular mechanosensing is pivotal for virtually all biological processes, and many molecular mechano-sensors and their way of function are being uncovered. In this work, we suggest that c-Src kinase acts as a direct mechano-sensor. c-Src is responsible for, among others, cell proliferation, and shows increased activity in stretched cells. In its native state, c-Src has little basal activity, because its kinase domain binds to an SH2 and SH3 domain. However, it is known that c-Src can bind to p130Cas, through which force can be transmitted to the membrane. Using molecular dynamics simulations, we show that force acting between the membrane-bound N-terminus of the SH3 domain and p130Cas induces partial SH3 unfolding, thereby impeding rebinding of the kinase domain onto SH2/SH3 and effectively enhancing kinase activity. Forces involved in this process are slightly lower or similar to the forces required to pull out c-Src from the membrane through the myristoyl linker, and key interactions involved in this anchoring are salt bridges between negative lipids and nearby basic residues in c-Src. Thus, c-Src appears to be a candidate for an intriguing mechanosensing mechanism of impaired kinase inhibition, which can be potentially tuned by membrane composition and other environmental factors.

Molecular simulations of proteins: From simplified physical interactions to complex biological phenomena

Molecular dynamics simulation is the most popular computational technique for investigating the structural and dynamical behaviour of proteins, in search of the molecular basis of their function. Far from being a completely settled field of research, simulations are still evolving to best capture the essential features of the atomic interactions that govern a protein's inner motions. Modern force fields are becoming increasingly accurate in providing a physical description adequate to this purpose, and allow us to model complex biological systems under fairly realistic conditions. Furthermore, the use of accelerated sampling techniques is improving our access to the observation of progressively larger molecular structures, longer time scales, and more hidden functional events. In this review, the basic principles of molecular dynamics simulations and a number of key applications in the area of protein science are summarized, and some of the most important results are discussed. Examples include the study of the structure, dynamics and binding properties of 'difficult' targets, such as intrinsically disordered proteins and membrane receptors, and the investigation of challenging phenomena like hydration-driven processes and protein aggregation. The findings described provide an overall picture of the current state of this research field, and indicate new perspectives on the road ahead to the upcoming future of molecular simulations.

Structural analysis of P450 AmphL from Streptomyces nodosus provides insights into substrate selectivity of polyene macrolide antibiotic biosynthetic P450s

AmphL is a cytochrome P450 enzyme that catalyzes the C8 oxidation of 8-deoxyamphotericin B to the polyene macrolide antibiotic, amphotericin B. To understand this substrate selectivity, we solved the crystal structure of AmphL to a resolution of 2.0 Å in complex with amphotericin B and performed molecular dynamics (MD) simulations. A detailed comparison with the closely related P450, PimD, which catalyzes the epoxidation of 4,5-desepoxypimaricin to the macrolide antibiotic, pimaricin, reveals key catalytic structural features responsible for stereo- and regio-selective oxidation. Both P450s have a similar access channel that runs parallel to the active site I helix over the surface of the heme. Molecular dynamics simulations of substrate binding reveal PimD can “pull” substrates further into the P450 access channel owing to additional electrostatic interactions between the protein and the carboxyl group attached to the hemiketal ring of 4,5-desepoxypimaricin. This substrate interaction is absent in AmphL although the additional substrate -OH groups in 8-deoxyamphotericin B help to correctly position the substrate for C8 oxidation. Simulations of the oxy-complex indicates that these -OH groups may also participate in a proton relay network required for O₂ activation as has been suggested for two other macrolide P450s, PimD and P450eryF. These findings provide experimentally testable models that can potentially contribute to a new generation of novel macrolide antibiotics with enhanced antifungal and/or antiprotozoal efficacy.

Thermodynamics of ion exchange coupled with swelling reactions in hydrated clay minerals

Crystalline hydrates of swelling clay minerals (smectites) exhibit a strong coupling between their ion exchange and hydration/dehydration reactions. The uptake or removal of water from smectite interlayers as a result of a change in the environmental conditions also leads to the partitioning of cations. Three factors, the solid ion composition, the solid basal spacing/water content, and the aqueous solution composition, are all implicated in controlling the thermodynamics of ion exchange. However, conventional approaches to measuring the exchange free energy cannot separate the influence of each of these individual factors. Here, we explore the energetics of the swelling and ion exchange reactions in montmorillonite using a potential of mean force approach and the thermodynamic integration method within molecular simulations. We investigate the influence of solution and clay composition on the spontaneity of the reactions, focusing on the 2 water-layer hydration state. The swelling simulations provide the equilibrium water content, interlayer water structure, and basal spacings, while thermodynamic integration of sodium-potassium exchange in the aqueous solution and solid phase are combined to calculate ion exchange free energies as a function of solution composition. Results confirm the tendency of the clay to collapse to lower hydration states as the concentration of the solution increases. Changes to the equilibrium water content, even at fixed hydration states, and the composition of the mixed electrolyte solution play a critical role in driving ion exchange and the selectivities of the clay to the exchanged cation, while the composition of the solid phase is shown to be insignificant. These findings underscore the extreme sensitivity of clay swelling and ion exchange thermodynamics to small (tenths of an Angstrom) deviations in layer spacing.

A modified Jarzynski free-energy estimator to eliminate non-conservative forces and its application in nanoparticle-membrane interactions

Computational methods to understand interactions in bio-complex systems are however limited to time-scales typically much shorter than in Nature. For example, on the nanoscale level, interactions between nanoparticles (NPs)/molecules/peptides and membranes are central in complex biomolecular processes such as membrane-coated NPs or cellular uptake. This can be remedied by the application of e.g. Jarzynski's equality where thermodynamic properties are extracted from non-equilibrium simulations. Although, the out of equilibrium work leads to non-conservative forces. We here propose a correction Pair Forces method, that removes these forces. Our proposed method is based on the calculation of pulling forces in backward and forward directions for the Jarzynski free-energy estimator using steered molecular dynamics simulation. Our results show that this leads to much improvement for NP-membrane translocation free energies. Although here we have demonstrated the application of the method in molecular dynamics simulation, it could be applied for experimental approaches.

Free Energy Profiles of Proton Transfer Reactions: Density Functional Benchmark from Biased Ab Initio Dynamics

By coupling an enhanced sampling algorithm with an orbital-localized variant of Car-Parrinello molecular dynamics, the so-called atomic centered density matrix propagation model, we reconstruct the free energy profiles along reaction pathways using different density functional approximations (DFAs) ranging from locals to hybrids. In particular, we compare the computed free energy barrier height of proton transfer (PT) reactions to those obtained by a more traditional static approach, based on the intrinsic reaction coordinate (IRC), for two case systems, namely malonaldehyde and formic acid dimer. The obtained results show that both the IRC profiles and the potentials of mean force, derived from biased dynamic trajectories, are very sensitive to the density functional approximation applied. More precisely, we observe that, with the notable exception of M06-L, local density functionals always strongly underestimate the reaction barrier heights. More generally, we find that also the shape of the free energy profile is very sensitive to the density functional choice, thus highlighting the effect, often neglected, that the choice of DFA has also in the case of dynamics simulations.

Molecular dynamics study of the internalization of cell-penetrating peptides containing unnatural amino acids across membranes

Peptide-based delivery systems that deliver target molecules into cells have been gaining traction. These systems need cell-penetrating peptides (CPPs), which have the remarkable ability to penetrate into biological membranes and help internalize different cargoes into cells through the cell membranes. The molecular internalization mechanism and structure-function relationships of CPPs are not clear, although the incorporation of nonproteinogenic amino acids such as α-aminoisobutyric acid (Aib) has been reported to increase their helicity, biostability and penetration efficiencies. Here, we used molecular dynamics to study two Aib-containing CPPs, poly(LysAibAla)₃ (KAibA) and poly(LysAibGly)₃ (KAibG), that previously showed high cell internalization efficiency. KAibA and KAibG displayed the lowest internalization energies among the studied CPPs, showing distinct internalization mechanisms depending on the lipid composition of the model membranes. The presence of Aib residues allows these CPPs to adopt amphipathic folding to efficiently penetrate through the membranes. Elucidating how Aib incorporation affects CPP-membrane binding and interactions is beneficial for the design of CPPs for efficient intracellular delivery.

Impact of PROTAC Linker Plasticity on the Solution Conformations and Dissociation of the Ternary Complex

The conformational behavior of a small molecule free in solution is important to understand the free energy of binding to its target. This could be of special interest for proteolysis-targeting chimeras (PROTACs) due to their often flexible and lengthy linkers and the need to induce a ternary complex. Here, we report on the molecular dynamics (MD) simulations of two PROTACs, MZ1 and dBET6, revealing different linker conformational behaviors. The simulation of MZ1 in dimethyl sulfoxide (DMSO) agrees well with the nuclear magnetic resonance study, providing strong support for the relevance of our simulations. To further understand the role of linker plasticity in the formation of a ternary complex, the dissociation of the complex von Hippel-Lindau-MZ1-BRD4 is investigated in detail by steered simulations and is shown to follow a two-step pathway. Interestingly, both MZ1 and dBET6 display in water, a tendency toward an intramolecular lipophilic interaction between the two warheads. The hydrophobic contact of the two warheads would prevent them from binding to their respective proteins and might have an effect on the efficacy of induced cellular protein degradation. However, conformations featuring this hydrophobic contact of the two warheads are calculated to be marginally more favorable.

A computational study of competing conformational selection and induced fit in an abiotic system

Host-guest complexations can be described by two competing mechanisms, conformational selection (CS) and induced fit (IF). In this work, we used a combination of nudged elastic band (NEB), adaptive steered molecular dynamics (ASMD), and density functional theory (DFT, with a correction for dispersion) to study the dynamics of the pathways (IF/CS) by which two conformers of basket B(+) and B(-) interconvert and trap CX₄ guests (X = Cl and Br). While the results from NEB/DFT studies disclosed host-guest noncovalent contacts reducing the basket's conformational dynamics, ASMD methodology suggested an associative mechanism for the guest complexation. With theory in excellent agreement with experiments, NEB and ASMD emerge as the methods of choice for studying dynamics of supramolecular systems.

The role of interface on the toughening and failure mechanisms of thermoplastic nanocomposites reinforced with nanofibrillated rubber

The interface plays a crucial role in the physical and functional properties of polymer nanocomposites, yet its effects have not been fully recognized in the setting of classical continuum-based modeling. In the present study, we investigate the roles of interface and interfiber interactions on the toughening effects of rubber nanofibers embodied in thermoplastic-based materials. Emphasis is placed on establishing comprehensive theoretical and atomistic descriptions of the nanocomposite systems subjected to pull-out and uniaxial extension in the longitudinal and transverse directions. Using the framework of molecular dynamics, the annealed melt-drawn nanofibers were spontaneously formed via the proposed four-step methodology. The generated nanofibers were then crosslinked using the proposed robust topology-matching algorithm, through which the chemical reactions arising in the crosslinking were closely assimilated. The interfiber interactions were also examined with respect to separation distances and nanofiber radius via a nanofiber-pair atomistic scheme, and the obtained results were subsequently incorporated into the pull-out and uniaxial test simulations. The results indicate that the compatibilizer grafting results in enhanced interfacial shear strength by introducing extra chemical interactions at the interface. In particular, it was found that the compatibilizer restricts the formation and coalescence of nanovoids, resulting in enhanced toughening effects. Together, we have shown that the presence of a small amount of well-dispersed rubber nanofibrillar network whose surfaces are grafted with maleic anhydride compatibilizer can dramatically increase the toughness and alter the failure mechanisms of the nanocomposites without any deterioration in the stiffness, which is also consistent with the recent experimental observations in our lab. The interfacial failure mechanism was also investigated by monitoring the changes in the atomic concentration profiles, mean square displacement and fractional free volume. The results obtained may serve as a promising alternative for the continuum-based modeling and analysis of interfaces.

Searching and designing potential inhibitors for SARS-CoV-2 Mpro from natural sources using atomistic and deep-learning calculations

The spread of severe acute respiratory syndrome coronavirus 2 novel coronavirus (SARS-CoV-2) worldwide has caused the coronavirus disease 2019 (COVID-19) pandemic. A hundred million people were infected, resulting in several millions of death worldwide. In order to prevent viral replication, scientists have been aiming to prevent the biological activity of the SARS-CoV-2 main protease (3CL pro or Mpro). In this work, we demonstrate that using a reasonable combination of deep-learning calculations and atomistic simulations could lead to a new approach for developing SARS-CoV-2 main protease (Mpro) inhibitors. Initially, the binding affinities of the natural compounds to SARS-CoV-2 Mpro were estimated via atomistic simulations. The compound tomatine, thevetine, and tribuloside could bind to SARS-CoV-2 Mpro with nanomolar/high-nanomolar affinities. Secondly, the deep-learning (DL) calculations were performed to chemically alter the top-lead natural compounds to improve ligand-binding affinity. The obtained results were then validated by free energy calculations using atomistic simulations. The outcome of the research will probably boost COVID-19 therapy.

Direct Measurement of Intermolecular Mechanical Force for Nonspecific Interactions between Small Molecules

Mechanical force can evaluate intramolecular interactions in macromolecules. Because of the rapid motion of small molecules, it is extremely challenging to measure mechanical forces of nonspecific intermolecular interactions. Here, we used optical tweezers to directly examine the intermolecular mechanical force (IMMF) of nonspecific interactions between two cholesterols. We found that IMMFs of dimeric cholesterol complexes were dependent on the orientation of the interaction. The surprisingly high IMMF in cholesterol dimers (∼30 pN) is comparable to the mechanical stability of DNA secondary structures. Using Hess-like cycles, we quantified that changes in free energy of solubilizing cholesterol (ΔGsolubility) by β-cyclodextrin (βCD) and methylated βCD (Me-βCD) were as low as -16 and -27 kcal/mol, respectively. Compared to the ΔGsolubility of cholesterols in water (5.1 kcal/mol), these values indicated that cyclodextrins can easily solubilize cholesterols. Our results demonstrated that the IMMF can serve as a generic and multipurpose variable to dissect nonspecific intermolecular interactions among small molecules into orientational components.

Role of salt-bridging interactions in recognition of viral RNA by arginine-rich peptides

Interactions between RNA molecules and proteins are critical to many cellular processes and are implicated in various diseases. The RNA-peptide complexes are good model systems to probe the recognition mechanism of RNA by proteins. In this work, we report studies on the binding-unbinding process of a helical peptide from a viral RNA element using nonequilibrium molecular dynamics simulations. We explored the existence of various dissociation pathways with distinct free-energy profiles that reveal metastable states and distinct barriers to peptide dissociation. We also report the free-energy differences for each of the four pathways to be 96.47 ± 12.63, 96.1 ± 10.95, 91.83 ± 9.81, and 92 ± 11.32 kcal/mol. Based on the free-energy analysis, we further propose the preferred pathway and the mechanism of peptide dissociation. The preferred pathway is characterized by the formation of sequential hydrogen-bonding and salt-bridging interactions between several key arginine amino acids and the viral RNA nucleotides. Specifically, we identified one arginine amino acid (R8) of the peptide to play a significant role in the recognition mechanism of the peptide by the viral RNA molecule.

Reweighted Jarzynski Sampling: Acceleration of Rare Events and Free Energy Calculation with a Bias Potential Learned from Nonequilibrium Work

We introduce a simple enhanced sampling approach for the calculation of free energy differences and barriers along a one-dimensional reaction coordinate. First, a small number of short nonequilibrium simulations are carried out along the reaction coordinate, and the Jarzynski equality is used to learn an approximate free energy surface from the nonequilibrium work distribution. This free energy estimate is represented in a compact form as an artificial neural network and used as an external bias potential to accelerate rare events in a subsequent molecular dynamics simulation. The final free energy estimate is then obtained by reweighting the equilibrium probability distribution of the reaction coordinate sampled under the influence of the external bias. We apply our reweighted Jarzynski sampling recipe to four processes of varying scales and complexities─spanning chemical reaction in the gas phase, pair association in solution, and droplet nucleation in supersaturated vapor. In all cases, we find reweighted Jarzynski sampling to be a very efficient strategy, resulting in rapid convergence of the free energy to high precision.

Measurement of Single-Molecule Forces in Cholesterol and Cyclodextrin Host-Guest Complexes

Biological host molecules such as β-cyclodextrins (β-CDs) have been used to remove cholesterol guests from membranes and artery plaques. In this work, we calibrated the host-guest intermolecular mechanical forces (IMMFs) between cholesterol and cyclodextrin complexes by combining single-molecule force spectroscopy in optical tweezers and computational molecular simulations for the first time. Compared to native β-CD, methylated beta cyclodextrins complexed with cholesterols demonstrated higher mechanical stabilities due to the loss of more high-energy water molecules inside the methylated β-CD cavities. This result is consistent with the finding that methylated β-CD is more potent at solubilizing cholesterols than β-CD, suggesting that the IMMF can serve as a novel indicator to evaluate the solubility of small molecules such as cholesterols. Importantly, we found that the force spectroscopy measured in such biological host-guest complexes is direction-dependent: pulling from the alkyl end of the cholesterol molecule resulted in a larger IMMF than that from the hydroxyl end of the cholesterol molecule. Molecular dynamics coupled with umbrella sampling simulations further revealed that cholesterol molecules tend to enter or leave from the wide opening of cyclodextrins. Such an orientation rationalizes that cyclodextrins are rather efficient at extracting cholesterols from the phospholipid bilayer in which hydroxyl groups of cholesterols are readily exposed to the hydrophobic cavities of cyclodextrins. We anticipate that the IMMF measured by both experimental and computational force spectroscopy measurements help elucidate solubility mechanisms not only for cholesterols in different environments but also to host-guest systems in general, which have been widely exploited for their solubilization properties in drug delivery, for example.

Accelerating Reaction Rates of Biomolecules by Using Shear Stress in Artificial Capillary Systems

Biomimetics is a design principle within chemistry, biology, and engineering, but chemistry biomimetic approaches have been generally limited to emulating nature's chemical toolkit while emulation of nature's physical toolkit has remained largely unexplored. To begin to explore this, we designed biophysically mimetic microfluidic reactors with characteristic length scales and shear stresses observed within capillaries. We modeled the effect of shear with molecular dynamics studies and showed that this induces specific normally buried residues to become solvent accessible. We then showed using kinetics experiments that rates of reaction of these specific residues in fact increase in a shear-dependent fashion. We applied our results in the creation of a new microfluidic approach for the multidimensional study of cysteine biomarkers. Finally, we used our approach to establish dissociation of the therapeutic antibody trastuzumab in a reducing environment. Our results have implications for the efficacy of existing therapeutic antibodies in blood plasma as well as suggesting in general that biophysically mimetic chemistry is exploited in biology and should be explored as a research area.

Structural Basis of Salicylic Acid Decarboxylase Reveals a Unique Substrate Recognition Mode and Access Channel

Salicylic acid (SA) decarboxylase from Trichosporon moniliiforme (TmSdc), which reversibly catalyses the decarboxylation of SA to yield phenol, is of significant interest because of its potential for the production of benzoic acid derivatives under environmentally friendly reaction conditions. TmSdc showed a preference for C2 hydroxybenzoate derivatives, with k_cat/K_m of SA being 3.2 × 10³ M^-1 s^-1. Here, we presented the first crystal structures of TmSdc, including a complex with SA. The three conserved residues of Glu8, His169, and Asp298 are the catalytic residues within the TIM-barrel domain of TmSdc. Trp239 forms a unique hydrophobic recognition site by interacting with the phenyl ring of SA, while Arg235 is responsible for recognizing the hydroxyl group at the C2 of SA via hydrogen bond interactions. Using a semi-rational combinatorial active-site saturation test, we obtained the TmSdc mutant MT3 (Y64T/P191G/F195V/E302D), which exhibited a 26.4-fold increase in k_cat/K_m with SA, reaching 8.4 × 10⁴ M^-1 s^-1. Steered molecular dynamics simulations suggested that the structural changes in MT3 relieved the steric hindrance within the substrate access channel and enlarged the substrate-binding pocket, leading to the increased activity by improving substrate access. Our data elucidate the unique substrate recognition mode and the substrate entrance tunnel of SA decarboxylase.

Membrane-Facilitated Receptor Access and Binding Mechanisms of Long-Acting β2-Adrenergic Receptor Agonists

The drugs salmeterol, formoterol, and salbutamol constitute the frontline treatment of asthma and other chronic pulmonary diseases. These drugs activate the β2-adrenergic receptors (β2-AR), a class A G protein-coupled receptor (GPCR), and differ significantly in their clinical onset and duration of actions. According to the microkinetic model, the long duration of action of salmeterol and formoterol compared with salbutamol were attributed, at least in part, to their high lipophilicity and increased local concentrations in the membrane near the receptor. However, the structural and molecular bases of how the lipophilic drugs reach the binding site of the receptor from the surrounding membrane remain unknown. Using a variety of classic and enhanced molecular dynamics simulation techniques, we investigated the membrane partitioning characteristics, binding, and unbinding mechanisms of the ligands. The obtained results offer remarkable insight into the functional role of membrane lipids in the ligand association process. Strikingly, salmeterol entered the binding site from the bilayer through transmembrane helices 1 and 7. The entry was preceded by membrane-facilitated rearrangement and presentation of its phenyl-alkoxy-alkyl tail as a passkey to an access route gated by F193, a residue known to be critical for salmeterol's affinity. Formoterol's access is through the aqueous path shared by other β2-AR agents. We observed a novel secondary path for salbutamol that is distinct from its primary route. Our study offers a mechanistic description for the membrane-facilitated access and binding of ligands to a membrane protein and establishes a groundwork for recognizing membrane lipids as an integral component in the molecular recognition process. SIGNIFICANCE STATEMENT: The cell membrane's functional role behind the duration of action of long-acting β2-adrenergic receptor (β2-AR) agonists such as salmeterol has been a subject of debate for a long time. This study investigated the binding and unbinding mechanisms of the three commonly used β2-AR agonists, salmeterol, formoterol, and salbutamol, using advanced simulation techniques. The obtained results offer unprecedented insights into the active role of membrane lipids in facilitating access and binding of the ligands, affecting the molecular recognition process and thus their pharmacology.

Encapsulation of an anticancer drug Isatin inside a host nano-vehicle SWCNT: a molecular dynamics simulation

The use of carbon nanotubes as anticancer drug delivery cargo systems is a promising modality as they are able to perforate cellular membranes and transport the carried therapeutic molecules into the cellular components. Our work describes the encapsulation process of a common anticancer drug, Isatin (1H-indole-2,3-dione) as a guest molecule, in a capped single-walled carbon nanotube (SWCNT) host with chirality of (10,10). The encapsulation process was modelled, considering an aqueous solution, by a molecular dynamics (MD) simulation under a canonical NVT ensemble. The interactions between the atoms of Isatin were obtained from the DREIDING force filed. The storage capacity of the capped SWCNT host was evaluated to quantify its capacity to host multiple Isatin molecules. Our results show that the Isatin can be readily trapped inside the volume cavity of the capped SWCNT and it remained stable, as featured by a reduction in the van der Waals forces between Isatin guest and the SWCNT host (at approximately - 30 kcal mol^-1) at the end of the MD simulation (15 ns). Moreover, the free energy of encapsulation was found to be - 34 kcal mol^-1 suggesting that the Isatin insertion procedure into the SWCNT occurred spontaneously. As calculated, a capped SWCNT (10,10) with a length of 30 Å, was able to host eleven (11) molecules of Isatin, that all remained steadily encapsulated inside the SWCNT volume cavity, showing a potential for the use of carbon nanotubes as drug delivery cargo systems.

On the calculation of free energies over Hamiltonian and order parameters via perturbation and thermodynamic integration

In this work, complementary formulas are presented to compute free-energy differences via perturbation (FEP) methods and thermodynamic integration (TI). These formulas are derived by selecting only the most statistically significant data from the information extractable from the simulated points involved. On the one hand, commonly used FEP techniques based on overlap sampling leverage the full information contained in the overlapping macrostate probability distributions. On the other hand, conventional TI methods only use information on the first moments of those distributions, as embodied by the first derivatives of the free energy. Since the accuracy of simulation data degrades considerably for high-order moments (for FEP) or free-energy derivatives (for TI), it is proposed to consider, consistently for both methods, data up to second-order moments/derivatives. This provides a compromise between the limiting strategies embodied by common FEP and TI and leads to simple, optimized expressions to evaluate free-energy differences. The proposed formulas are validated with an analytically solvable harmonic Hamiltonian (for assessing systematic errors), an atomistic system (for computing the potential of mean force with coordinate-dependent order parameters), and a binary-component coarse-grained model (for tracing a solid-liquid phase diagram in an ensemble sampled through alchemical transformations). It is shown that the proposed FEP and TI formulas are straightforward to implement, perform similarly well, and allow robust estimation of free-energy differences even when the spacing of successive points does not guarantee them to have proper overlapping in phase space.

Steered molecular dynamics and stability analysis on PAH dimerisation and condensation on fullerene and soot surfaces

The mechanism of how a soot nucleus is impacted by polycyclic aromatic hydrocarbons (PAHs) and then grows through PAH condensation remains unclear. Using steered molecular dynamics (SMD), the non-bonding interaction between PAHs and soot was quantitatively studied using the free energy distribution during the dimerisation and condensation. The results showed that only two dimers (A₇-A₁₀ and 2 A₁₀) remained stable at 1000 K. The simulations showed that PAH condensation on a fullerene should not be ignored in soot mass growth. For fullerenes with a diameter not less than 1.8 nm (C₅₄₀), even A₄ condenses at temperatures of 1500 K, and A10 condenses stably on the surface of fullerenes even at 2000 K. The effects of multilayers and hydrogenated fullerenes on the free energy of PAH condensation are different. The stability of PAH dimers and PAH condensation pairs was discussed through free energy and chemical equilibrium. The results show that larger dimers are more stable than small ones at flame temperatures. Condensation is far more important than nucleation in mass growth at flame temperatures. Furthermore, the larger the PAH is, the higher the transformation ratio of the PAH in condensation on soot and thus the more stable the condensation product is. Finally, both the stability analysis of an upper temperature limit for condensation and simulation results of ReaxFF-MD cross-confirm that pyrene stably condensates on a simplified nascent soot (C₅₄₀) and a simulated soot (C₆₅₈H₃₁₉O₉), respectively, at 1500 K, but not at higher temperatures over 1800 K.

Theoretical Encapsulation of Fluorouracil (5-FU) Anti-Cancer Chemotherapy Drug into Carbon Nanotubes (CNT) and Boron Nitride Nanotubes (BNNT)

INTRODUCTION

Chemotherapy with anti-cancer drugs is considered the most common approach for killing cancer cells in the human body. However, some barriers such as toxicity and side effects would limit its usage. In this regard, nano-based drug delivery systems have emerged as cost-effective and efficient for sustained and targeted drug delivery. Nanotubes such as carbon nanotubes (CNT) and boron nitride nanotubes (BNNT) are promising nanocarriers that provide the cargo with a large inner volume for encapsulation. However, understanding the insertion process of the anti-cancer drugs into the nanotubes and demonstrating drug-nanotube interactions starts with theoretical analysis.

METHODS

First, interactions parameters of the atoms of 5-FU were quantified from the DREIDING force field. Second, the storage capacity of BNNT (8,8) was simulated to count the number of drugs 5-FU encapsulated inside the cavity of the nanotubes. In terms of the encapsulation process of the one drug 5-FU into nanotubes, it was clarified that the drug 5-FU was more rapidly adsorbed into the cavity of the BNNT compared with the CNT due to the higher van der Waals (vdW) interaction energy between the drug and the BNNT.

RESULTS

The obtained values of free energy confirmed that the encapsulation process of the drug inside the CNT and BNNT occurred spontaneously with the free energies of -14 and -25 kcal·mol-1, respectively.

DISCUSSION

However, the lower value of the free energy in the system containing the BNNT unraveled more stability of the encapsulated drug inside the cavity of the BNNT comparing the system having CNT. The encapsulation of Fluorouracil (5-FU) anti-cancer chemotherapy drug (commercial name: Adrucil®) into CNT (8,8) and BNNT (8,8) with the length of 20 Å in an aqueous solution was discussed herein applying molecular dynamics (MD) simulation.

Wetting characteristics of polymer adhesives with different chain bending stiffness

Polymer adhesives are widely used in daily applications and in industry owing to their flexibility and overall non-toxicity, particularly in interfacial adhesion. The spreading of polymer adhesives on adherend is one of the essential considerations for the interfacial adhesion of polymer adhesives, which is strongly related to their wetting behaviors. While relationships between polymer microstructure and adhesion have been investigated in previous studies, it remains challenging to unveil the effect of polymer microstructure on wettability. To address this issue, here we utilize coarse-grained molecular dynamics (CGMD) simulations to systematically elucidate how the wettability of a polymer adhesive droplet on a surface depends on bending stiffness. The wetting dynamics and the contact angle are studied to show the evolution of morphology of droplets during the wetting process. The results indicate the wettability is weakened by the increase of bending stiffness of polymer chain. Detailed thermodynamic property analysis is further conducted, revealing that the adhesion between the polymer droplet and substrate deteriorates due to the decline of wettability. Interestingly, we observe such deterioration becomes more significant by both increasing the temperature and decreasing the bending stiffness. Our study sheds light on the dependence of chain bending stiffness and temperature on the wetting behavior of polymer adhesive droplets, and offers insights, which, upon experimental validation can then be used for the design of adhesives or hydrogels.

Unravelling the inhibitory zinc ion binding site and the metal exchange mechanism in human DPP III

Dipeptidyl peptidase III (DPP III), a zinc-dependent exopeptidase, is widely distributed in organisms and present in almost all human tissues. In addition to its involvement in protein catabolism, it plays a role in oxidative stress and blood pressure regulation, and there is evidence of its involvement in pain modulation. Excess zinc ions have been found to inhibit its hydrolytic activity, but the binding affinity, binding site geometry, and mechanism of inhibitory activity have been unknown. Using several different computational approaches, we determined the inhibitory zinc ion binding site, its coordination and relative binding affinity. During some simulations the translocation of the zinc ion from the inhibitory to the catalytic binding site was observed, accompanied by movement of the catalytic zinc ion toward the exit of the substrate binding site. The traced behavior suggests an associative type of metal ion exchange, in which the formation of the ternary complex between enzyme and two metal ions precedes the exit of the catalytic metal ion. Differently from our previous findings that binding of a peptide induces partial opening of hDPP III, the globularity of the protein did not change in MD simulations of the hermorphin-like peptide bound to hDPP III with two zinc ions. However, the entrance to the interdomain cleft widens during Zn diffusion into the protein and was found to be the highest energy barrier in the process of metal translocation from the solvent to the active site. Finally, we discuss why excess zinc reduces enzyme activity.

Data-Driven Langevin Modeling of Nonequilibrium Processes

Given nonstationary data from molecular dynamics simulations, a Markovian Langevin model is constructed that aims to reproduce the time evolution of the underlying process. While at equilibrium the free energy landscape is sampled, nonequilibrium processes can be associated with a biased energy landscape, which accounts for finite sampling effects and external driving. When the data-driven Langevin equation (dLE) approach [Phys. Rev. Lett. 2015, 115, 050602] is extended to the modeling of nonequilibrium processes, an efficient way to calculate multidimensional Langevin fields is outlined. The dLE is shown to correctly account for various nonequilibrium processes, including the enforced dissociation of sodium chloride in water, the pressure-jump induced nucleation of a liquid of hard spheres, and the conformational dynamics of a helical peptide sampled from nonstationary short trajectories.

Review on Interfacial Bonding Mechanism of Functional Polymer Coating on Glass in Atomistic Modeling Perspective

Atomistic modeling methods are successfully applied to understand interfacial interaction in nanoscale size and analyze adhesion mechanism in the organic-inorganic interface. In this paper, we review recent representative atomistic simulation works, focusing on the interfacial bonding, adhesion strength, and failure behavior between polymer film and silicate glass. The simulation works are described under two categories, namely non-bonded and bonded interaction. In the works for non-bonded interaction, three main interactions, namely van der Waals interaction, polar interaction, and hydrogen bonds, are investigated, and the contributions to interfacial adhesion energy are analyzed. It is revealed that the most dominant interaction for adhesion is hydrogen bonding, but flexibility of the polymer film and modes of adhesion measurement test do affect adhesion and failure behavior. In the case of bonded interactions, the mechanism of covalent silane bond formation through condensation and hydrolysis process is reviewed, and surface reactivity, molecular density, and adhesion properties are calculated with an example of silane functionalized polymer. Besides interfacial interactions, effects of external conditions, such as surface morphology of the glass substrate and relative humidity on the adhesion and failure behavior, are presented, and modeling techniques developed for building interfacial system and calculating adhesion strengths are briefly introduced.

The transport of a charged peptide through carbon nanotubes under an external electric field: a molecular dynamics simulation

The study of interactions between biomolecules and carbon nanotubes (CNTs) is of great importance in CNT-based drug delivery systems and biomedical devices. In this work, the transport of polyarginine (R8) peptide through CNTs under an external electric field was investigated via all-atom molecular dynamics (AAMD) simulation. It was found that the electric field can assist the R8 peptide to overcome the resistance and make the transport smooth. Moreover, the efficiency of transport was improved with the increasing intensity of the electric field in a suitable range. In addition, we also investigated the effects of different types of CNTs on the transport of the R8 peptide and found that the single-walled carbon nanotube (SWCNT) was more suitable for transporting the R8 peptide than the double-walled carbon nanotube (DWCNT) due to its lower energy barrier to the R8 peptide. All these findings shed light on the role of the electric field on the transport of the R8 peptide through CNTs and also gave some valuable insights into the effects of CNT types on the transport process of the peptide.

Electrochemical adsorption of hydrogen on mixed Pd2Pt nanostructures

In the present contribution we have focused on the electrochemical adsorption of a proton from the solution-the Volmer reaction-on a variety of systems based on bimetallic nanostructures-clusters and wires-of Pd and Pt deposited on a surface of Au(111). We have calculated the free energy surface for the electron transfer step by a combination of DFT calculations, MD simulations and the theory of electrocatalysis. We analyze in detail the interaction of the metal d band with the valence orbital of the hydrogen and its effect on the catalytic activity as well as several aspects that influence the electrode reactivity such as spatial arrangements of the nanostructures, the solvation shell and chemical factors. We found that the mixed Pd₂Pt wire interacts strongly with hydrogen, and retains an almost complete solvation shell, which is reflected in a substantially reduced activation energy for the Volmer step. Thus, Pd₂Pt wires on Au(111) are predicted to be efficient electrocatalysts for the reaction.

α-Helical Antimicrobial Peptide Encapsulation and Release from Boron Nitride Nanotubes: A Computational Study

Introduction

Antimicrobial peptides are potential therapeutics as anti-bacteria, anti-viruses, anti-fungi, or anticancers. However, they suffer from a short half-life and drug resistance which limit their long-term clinical usage.

Methods

Herein, we captured the encapsulation of antimicrobial peptide HA-FD-13 into boron nitride nanotube (BNNT) (20,20) and its release due to subsequent insertion of BNNT (14,14) with molecular dynamics simulation.

Results

The peptide-BNNT (20,20) van der Waals (vdW) interaction energy decreased to −270 kcal·mol⁻¹ at the end of the simulation (15 ns). However, during the period of 0.2–1.8 ns, when half of the peptide was inside the nanotube, the encapsulation was paused due to an energy barrier in the vicinity of BNNT and subsequently the external intervention, such that the self-adjustment of the peptide allowed full insertion. The free energy of the encapsulation process was −200.12 kcal·mol⁻¹, suggesting that the insertion procedure occurred spontaneously.

Discussion

Once the BNNT (14,14) entered into the BNNT (20,20), the peptide was completely released after 83.8 ps. This revealed that the vdW interaction between the BNNT (14,14) and BNNT (20,20) was stronger than between BNNT (20,20) and the peptide; therefore, the BNNT (14,14) could act as a piston pushing the peptide outside the BNNT (20,20). Moreover, the sudden drop in the vdW energy between nanotubes to the value of the −1300 Kcal·mol⁻¹ confirmed the self-insertion of the BNNT (14,14) into the BNNT (20,20) and correspondingly the release of the peptide.

Computational methods for exploring protein conformations

Proteins are dynamic molecules that can transition between a potentially wide range of structures comprising their conformational ensemble. The nature of these conformations and their relative probabilities are described by a high-dimensional free energy landscape. While computer simulation techniques such as molecular dynamics simulations allow characterisation of the metastable conformational states and the transitions between them, and thus free energy landscapes, to be characterised, the barriers between states can be high, precluding efficient sampling without substantial computational resources. Over the past decades, a dizzying array of methods have emerged for enhancing conformational sampling, and for projecting the free energy landscape onto a reduced set of dimensions that allow conformational states to be distinguished, known as collective variables (CVs), along which sampling may be directed. Here, a brief description of what biomolecular simulation entails is followed by a more detailed exposition of the nature of CVs and methods for determining these, and, lastly, an overview of the myriad different approaches for enhancing conformational sampling, most of which rely upon CVs, including new advances in both CV determination and conformational sampling due to machine learning.

Benchmark of Popular Free Energy Approaches Revealing the Inhibitors Binding to SARS-CoV-2 Mpro

The COVID-19 pandemic has killed millions of people worldwide since its outbreak in December 2019. The pandemic is caused by the SARS-CoV-2 virus whose main protease (Mpro) is a promising drug target since it plays a key role in viral proliferation and replication. Currently, developing an effective therapy is an urgent task, which requires accurately estimating the ligand-binding free energy to SARS-CoV-2 Mpro. However, it should be noted that the accuracy of a free energy method probably depends on the protein target. A highly accurate approach for some targets may fail to produce a reasonable correlation with the experiment when a novel enzyme is considered as a drug target. Therefore, in this context, the ligand-binding affinity to SARS-CoV-2 Mpro was calculated via various approaches. The molecular docking approach was manipulated using Autodock Vina (Vina) and Autodock4 (AD4) protocols to preliminarily investigate the ligand-binding affinity and pose to SARS-CoV-2 Mpro. The binding free energy was then refined using the fast pulling of ligand (FPL), linear interaction energy (LIE), molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA), and free energy perturbation (FEP) methods. The benchmark results indicated that for docking calculations, Vina is more accurate than AD4, and for free energy methods, FEP is the most accurate method, followed by LIE, FPL, and MM-PBSA (FEP > LIE ≈ FPL > MM-PBSA). Moreover, atomistic simulations revealed that the van der Waals interaction is the dominant factor. The residues Thr26, His41, Ser46, Asn142, Gly143, Cys145, His164, Glu166, and Gln189 are essential elements affecting the binding process. Our benchmark provides guidelines for further investigations using computational approaches.

Helix-coil transition and conformational deformity in Aβ42-monomer: a case study using the Zn2+ cation

The metal ions (like Fe2+, Zn²⁺, Cu2+) are known to influence the amyloid beta (Aβ) aggregation. In this study, we have examined the conformational and dynamical changes during the coordination of Aβ-monomer with the Zn2+ ion using all-atom molecular dynamics (MD) simulations using explicit solvent models. We have probed the unfolding of the full-length Aβ₄₂ monomer both inclusive and exclusive of the Zn2+ cation, with 1:1 ratio of the peptide and the Zn2+ cation. The inclusion of the Zn2+ cation shows differential intra-peptide interactions which has been probed using various analyses. The Helix - Coil transition of the wild type Aβ42 monomer is studied using the steered molecular dynamics simulations by taking the end-to-end C-α distance across the peptide. This gives an idea of the unequal intra - peptide and peptide - water interactions being found across the length of the Aβ monomer. The transition of an α-helix dominated wild-type (WT) Aβ structure to the unfolded coil structure gives significant evidence of the intra-peptide hydrogen bonding shifts in the presence of the Zn2+ cation. This accounts for the structural and the dynamical variations that take place in the Aβ monomer in the presence of the Zn2+ cation to mimic the conditions/environment at the onset of fibrillation.Communicated by Ramaswamy H. Sarma.

Energetics and structure of alanine-rich α-helices via adaptive steered molecular dynamics

The energetics and hydrogen bonding profiles of the helix-to-coil transition were found to be an additive property and to increase linearly with chain length, respectively, in alanine-rich α-helical peptides. A model system of polyalanine repeats was used to establish this hypothesis for the energetic trends and hydrogen bonding profiles. Numerical measurements of a synthesized polypeptide Ac-Y(AEAAKA)kF-NH2 and a natural α-helical peptide a2N (1-17) provide evidence of the hypothesis's generality. Adaptive steered molecular dynamics was employed to investigate the mechanical unfolding of all of these alanine-rich polypeptides. We found that the helix-to-coil transition is primarily dependent on the breaking of the intramolecular backbone hydrogen bonds and independent of specific side-chain interactions and chain length. The mechanical unfolding of the α-helical peptides results in a turnover mechanism in which a 310-helical structure forms during the unfolding, remaining at a near constant population and thereby maintaining additivity in the free energy. The intermediate partially unfolded structures exhibited polyproline II helical structure as previously seen by others. In summary, we found that the average force required to pull alanine-rich α-helical peptides in between the endpoints-namely the native structure and free coil-is nearly independent of the length or the specific primary structure.

Further analyses of APRIL/APRIL-Receptor/Glycosaminoglycan interactions by biochemical assays linked to computational studies

A proliferation-inducing ligand (APRIL) is a member of the tumor necrosis factor superfamily. APRIL is quite unique in this superfamily for at least for two reasons: i) it binds to glycosaminoglycans (GAGs) via its positively charged N-terminus; ii) one of its signaling receptor, the transmembrane activator CAML interactor (TACI) was also reported to bind GAGs. Here, as provided by biochemical evidences with the use of an APRIL deletion mutant linked to computational studies, APRIL-GAG interaction involved other regions than the APRIL N-terminus. Preferential interaction of APRIL with heparin followed by chondroitin sulfate E were confirmed by in silico analysis. Both computational and experimental approaches did not reveal heparan sulfate binding to TACI. Together, computational results corroborated experiments contributing with atomistic details to the knowledge on this biologically relevant trimolecular system. Additionally, a high-throughput rigorous analysis of the free energy calculations data was performed to critically evaluate the applied computational methodologies.

Boron Nitride Nanotube as an Antimicrobial Peptide Carrier: A Theoretical Insight

INTRODUCTION

Nanotube-based drug delivery systems have received considerable attention because of their large internal volume to encapsulate the drug and the ability to penetrate tissues, cells, and bacteria. In this regard, understanding the interaction between the drug and the nanotube to evaluate the encapsulation behavior of the drug in the nanotube is of crucial importance.

METHODS

In this work, the encapsulation process of the cationic antimicrobial peptide named cRW3 in the biocompatible boron nitride nanotube (BNNT) was investigated under the Canonical ensemble (NVT) by molecular dynamics (MD) simulation.

RESULTS

The peptide was absorbed into the BNNT by van der Waals (vdW) interaction between cRW3 and the BNNT, in which the vdW interaction decreased during the simulation process and reached the value of -142.7 kcal·mol-1 at 4 ns.

DISCUSSION

The increase in the potential mean force profile of the encapsulated peptide during the pulling process of cRW3 out of the nanotube showed that its insertion into the BNNT occurred spontaneously and that the inserted peptide had the desired stability. The energy barrier at the entrance of the BNNT caused a pause of 0.45 ns when half of the peptide was inside the BNNT during the encapsulation process. Therefore, during this period, the peptide experienced the weakest movement and the smallest conformational changes.

Disaggregation mechanism of prion amyloid for tweezer inhibitor

The aggregation of amyloid has been an important event in the pathology of amyloidogenicity. A number of small molecules have been designed for Amyloidosis treatment. Molecular tweezer CLR01, a potential drug for misfolded β-amyloids inhibition, was reportedly bind directly to Lysine residues and interrupt oligomerization. However, the disaggregation mechanism of amyloid for this inhibitor is unclear. Here we used long timescale of molecular dynamic simulation to reveal the mechanism of disaggregation for pentamer prion amyloid. Molecular docking and molecular dynamics simulation demonstrate that CLR01 is attached with Lysine₂₂₂ nitrogen by π-cation interaction of its nine aromatic rings and formation of salt bridge/hydrogen bond of one of the two rotatable peripheral anionic phosphate groups. Upon CLR01 binding, we found a major shifting occurs in initial conformation of the oligomer and stretch out the N-terminal chain A from the rest of the amyloid which seems to be the first stage of disaggregated the fibrils slowly yet efficiently. Moreover, the CLR01 remodelled the pentamer Prion_220-272 into a compact structure which might be the resistant conformation for further oligomerization. Our work will contribute to better understand the interaction and deterioration mechanism of molecular tweezer for prions and similar amyloids, and offer significant insights into therapeutic development for Amyloidosis treatment.

Potential of mean force and umbrella sampling simulation for the transport of 5-oxazolidinone in heterotetrameric sarcosine oxidase

The structure of heterotetrameric sarcosine oxidase (HSO) contains a highly complex system composed of a large cavity and tunnels, which are essential for the reaction and migration of the reactants, products, and intermediates. Previous geometrical analysis using the CAVER program has predicted that there are three possible tunnels, T1, T2, and T3, for the exit pathway of the iminium intermediate, 5-oxazolidinone (5-OXA), of the enzyme reaction. Previous molecular dynamics (MD) simulation of HSO has identified the regions containing the water channels from the density distribution of water. The simulation indicated that tunnel T3 is the most probable exit pathway of 5-OXA. In the present study, the potential of mean force (PMF) for the transport of 5-OXA through tunnels T1, T2, and T3 was calculated using umbrella sampling (US) MD simulations and the weighted histogram analysis method. The PMF profiles for the three tunnels support the notion that tunnel T3 is the exit pathway of 5-OXA, and that 5-OXA tends to stay at the middle of the tunnel. The maximum errors of the calculated PMF for the predicted exit pathway, tunnel T3, were estimated by repeating the US simulations using different sets of initial positions. The PMF profile was also calculated for the transport of glycine within T3. The PMF profiles from the US simulations were in good agreement with the previous predictions that 5-OXA escape through tunnel T3 and how glycine is released to the outside of HSO was discussed.