Hepatic SEL1L-HRD1 ER-associated degradation regulates systemic iron homeostasis via ceruloplasmin.

Iron homeostasis is critical for cellular and organismal function and is tightly regulated to prevent toxicity or anemia due to iron excess or deficiency, respectively. However, subcellular regulatory mechanisms of iron remain largely unexplored. Here, we report that SEL1L-HRD1 protein complex of endoplasmic reticulum (ER)-associated degradation (ERAD) in hepatocytes controls systemic iron homeostasis in a ceruloplasmin (CP)-dependent, and ER stress-independent, manner. Mice with hepatocyte-specific Sel1L deficiency exhibit altered basal iron homeostasis and are sensitized to iron deficiency while resistant to iron overload. Proteomics screening for a factor linking ERAD deficiency to altered iron homeostasis identifies CP, a key ferroxidase involved in systemic iron distribution by catalyzing iron oxidation and efflux from tissues. Indeed, CP is highly unstable and a bona fide substrate of SEL1L-HRD1 ERAD. In the absence of ERAD, CP protein accumulates in the ER and is shunted to refolding, leading to elevated secretion. Providing clinical relevance of these findings, SEL1L-HRD1 ERAD is responsible for the degradation of a subset of disease-causing CP mutants, thereby attenuating their pathogenicity. Together, this study uncovers the role of SEL1L-HRD1 ERAD in systemic iron homeostasis and provides insights into protein misfolding-associated proteotoxicity.

Mechanistic insights into ligand dissociation from the SARS-CoV-2 spike glycoprotein.

The COVID-19 pandemic, driven by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has spurred an urgent need for effective therapeutic interventions. The spike glycoprotein of the SARS-CoV-2 is crucial for infiltrating host cells, rendering it a key candidate for drug development. By interacting with the human angiotensin-converting enzyme 2 (ACE2) receptor, the spike initiates the infection of SARS-CoV-2. Linoleate is known to bind the spike glycoprotein, subsequently reducing its interaction with ACE2. However, the detailed mechanisms underlying the protein-ligand interaction remain unclear. In this study, we characterized the pathways of ligand dissociation and the conformational changes associated with the spike glycoprotein by using ligand Gaussian accelerated molecular dynamics (LiGaMD). Our simulations resulted in eight complete ligand dissociation trajectories, unveiling two distinct ligand unbinding pathways. The preference between these two pathways depends on the gate distance between two α-helices in the receptor binding domain (RBD) and the position of the N-linked glycan at N343. Our study also highlights the essential contributions of K417, N121 glycan, and N165 glycan in ligand unbinding, which are equally crucial in enhancing spike-ACE2 binding. We suggest that the presence of the ligand influences the motions of these residues and glycans, consequently reducing accessibility for spike-ACE2 binding. These findings enhance our understanding of ligand dissociation from the spike glycoprotein and offer significant implications for drug design strategies in the battle against COVID-19.

Molecular dynamics study reveals key disruptors of MEIG1-PACRG interaction.

Interactions between the meiosis-expressed gene 1 (MEIG1) and Parkin co-regulated gene (PACRG) protein are critical in the formation of mature sperm cells. Targeting either MEIG1 or PACRG protein could be a contraceptive strategy. The W50A and Y68A mutations on MEIG1 are known to interrupt the MEIG1-PACRG interactions resulting in defective sperm cells. However, the details about how the mutants disrupt the protein-protein binding are not clear. In this study, we reveal insights on MEIG1 and PACRG protein dynamics by applying Gaussian-accelerated molecular dynamics (GaMD) simulations and post-GaMD analysis. Our results show that the mutations destabilize the protein-protein interfacial interaction. The effect of the Y68A mutation is more significant than W50A as Y68 forms stronger polar interactions with PACRG. Because both human and mouse models demonstrate similar dynamic properties, the findings from mouse proteins can be applied to the human system. Moreover, we report a potential ligand binding pocket on the MEIG1 and PACRG interaction surface that could be a target for future drug design to inhibit the MEIG1-PACRG interaction. PACRG shows more qualified pockets along the protein-protein interface, implying that it is a better target than MEIG1. Our work provides a fundamental understanding of MEIG1 and PACRG protein dynamics, paving the way for drug discovery in male-based contraception.

Multiscale computational study of ligand binding pathways: Case of p38 MAP kinase and its inhibitors.

Protein kinases are one of the most important drug targets in the past 10 years. Understanding the inhibitor association processes will profoundly impact new binder designs with preferred binding kinetics. However, after more than a decade of effort, a complete atomistic-level study of kinase inhibitor binding pathways is still lacking. As all kinases share a similar scaffold, we used p38 kinase as a model system to investigate the conformational dynamics and free energy transition of inhibitor binding toward kinases. Two major kinase conformations, Asp-Phe-Gly (DFG)-in and DFG-out, and three types of inhibitors, type I, II, and III, were thoroughly investigated in this work. We performed Brownian dynamics simulations and up to 340 µs Gaussian-accelerated molecular dynamics simulations to capture the inhibitor binding paths and a series of conformational transitions of the p38 kinase from its apo to inhibitor-bound form. Eighteen successful binding trajectories, including all types of inhibitors, are reported herein. Our simulations suggest a mechanism of inhibitor recruitment, a faster ligand association step to a pre-existing DFG-in/DFG-out p38 protein, followed by a slower molecular rearrangement step to adjust the protein-ligand conformation followed by a shift in the energy landscape to reach the final bound state. The ligand association processes also reflect the energetic favor of type I and type II/III inhibitor binding through ATP and allosteric channels, respectively. These different binding routes are directly responsible for the fast (type I binders) and slow (type II/III binders) kinetics of different types of p38 inhibitors. Our findings also echo the recent study of p38 inhibitor dissociation, implying that ligand unbinding could undergo a reverse path of binding, and both processes share similar metastates. This study deepens the understanding of molecular and energetic features of kinase inhibitor-binding processes and will inspire future drug development from a kinetic point of view.

Enhanced Disrupting Effect of Benzophenone-1 Chlorination Byproducts to the Androgen Receptor: Cell-Based Assays and Gaussian Accelerated Molecular Dynamics Simulations.

Benzophenone-1 (BP-1), one of the commonly used ultraviolet filters, has caused increasing public concern due to frequently detected residues in environmental and recreational waters. Its susceptibility to residual chlorine and the potential to subsequently trigger endocrine disruption remain unknown. We herein investigated the chlorination of BP-1 in swimming pool water and evaluated the endocrine disruption toward the human androgen receptor (AR). The structures of monochlorinated (P1) and dichlorinated (P2) products were separated and characterized by mass spectrometry and 1H-1H NMR correlation spectroscopy. P1 and P2 exhibited significantly higher antiandrogenic activity in yeast two-hybrid assays (EC50, 6.13 µM and 9.30 µM) than did BP-1 (12.89 µM). Our 350 ns Gaussian accelerated molecular dynamics simulations showed the protein dynamics in a long-time scale equilibrium, and further energy calculations revealed that although increased hydrophobic interactions are primarily responsible for enhanced binding affinities between chlorinated products and the AR ligand binding domain, the second chloride in P2 still hinders the complex motion because of the solvation penalty. The mixture of BP-1-P1-P2 elicited additive antiandrogenic activity, well fitted by the concentration addition model. P1 and P2 at 1 µM consequently downregulated the mRNA expression of AR-regulated genes, NKX3.1 and KLK3, by 1.7-9.1-fold in androgen-activated LNCaP cells. Because chlorination of BP-1 occurs naturally by residual chlorine in aquatic environments, our results regarding enhanced antiandrogenic activity and disturbed AR signaling provided evidence linking the use of personal care products with potential health risks.

Endosidin2 targets conserved exocyst complex subunit EXO70 to inhibit exocytosis.

The exocyst complex regulates the last steps of exocytosis, which is essential to organisms across kingdoms. In humans, its dysfunction is correlated with several significant diseases, such as diabetes and cancer progression. Investigation of the dynamic regulation of the evolutionarily conserved exocyst-related processes using mutants in genetically tractable organisms such as Arabidopsis thaliana is limited by the lethality or the severity of phenotypes. We discovered that the small molecule Endosidin2 (ES2) binds to the EXO70 (exocyst component of 70 kDa) subunit of the exocyst complex, resulting in inhibition of exocytosis and endosomal recycling in both plant and human cells and enhancement of plant vacuolar trafficking. An EXO70 protein with a C-terminal truncation results in dominant ES2 resistance, uncovering possible distinct regulatory roles for the N terminus of the protein. This study not only provides a valuable tool in studying exocytosis regulation but also offers a potentially new target for drugs aimed at addressing human disease.

Ligand Binding Pathways and Conformational Transitions of the HIV Protease.

It is important to determine the binding pathways and mechanisms of ligand molecules to target proteins to effectively design therapeutic drugs. Molecular dynamics (MD) is a promising computational tool that allows us to simulate protein-drug binding at an atomistic level. However, the gap between the time scales of current simulations and those of many drug binding processes has limited the usage of conventional MD, which has been reflected in studies of the HIV protease. Here, we have applied a robust enhanced simulation method, Gaussian accelerated molecular dynamics (GaMD), to sample binding pathways of the XK263 ligand and associated protein conformational changes in the HIV protease. During two of 10 independent GaMD simulations performed over 500-2500 ns, the ligand was observed to successfully bind to the protein active site. Although GaMD-derived free energy profiles were not fully converged because of insufficient sampling of the complex system, the simulations still allowed us to identify relatively low-energy intermediate conformational states during binding of the ligand to the HIV protease. Relative to the X-ray crystal structure, the XK263 ligand reached a minimum root-mean-square deviation (RMSD) of 2.26 Å during 2.5 µs of GaMD simulation. In comparison, the ligand RMSD reached a minimum of only ∼5.73 Å during an earlier 14 µs conventional MD simulation. This work highlights the enhanced sampling power of the GaMD approach and demonstrates its wide applicability to studies of drug-receptor interactions for the HIV protease and by extension many other target proteins.

Mechanism of the Association Pathways for a Pair of Fast and Slow Binding Ligands of HIV-1 Protease.

Equilibrium constants, together with kinetic rate constants of binding, are key factors in the efficacy and safety of drug compounds, informing drug design. However, the association pathways of protein-ligand binding, which contribute to their kinetic behaviors, are little understood. In this work, we used unbiased all-atom molecular dynamics (MD) simulations with an explicit solvent model to study the association processes of protein-ligand binding. Using the HIV protease (HIVp)-xk263 and HIVp-ritonavir protein-ligand systems as cases, we observed that ligand association is a multistep process involving diffusion, localization, and conformational rearrangements of the protein, ligand, and water molecules. Moreover, these two ligands preferred different routes of binding, which reflect two well-known binding mechanisms: induced-fit and conformation selection models. Our study shows that xk263 has a stronger capacity for desolvating surrounding water molecules, thereby inducing a semiopen conformation of the HIVp flaps (induced-fit model). In contrast, the slow dehydration characteristic of ritonavir allows for gradual association with the binding pocket of HIVp when the protein's flap conformation is fully open (conformation selection model). By studying the mechanism of ligand association and understanding the role of solvent molecules during the binding event, we can obtain a different perspective on the mechanism of macromolecule recognition, providing insights into drug discovery.

Visualizing the tunnel in tryptophan synthase with crystallography: Insights into a selective filter for accommodating indole and rejecting water.

Four new X-ray structures of tryptophan synthase (TS) crystallized with varying numbers of the amphipathic N-(4'-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F6) molecule are presented. These structures show one of the F6 ligands threaded into the tunnel from the ß-site and reveal a distinct hydrophobic region. Over this expanse, the interactions between F6 and the tunnel are primarily nonpolar, while the F6 phosphoryl group fits into a polar pocket of the ß-subunit active site. Further examination of TS structures reveals that one portion of the tunnel (T1) binds clusters of water molecules, whereas waters are not observed in the nonpolar F6 binding region of the tunnel (T2). MD simulation of another TS structure with an unobstructed tunnel also indicates the T2 region of the tunnel excludes water, consistent with a dewetted state that presents a significant barrier to the transfer of water into the closed ß-site. We conclude that hydrophobic molecules can freely diffuse between the α- and ß-sites via the tunnel, while water does not. We propose that exclusion of water serves to inhibit reaction of water with the α-aminoacrylate intermediate to form ammonium ion and pyruvate, a deleterious side reaction in the αß-catalytic cycle. Finally, while most TS structures show ßPhe280 partially blocking the tunnel between the α- and ß-sites, new structures show an open tunnel, suggesting the flexibility of the ßPhe280 side chain. Flexible docking studies and MD simulations confirm that the dynamic behavior of ßPhe280 allows unhindered transfer of indole through the tunnel, therefore excluding a gating role for this residue.

Characterization of Promiscuous Binding of Phosphor Ligands to Breast-Cancer-Gene 1 (BRCA1) C-Terminal (BRCT): Molecular Dynamics, Free Energy, Entropy and Inhibitor Design.

Inhibition of the protein-protein interaction (PPI) mediated by breast-cancer-gene 1 C-terminal (BRCT) is an attractive strategy to sensitize breast and ovarian cancers to chemotherapeutic agents that induce DNA damage. Such inhibitors could also be used for studies to understand the role of this PPI in DNA damage response. However, design of BRCT inhibitors is challenging because of the inherent flexibility associated with this domain. Several studies identified short phosphopeptides as tight BRCT binders. Here we investigated the thermodynamic properties of 18 phosphopeptides or peptide with phosphate mimic and three compounds with phosphate groups binding to BRCT to understand promiscuous molecular recognition and guide inhibitor design. We performed molecular dynamics (MD) simulations to investigate the interactions between inhibitors and BRCT and their dynamic behavior in the free and bound states. MD simulations revealed the key role of loops in altering the shape and size of the binding site to fit various ligands. The mining minima (M2) method was used for calculating binding free energy to explore the driving forces and the fine balance between configuration entropy loss and enthalpy gain. We designed a rigidified ligand, which showed unfavorable experimental binding affinity due to weakened enthalpy. This was because it lacked the ability to rearrange itself upon binding. Investigation of another phosphate group containing compound, C1, suggested that the entropy loss can be reduced by preventing significant narrowing of the energy well and introducing multiple new compound conformations in the bound states. From our computations, we designed an analog of C1 that introduced new intermolecular interactions to strengthen attractions while maintaining small entropic penalty. This study shows that flexible compounds do not always encounter larger entropy penalty, compared with other more rigid binders, and highlights a new strategy for inhibitor design.

Protonation states of the tryptophan synthase internal aldimine active site from solid-state NMR spectroscopy: direct observation of the protonated Schiff base linkage to pyridoxal-5'-phosphate.

The acid-base chemistry that drives catalysis in pyridoxal-5'-phosphate (PLP)-dependent enzymes has been the subject of intense interest and investigation since the initial identification of PLP's role as a coenzyme in this extensive class of enzymes. It was first proposed over 50 years ago that the initial step in the catalytic cycle is facilitated by a protonated Schiff base form of the holoenzyme in which the linking lysine ε-imine nitrogen, which covalently binds the coenzyme, is protonated. Here we provide the first (15)N NMR chemical shift measurements of such a Schiff base linkage in the resting holoenzyme form, the internal aldimine state of tryptophan synthase. Double-resonance experiments confirm the assignment of the Schiff base nitrogen, and additional (13)C, (15)N, and (31)P chemical shift measurements of sites on the PLP coenzyme allow a detailed model of coenzyme protonation states to be established.

Switches of hydrogen bonds during ligand-protein association processes determine binding kinetics.

Revealing the processes of ligand-protein associations deepens our understanding of molecular recognition and binding kinetics. Hydrogen bonds (H-bonds) play a crucial role in optimizing ligand-protein interactions and ligand specificity. In addition to the formation of stable H-bonds in the final bound state, the formation of transient H-bonds during binding processes contributes binding kinetics that define a ligand as a fast or slow binder, which also affects drug action. However, the effect of forming the transient H-bonds on the kinetic properties is little understood. Guided by results from coarse-grained Brownian dynamics simulations, we used classical molecular dynamics simulations in an implicit solvent model and accelerated molecular dynamics simulations in explicit waters to show that the position and distribution of the H-bond donor or acceptor of a drug result in switching intermolecular and intramolecular H-bond pairs during ligand recognition processes. We studied two major types of HIV-1 protease ligands: a fast binder, xk263, and a slow binder, ritonavir. The slow association rate in ritonavir can be attributed to increased flexibility of ritonavir, which yields multistep transitions and stepwise entering patterns and the formation and breaking of complex H-bond pairs during the binding process. This model suggests the importance of conversions of spatiotemporal H-bonds during the association of ligands and proteins, which helps in designing inhibitors with preferred binding kinetics.

Protein binding for detection of small changes on a nanoparticle surface.

Protein adsorption on nanoparticles is closely associated with the physicochemical properties of particles, in particular, their surface properties. We synthesized two batches of polyacrylic acid-coated nanoparticles under almost identical conditions except for the heating duration and found differences in the head-group structure of the polyacrylic acid. The structure change was confirmed by NMR and MS. The two batches of particles had varied binding affinities to a selected group of proteins. Computational work confirmed that the head group of the polymer on the surface of a nanoparticle could directly interact with a protein, and small structural changes in the head group were sufficient to result in a significant difference in the free energy of binding. Our results demonstrate that protein adsorption is so sensitive to the surface properties of particles that it can reveal even small variations in the structure of a nanoparticle surface ligand, and should be useful for quick assessment of nanoparticle properties.

Multiple Parameter Replica Exchange Gaussian Accelerated Molecular Dynamics for Enhanced Sampling and Free Energy Calculation of Biomolecular Systems.

This study introduces a novel method named multiple parameter replica exchange Gaussian accelerated molecular dynamics (MP-Rex-GaMD), building on the Gaussian accelerated molecular dynamics (GaMD) algorithm. GaMD enhances sampling and retrieves free energy information for biomolecular systems by adding a harmonic boost potential to smooth the potential energy surface without the need for predefined reaction coordinates. Our innovative approach advances the acceleration power and energetic reweighting accuracy of GaMD by incorporating a replica exchange algorithm that enables the exchange of multiple parameters, including the GaMD boost parameters of force constant and energy threshold, as well as temperature. Applying MP-Rex-GaMD to the three model systems of dialanine, chignolin, and HIV protease, we demonstrate its superior capability over conventional molecular dynamics and GaMD simulations in exploring protein conformations and effectively navigating various biomolecular states across energy barriers. MP-Rex-GaMD allows users to accurately map free energy landscapes through energetic reweighting, capturing the ensemble of biomolecular states from low-energy conformations to rare high-energy transitions within practical computational time scales.

Exploring additives beyond phthalates: Release from plastic mulching films, biodegradation and occurrence in agricultural soils.

It is widely recognized that applications of plastic films result in plastic pollution in agroecosystems. However, there is limited knowledge on the release and occurrence of additives beyond phthalates in agricultural soil. In this study, the rates of release and biodegradation of various additives, including phthalates, bisphenols, organophosphate esters, phenolic antioxidants, and ultraviolet absorbents from mulching films in soil were quantified by laboratory incubation. The rates of release and biodegradation ranged from 0.069 d-1 to 5.893 d-1 and from 1.43 × 10-3 d-1 to 0.600 d-1, respectively. Both of these rates were affected by temperature, flooding, and the properties of additives, films, and soils. An estimated 4000 metric tons of these additives were released into soil annually in China exclusively. The total concentrations of these additives in 80 agricultural soils varied between 228 and 3455 µg kg-1, with phenolic antioxidants, phthalates, and bisphenols accounting for 54.1%, 25.2%, and 17.9% of the total concentrations, respectively. A preliminary risk assessment suggested that the current levels of these additives could potentially present moderate hazards to the soil ecosystem.

Insights from free-energy calculations: protein conformational equilibrium, driving forces, and ligand-binding modes.

Accurate free-energy calculations provide mechanistic insights into molecular recognition and conformational equilibrium. In this work, we performed free-energy calculations to study the thermodynamic properties of different states of molecular systems in their equilibrium basin, and obtained accurate absolute binding free-energy calculations for protein-ligand binding using a newly developed M2 algorithm. We used a range of Asp-Phe-Gly (DFG)-in/out p38α mitogen-activated protein kinase inhibitors as our test cases. We also focused on the flexible DFG motif, which is closely connected to kinase activation and inhibitor binding. Our calculations explain the coexistence of DFG-in and DFG-out states of the loop and reveal different components (e.g., configurational entropy and enthalpy) that stabilize the apo p38α conformations. To study novel ligand-binding modes and the key driving forces behind them, we computed the absolute binding free energies of 30 p38α inhibitors, including analogs with unavailable experimental structures. The calculations revealed multiple stable, complex conformations and changes in p38α and inhibitor conformations, as well as balance in several energetic terms and configurational entropy loss. The results provide relevant physics that can aid in designing inhibitors and understanding protein conformational equilibrium. Our approach is fast for use with proteins that contain flexible regions for structure-based drug design.

Molecular Modeling of ABHD5 Structure and Ligand Recognition.

Alpha/beta hydrolase domain-containing 5 (ABHD5), also termed CGI-58, is the key upstream activator of adipose triglyceride lipase (ATGL), which plays an essential role in lipid metabolism and energy storage. Mutations in ABHD5 disrupt lipolysis and are known to cause the Chanarin-Dorfman syndrome. Despite its importance, the structure of ABHD5 remains unknown. In this work, we combine computational and experimental methods to build a 3D structure of ABHD5. Multiple comparative and machine learning-based homology modeling methods are used to obtain possible models of ABHD5. The results from Gaussian accelerated molecular dynamics and experimental data of the apo models and their mutants are used to select the most likely model. Moreover, ensemble docking is performed on representative conformations of ABHD5 to reveal the binding mechanism of ABHD5 and a series of synthetic ligands. Our study suggests that the ABHD5 models created by deep learning-based methods are the best candidate structures for the ABHD5 protein. The mutations of E41, R116, and G328 disturb the hydrogen bonding network with nearby residues and suppress membrane targeting or ATGL activation. The simulations also reveal that the hydrophobic interactions are responsible for binding sulfonyl piperazine ligands to ABHD5. Our work provides fundamental insight into the structure of ABHD5 and its ligand-binding mode, which can be further applied to develop ABHD5 as a therapeutic target for metabolic disease and cancer.

An Integrated Pharmacological, Structural, and Genetic Analysis of Extracellular Versus Intracellular ROS Production in Neutrophils.

The neutrophil NADPH oxidase produces both intracellular and extracellular reactive oxygen species (ROS). Although oxidase activity is essential for microbial killing, and ROS can act as signaling molecules in the inflammatory process, excessive extracellular ROS directly contributes to inflammatory tissue damage, as well as to cancer progression and immune dysregulation in the tumor microenvironment. How specific signaling pathways contribute to ROS localization is unclear. Here we used a systems pharmacology approach to identify the specific Class I PI3-K isoform p110ß, and PLD1, but not PLD2, as critical regulators of extracellular, but not intracellular ROS production in primary neutrophils. Combined crystallographic and molecular dynamics analysis of the PX domain of the oxidase component p47phox, which binds the lipid products of PI 3-K and PLD, was used to clarify the membrane-binding mechanism and guide the design of mutant mice whose p47phox is unable to bind 3-phosphorylated inositol phospholipids. Neutrophils from these K43A mutant animals were specifically deficient in extracellular, but not intracellular, ROS production, and showed increased dependency on signaling through the remaining PLD1 arm. These findings identify the PX domain of p47phox as a critical integrator of PLD1 and p110ß signaling for extracellular ROS production, and as a potential therapeutic target for modulating tissue damage and extracellular signaling during inflammation.

Gaussian accelerated molecular dynamics (GaMD): principles and applications.

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

Elucidating degradation mechanisms of florfenicol in soil by stable-isotope assisted nontarget screening.

Antibiotics in soil environments are a growing concern. Identifying transformation products is key to elucidating degradation pathways and mechanisms of antibiotics and other organic micropollutants. The primary challenge of transformation product identification is the interference of matrices. In this study, stable-isotope assisted nontarget screening was used to identify biodegradation products of florfenicol in soil. A total of 74 candidates were prioritized from thousands of mass features observed by a tiered peak filtering approach. Moreover, with the support of in silico prediction tools, the structures of 12 transformation products were elucidated, and 9 of them were reported for the first time. A biodegradation map of florfenicol consisting of amide hydrolysis, dechlorination, dehydration, defluorination, and sulfone reduction was established based on these identified products. A total of 8 products were also found in 6 field soil samples with manure application. Because of the structural similarity to florfenicol, some transformation products might still keep antimicrobial activity toward a variety of bacterial species. The strategies demonstrated in this study provide a basis for efficient identification of transformation products of other organic micropollutants in a variety of environmental matrices. The results also shed light on the degradation mechanisms, risk assessments, and regulations of these compounds.