Mitigation of acute kidney injury by cell-cycle inhibitors that suppress both CDK4/6 and OCT2 functions

Acute kidney injury (AKI) is a potentially fatal syndrome characterized by a rapid decline in kidney function caused by ischemic or toxic injury to renal tubular cells. The widely used chemotherapy drug cisplatin accumulates preferentially in the renal tubular cells and is a frequent cause of drug-induced AKI. During the development of AKI the quiescent tubular cells reenter the cell cycle. Strategies that block cell-cycle progression ameliorate kidney injury, possibly by averting cell division in the presence of extensive DNA damage. However, the early signaling events that lead to cell-cycle activation during AKI are not known. In the current study, using mouse models of cisplatin nephrotoxicity, we show that the G1/S-regulating cyclin-dependent kinase 4/6 (CDK4/6) pathway is activated in parallel with renal cell-cycle entry but before the development of AKI. Targeted inhibition of CDK4/6 pathway by small-molecule inhibitors palbociclib (PD-0332991) and ribociclib (LEE011) resulted in inhibition of cell-cycle progression, amelioration of kidney injury, and improved overall survival. Of additional significance, these compounds were found to be potent inhibitors of organic cation transporter 2 (OCT2), which contributes to the cellular accumulation of cisplatin and subsequent kidney injury. The unique cell-cycle and OCT2-targeting activities of palbociclib and LEE011, combined with their potential for clinical translation, support their further exploration as therapeutic candidates for prevention of AKI.Cell division is a fundamental biological process that is tightly regulated by evolutionarily conserved signaling pathways (1, 2). The initial decision to start cell division, the fidelity of subsequent DNA replication, and the final formation of daughter cells is monitored and regulated by these essential pathways (2–6). The cyclin-dependent kinases (CDKs) are the central players that orchestrate this orderly progression through the cell cycle (1, 2, 6, 7). The enzymatic activity of CDKs is regulated by complex mechanisms that include posttranslational modifications and expression of activating and inhibitory proteins (1, 2, 6, 7). The spatial and temporal changes in the activity of these CDK complexes are thought to generate the distinct substrate specificities that lead to sequential and unidirectional progression of the cell cycle (1, 8, 9).Cell-cycle deregulation is a universal feature of human cancer and a long-sought-after target for anticancer therapy (1, 10–13). Frequent genetic or epigenetic changes in mitogenic pathways, CDKs, cyclins, or CDK inhibitors are observed in various human cancers (1, 4, 11). In particular, the G1/S-regulating CDK4/6–cyclin D–inhibitors of CDK4 (INK4)–retinoblastoma (Rb) protein pathway frequently is disrupted in cancer cells (11, 14). These observations provided an impetus to develop CDK inhibitors as anticancer drugs. However, the earlier class of CDK inhibitors had limited specificity, inadequate clinical activity, poor pharmacokinetic properties, and unacceptable toxicity profiles (10, 11, 14, 15). These disappointing initial efforts now have been followed by the development of the specific CDK4/6 inhibitors palbociclib (PD0332991), ribociclib (LEE011), and abemaciclib (LY2835219), which have demonstrated manageable toxicities, improved pharmacokinetic properties, and impressive antitumor activity, especially in certain forms of breast cancer (14, 16). Successful early clinical trials with these three CDK4/6 inhibitors have generated cautious enthusiasm that these drugs may emerge as a new class of anticancer agents (14, 17). Palbociclib recently was approved by Food and Drug Administration for the treatment of metastatic breast cancer and became the first CDK4/6 inhibitor approved for anticancer therapy (18).In addition to its potential as an anticancer strategy, CDK4/6 inhibition in normal tissues could be exploited therapeutically for wide-ranging clinical conditions. For example, radiation-induced myelosuppression, caused by cell death of proliferating hematopoietic stem/progenitor cells, can be rescued by palbociclib (19, 20). Furthermore, cytotoxic anticancer agents cause significant toxicities to normal proliferating cells, which possibly could be mitigated by the concomitant use of CDK4/6 inhibitors (20, 21). More broadly, cell-cycle inhibition could have beneficial effects in disorders in which maladaptive proliferation of normal cells contributes to the disease pathology, as observed in vascular proliferative diseases, hyperproliferative skin diseases, and autoimmune disorders (22, 23). In support of this possibility, palbociclib treatment recently was reported to ameliorate disease progression in animal models of rheumatoid arthritis through cell-cycle inhibition of synovial fibroblasts (24).Abnormal cellular proliferation also is a hallmark of various kidney diseases (25), and cell-cycle inhibition has been shown to ameliorate significantly the pathogenesis of polycystic kidney disease (26), nephritis (27), and acute kidney injury (AKI) (28). Remarkably, during AKI, the normally quiescent renal tubular cells reenter the cell cycle (29–34), and blocking cell-cycle progression can reduce renal injury (28). Here, we provide evidence that the CDK4/6 pathway is activated early during AKI and demonstrate significant protective effects of CDK4/6 inhibitors in animal models of cisplatin-induced AKI. In addition, we found that the CDK4/6 inhibitors palbociclib and LEE011 are potent inhibitors of organic cation transporter 2 (OCT2), a cisplatin uptake transporter highly expressed in renal tubular cells (35–37). Our findings provide a rationale for the clinical development of palbociclib and LEE011 for the prevention and treatment of AKI.     

Exchange protein directly activated by cAMP plays a critical role in bacterial invasion during fatal rickettsioses

Rickettsiae are responsible for some of the most devastating human infections. A high infectivity and severe illness after inhalation make some rickettsiae bioterrorism threats. We report that deletion of the exchange protein directly activated by cAMP (Epac) gene, Epac1, in mice protects them from an ordinarily lethal dose of rickettsiae. Inhibition of Epac1 suppresses bacterial adhesion and invasion. Most importantly, pharmacological inhibition of Epac1 in vivo using an Epac-specific small-molecule inhibitor, ESI-09, completely recapitulates the Epac1 knockout phenotype. ESI-09 treatment dramatically decreases the morbidity and mortality associated with fatal spotted fever rickettsiosis. Our results demonstrate that Epac1-mediated signaling represents a mechanism for host–pathogen interactions and that Epac1 is a potential target for the prevention and treatment of fatal rickettsioses.Rickettsiae are responsible for some of the most devastating human infections (1–4). It has been forecasted that temperature increases attributable to global climate change will lead to more widespread distribution of rickettsioses (5). These tick-borne diseases are caused by obligately intracellular bacteria of the genus Rickettsia, including Rickettsia rickettsii, the causative agent of Rocky Mountain spotted fever (RMSF) in the United States and Latin America (2, 3), and Rickettsia conorii, the causative agent of Mediterranean spotted fever endemic to southern Europe, North Africa, and India (6). A high infectivity and severe illness after inhalation make some rickettsiae (including Rickettsia prowazekii, R. rickettsii, Rickettsia typhi, and R. conorii) bioterrorism threats (7). Although the majority of rickettsial infections can be controlled by appropriate broad-spectrum antibiotic therapy if diagnosed early, up to 20% of misdiagnosed or untreated (1, 3) and 5% of treated RMSF cases (8) result in a fatal outcome caused by acute disseminated vascular endothelial infection and damage (9). Fatality rates as high as 32% have been reported in hospitalized patients diagnosed with Mediterranean spotted fever (10). In addition, strains of R. prowazekii resistant to tetracycline and chloramphenicol have been developed in laboratories (11). Disseminated endothelial infection and endothelial barrier disruption with increased microvascular permeability are the central features of SFG rickettsioses (1, 2, 9). The molecular mechanisms involved in rickettsial infection remain incompletely elucidated (9, 12). A comprehensive understanding of rickettsial pathogenesis and the development of novel mechanism-based treatment are urgently needed.Living organisms use intricate signaling networks for sensing and responding to changes in the external environment. cAMP, a ubiquitous second messenger, is an important molecular switch that translates environmental signals into regulatory effects in cells (13). As such, a number of microbial pathogens have evolved a set of diverse virulence-enhancing strategies that exploit the cAMP-signaling pathways of their hosts (14). The intracellular functions of cAMP are predominantly mediated by the classic cAMP receptor, protein kinase A (PKA), and the more recently discovered exchange protein directly activated by cAMP (Epac) (15). Thus, far, two isoforms, Epac1 and Epac2, have been identified in humans (16, 17). Epac proteins function by responding to increased intracellular cAMP levels and activating the Ras superfamily small GTPases Ras-proximate 1 and 2 (Rap1 and Rap2). Accumulating evidence demonstrates that the cAMP/Epac1 signaling axis plays key regulatory roles in controlling various cellular functions in endothelial cells in vitro, including cell adhesion (18–21), exocytosis (22), tissue plasminogen activator expression (23), suppressor of cytokine signaling 3 (SOCS-3) induction (24–27), microtubule dynamics (28, 29), cell–cell junctions, and permeability and barrier functions (30–37). Considering the critical importance of endothelial cells in rickettsioses, we examined the functional roles of Epac1 in rickettsial pathogenesis in vivo, taking advantage of the recently generated Epac1 knockout mouse (38) and Epac-specific inhibitors (39, 40) generated from our laboratory. Our studies demonstrate that Epac1 plays a key role in rickettsial infection and represents a therapeutic target for fatal rickettsioses.     

Kinetics of nucleotide-dependent structural transitions in the kinesin-1 hydrolysis cycle

To dissect the kinetics of structural transitions underlying the stepping cycle of kinesin-1 at physiological ATP, we used interferometric scattering microscopy to track the position of gold nanoparticles attached to individual motor domains in processively stepping dimers. Labeled heads resided stably at positions 16.4 nm apart, corresponding to a microtubule-bound state, and at a previously unseen intermediate position, corresponding to a tethered state. The chemical transitions underlying these structural transitions were identified by varying nucleotide conditions and carrying out parallel stopped-flow kinetics assays. At saturating ATP, kinesin-1 spends half of each stepping cycle with one head bound, specifying a structural state for each of two rate-limiting transitions. Analysis of stepping kinetics in varying nucleotides shows that ATP binding is required to properly enter the one-head–bound state, and hydrolysis is necessary to exit it at a physiological rate. These transitions differ from the standard model in which ATP binding drives full docking of the flexible neck linker domain of the motor. Thus, this work defines a consensus sequence of mechanochemical transitions that can be used to understand functional diversity across the kinesin superfamily.Kinesin-1 is a motor protein that steps processively toward microtubule plus-ends, tracking single protofilaments and hydrolyzing one ATP molecule per step (1–6). Step sizes corresponding to the tubulin dimer spacing of 8.2 nm are observed when the molecule is labeled by its C-terminal tail (7–10) and to a two-dimer spacing of 16.4 nm when a single motor domain is labeled (4, 11, 12), consistent with the motor walking in a hand-over-hand fashion. Kinesin has served as an important model system for advancing single-molecule techniques (7–10) and is clinically relevant for its role in neurodegenerative diseases (13), making dissection of its step a popular ongoing target of study.Despite decades of work, many essential components of the mechanochemical cycle remain disputed, including (i) how much time kinesin-1 spends in a one-head–bound (1HB) state when stepping at physiological ATP concentrations, (ii) whether the motor waits for ATP in a 1HB or two-heads–bound (2HB) state, and (iii) whether ATP hydrolysis occurs before or after tethered head attachment (4, 11, 14–20). These questions are important because they are fundamental to the mechanism by which kinesins harness nucleotide-dependent structural changes to generate mechanical force in a manner optimized for their specific cellular tasks. Addressing these questions requires characterizing a transient 1HB state in the stepping cycle in which the unattached head is located between successive binding sites on the microtubule. This 1HB intermediate is associated with the force-generating powerstroke of the motor and underlies the detachment pathway that limits motor processivity. Optical trapping (7, 19, 21, 22) and single-molecule tracking studies (4, 8–11) have failed to detect this 1HB state during stepping. Single-molecule fluorescence approaches have detected a 1HB intermediate at limiting ATP concentrations (11, 12, 14, 15), but apart from one study that used autocorrelation analysis to detect a 3-ms intermediate (17), the 1HB state has been undetectable at physiological ATP concentrations.Single-molecule microscopy is a powerful tool for studying the kinetics of structural changes in macromolecules (23). Tracking steps and potential substeps for kinesin-1 at saturating ATP has until now been hampered by the high stepping rates of the motor (up to 100 s⁻¹), which necessitates high frame rates, and the small step size (8.2 nm), which necessitates high spatial precision (7). Here, we apply interferometric scattering microscopy (iSCAT), a recently established single-molecule tool with high spatiotemporal resolution (24–27) to directly visualize the structural changes underlying kinesin stepping. By labeling one motor domain in a dimeric motor, we detect a 1HB intermediate state in which the tethered head resides over the bound head for half the duration of the stepping cycle at saturating ATP. We further show that at physiological stepping rates, ATP binding is required to enter this 1HB state and that ATP hydrolysis is required to exit it. This work leads to a significant revision of the sequence and kinetics of mechanochemical transitions that make up the kinesin-1 stepping cycle and provides a framework for understanding functional diversity across the kinesin superfamily.     

Drosophila seminal protein ovulin mediates ovulation through female octopamine neuronal signaling

Across animal taxa, seminal proteins are important regulators of female reproductive physiology and behavior. However, little is understood about the physiological or molecular mechanisms by which seminal proteins effect these changes. To investigate this topic, we studied the increase in Drosophila melanogaster ovulation behavior induced by mating. Ovulation requires octopamine (OA) signaling from the central nervous system to coordinate an egg’s release from the ovary and its passage into the oviduct. The seminal protein ovulin increases ovulation rates after mating. We tested whether ovulin acts through OA to increase ovulation behavior. Increasing OA neuronal excitability compensated for a lack of ovulin received during mating. Moreover, we identified a mating-dependent relaxation of oviduct musculature, for which ovulin is a necessary and sufficient male contribution. We report further that oviduct muscle relaxation can be induced by activating OA neurons, requires normal metabolic production of OA, and reflects ovulin’s increasing of OA neuronal signaling. Finally, we showed that as a result of ovulin exposure, there is subsequent growth of OA synaptic sites at the oviduct, demonstrating that seminal proteins can contribute to synaptic plasticity. Together, these results demonstrate that ovulin increases ovulation through OA neuronal signaling and, by extension, that seminal proteins can alter reproductive physiology by modulating known female pathways regulating reproduction.Throughout internally fertilizing animals, seminal proteins play important roles in regulating female fertility by altering female physiology and, in some cases, behavior after mating (reviewed in refs. 1–3). Despite this, little is understood about the physiological mechanisms by which seminal proteins induce postmating changes and how their actions are linked with known networks regulating female reproductive physiology.In Drosophila melanogaster, the suite of seminal proteins has been identified, as have many seminal protein-dependent postmating responses, including changes in egg production and laying, remating behavior, locomotion, feeding, and in ovulation rate (reviewed in refs. 2 and 3). For example, the Drosophila seminal protein ovulin elevates ovulation rate to maximal levels during the 24 h following mating (4, 5), and the seminal protein sex peptide (SP) suppresses female mating receptivity and increases egg-laying behavior for several days after mating (6–10). However, although a receptor for SP has been identified (11), along with elements of the neural circuit in which it is required (12–14), SP’s mechanism of action has not yet been linked to regulatory networks known to control postmating behaviors. Thus, a crucial question remains: how do male-derived seminal proteins interact with regulatory networks in females to trigger postmating responses?We addressed this question by examining the stimulation of Drosophila ovulation by the seminal protein ovulin. In insects, ovulation, defined here as the release of an egg from the ovary to the uterus, is among the best understood reproductive processes in terms of its physiology and neurogenetics (15–27). In D. melanogaster, ovulation requires input from neurons in the abdominal ganglia that release the catecholaminergic neuromodulators octopamine (OA) and tyramine (17, 18, 28). Drosophila ovulation also requires an OA receptor, OA receptor in mushroom bodies (OAMB) (19, 20). Moreover, it has been proposed that OA may integrate extrinsic factors to regulate ovulation rates (17). Noradrenaline, the vertebrate structural and functional equivalent to OA (29, 30), is important for mammalian ovulation, and its dysregulation has been associated with ovulation disorders (31–38). In this paper we investigate the role of neurons that release OA and tyramine in ovulin’s action. For simplicity, we refer to these neurons as “OA neurons” to reflect the well-established role of OA in ovulation behavior (16–20, 22).We investigated how action of the seminal protein ovulin relates to the conserved canonical neuromodulatory pathway that regulates ovulation physiology (39–41). We found that ovulin increases ovulation and egg laying through OA neuronal signaling. We also found that ovulin relaxes oviduct muscle tonus, a postmating process that is also mediated by OA neuronal signaling. Finally, subsequent to these effects we detected an ovulin-dependent increase in synaptic sites between OA motor neurons and oviduct muscle, suggesting that ovulin’s stimulation of OA neurons could have increased their synaptic activity. These results suggest that ovulin affects ovulation by manipulating the gain of a neuromodulatory pathway regulating ovulation physiology.     

Predictive Value of Brain Natriuretic Peptides in Patients on Peritoneal Dialysis: Results from the ADEMEX Trial

Background and objectives: Natriuretic peptides have been suggested to be of value in risk stratification in dialysis patients. Data in patients on peritoneal dialysis remain limited.Design, setting, participants, & measurements: Patients of the ADEMEX trial (ADEquacy of peritoneal dialysis in MEXico) were randomized to a control group [standard 4 × 2L continuous ambulatory peritoneal dialysis (CAPD); n = 484] and an intervention group (CAPD with a target creatinine clearance ≥60L/wk/1.73 m²; n = 481). Natriuretic peptides were measured at baseline and correlated with other parameters as well as evaluated for effects on patient outcomes.Results: Control group and intervention group were comparable at baseline with respect to all measured parameters. Baseline values of natriuretic peptides were elevated and correlated significantly with levels of residual renal function but not with body size or diabetes. Baseline values of N-terminal fragment of B-type natriuretic peptide (NT-proBNP) but not proANP(1–30), proANP(31–67), or proANP(1–98) were independently highly predictive of overall survival and cardiovascular mortality. Volume removal was also significantly correlated with patient survival.Conclusions. NT-proBNP have a significant predictive value for survival of CAPD patients and may be of value in guiding risk stratification and potentially targeted therapeutic interventions.Plasma levels of cardiac natriuretic peptides are elevated in patients with chronic kidney disease, owing to impairment of renal function, hypertension, hypervolemia, and/or concomitant heart disease (1–7). Atrial natriuretic peptide (ANP) and particularly brain natriuretic peptide (BNP) levels are linked independently to left ventricular mass (3–5,8–16) and function (3,6–17) and predict total and cardiovascular mortality (1,3,8,10,12,18) as well as cardiac events (12,19). ANP and BNP decrease significantly during hemodialysis treatment but increase again during the interdialytic interval (1,2,4,6,7,14,17,20–23). Levels in patients on peritoneal dialysis (PD) have been found to be lower than in patients on hemodialysis (11,24–26), but the correlations with left ventricular function and structure are maintained in both types of dialysis modalities (11,15,27,28).The high mortality of patients on peritoneal dialysis and the failure of dialytic interventions to alter this mortality (29,30) necessitate renewed attention into novel methods of stratification and identification of patients at highest risk to be targeted for specific interventions. Cardiac natriuretic peptides are increasingly considered to fulfill this role in nonrenal patients. Evaluations of cardiac natriuretic peptides in patients on PD have been limited by small numbers (3,9,11,12,15,24–26) and only one study examined correlations between natriuretic peptide levels and outcomes (12). The PD population enrolled in the ADEMEX trial offered us the opportunity to evaluate cardiac natriuretic peptides and their value in predicting outcomes in the largest clinical trial ever performed on PD (29,30). It is hoped that such an evaluation would identify patients at risk even in the absence of overt clinical disease and hence facilitate or encourage interventions with salutary outcomes.     

Self-assembly of alkynylplatinum(II) terpyridine amphiphiles into nanostructures via steric control and metal–metal interactions

A series of mono- and dinuclear alkynylplatinum(II) terpyridine complexes containing the hydrophilic oligo(para-phenylene ethynylene) with two 3,6,9-trioxadec-1-yloxy chains was designed and synthesized. The mononuclear alkynylplatinum(II) terpyridine complex was found to display a very strong tendency toward the formation of supramolecular structures. Interestingly, additional end-capping with another platinum(II) terpyridine moiety of various steric bulk at the terminal alkyne would lead to the formation of nanotubes or helical ribbons. These desirable nanostructures were found to be governed by the steric bulk on the platinum(II) terpyridine moieties, which modulates the directional metal−metal interactions and controls the formation of nanotubes or helical ribbons. Detailed analysis of temperature-dependent UV-visible absorption spectra of the nanostructured tubular aggregates also provided insights into the assembly mechanism and showed the role of metal−metal interactions in the cooperative supramolecular polymerization of the amphiphilic platinum(II) complexes.Square-planar d⁸ platinum(II) polypyridine complexes have long been known to exhibit intriguing spectroscopic and luminescence properties (1–54) as well as interesting solid-state polymorphism associated with metal−metal and π−π stacking interactions (1–14, 25). Earlier work by our group showed the first example, to our knowledge, of an alkynylplatinum(II) terpyridine system [Pt(tpy)(C ≡ CR)]⁺ that incorporates σ-donating and solubilizing alkynyl ligands together with the formation of Pt···Pt interactions to exhibit notable color changes and luminescence enhancements on solvent composition change (25) and polyelectrolyte addition (26). This approach has provided access to the alkynylplatinum(II) terpyridine and other related cyclometalated platinum(II) complexes, with functionalities that can self-assemble into metallogels (27–31), liquid crystals (32, 33), and other different molecular architectures, such as hairpin conformation (34), helices (35–38), nanostructures (39–45), and molecular tweezers (46, 47), as well as having a wide range of applications in molecular recognition (48–52), biomolecular labeling (48–52), and materials science (53, 54). Recently, metal-containing amphiphiles have also emerged as a building block for supramolecular architectures (42–44, 55–59). Their self-assembly has always been found to yield different molecular architectures with unprecedented complexity through the multiple noncovalent interactions on the introduction of external stimuli (42–44, 55–59).Helical architecture is one of the most exciting self-assembled morphologies because of the uniqueness for the functional and topological properties (60–69). Helical ribbons composed of amphiphiles, such as diacetylenic lipids, glutamates, and peptide-based amphiphiles, are often precursors for the growth of tubular structures on an increase in the width or the merging of the edges of ribbons (64, 65). Recently, the optimization of nanotube formation vs. helical nanostructures has aroused considerable interests and can be achieved through a fine interplay of the influence on the amphiphilic property of molecules (66), choice of counteranions (67, 68), or pH values of the media (69), which would govern the self-assembly of molecules into desirable aggregates of helical ribbons or nanotube scaffolds. However, a precise control of supramolecular morphology between helical ribbons and nanotubes remains challenging, particularly for the polycyclic aromatics in the field of molecular assembly (64–69). Oligo(para-phenylene ethynylene)s (OPEs) with solely π−π stacking interactions are well-recognized to self-assemble into supramolecular system of various nanostructures but rarely result in the formation of tubular scaffolds (70–73). In view of the rich photophysical properties of square-planar d⁸ platinum(II) systems and their propensity toward formation of directional Pt···Pt interactions in distinctive morphologies (27–31, 39–45), it is anticipated that such directional and noncovalent metal−metal interactions might be capable of directing or dictating molecular ordering and alignment to give desirable nanostructures of helical ribbons or nanotubes in a precise and controllable manner.Herein, we report the design and synthesis of mono- and dinuclear alkynylplatinum(II) terpyridine complexes containing hydrophilic OPEs with two 3,6,9-trioxadec-1-yloxy chains. The mononuclear alkynylplatinum(II) terpyridine complex with amphiphilic property is found to show a strong tendency toward the formation of supramolecular structures on diffusion of diethyl ether in dichloromethane or dimethyl sulfoxide (DMSO) solution. Interestingly, additional end-capping with another platinum(II) terpyridine moiety of various steric bulk at the terminal alkyne would result in nanotubes or helical ribbons in the self-assembly process. To the best of our knowledge, this finding represents the first example of the utilization of the steric bulk of the moieties, which modulates the formation of directional metal−metal interactions to precisely control the formation of nanotubes or helical ribbons in the self-assembly process. Application of the nucleation–elongation model into this assembly process by UV-visible (UV-vis) absorption spectroscopic studies has elucidated the nature of the molecular self-assembly, and more importantly, it has revealed the role of metal−metal interactions in the formation of these two types of nanostructures.     

Understanding the kinetic mechanism of RNA single base pair formation

RNA functions are intrinsically tied to folding kinetics. The most elementary step in RNA folding is the closing and opening of a base pair. Understanding this elementary rate process is the basis for RNA folding kinetics studies. Previous studies mostly focused on the unfolding of base pairs. Here, based on a hybrid approach, we investigate the folding process at level of single base pairing/stacking. The study, which integrates molecular dynamics simulation, kinetic Monte Carlo simulation, and master equation methods, uncovers two alternative dominant pathways: Starting from the unfolded state, the nucleotide backbone first folds to the native conformation, followed by subsequent adjustment of the base conformation. During the base conformational rearrangement, the backbone either retains the native conformation or switches to nonnative conformations in order to lower the kinetic barrier for base rearrangement. The method enables quantification of kinetic partitioning among the different pathways. Moreover, the simulation reveals several intriguing ion binding/dissociation signatures for the conformational changes. Our approach may be useful for developing a base pair opening/closing rate model.RNAs perform critical cellular functions at the level of gene expression and regulation (1–4). RNA functions are determined not only by RNA structure or structure motifs [e.g., tetraloop hairpins (5, 6)] but also by conformational distributions and dynamics and kinetics of conformational changes. For example, riboswitches can adopt different conformations in response to specific conditions of the cellular environment (7, 8). Understanding the kinetics, such as the rate and pathways for the conformational changes, is critical for deciphering the mechanism of RNA function (9–19). Extensive experimental and theoretical studies on RNA folding kinetics have provided significant insights into the kinetic mechanism of RNA functions (19–36). However, due to the complexity of the RNA folding energy landscape (37–46) and the limitations of experimental tools (47–55), many fundamental problems, including single base flipping and base pair formation and fraying, remain unresolved. These unsolved fundamental problems have hampered our ability to resolve other important issues, such as RNA hairpin and larger structure folding kinetics. Several key questions remain unanswered, such as whether the hairpin folding is rate-limited by the conformational search of the native base pairs, whose formation leads to fast downhill folding of the whole structure, or by the breaking of misfolded base pairs before refolding to the native structure (18, 19, 54–73).Motivated by the need to understand the basic steps of nucleic acids folding, Hagan et al. (74) performed forty-three 200-ps unfolding trajectories at 400 K and identified both on- and off-pathway intermediates and two dominant unfolding pathways for a terminal C-G base pair in a DNA duplex. In one of the pathways, base pairing and stacking interactions are broken concomitantly, whereas in the other pathway, base stacking is broken after base pairing is disrupted. Furthermore, the unfolding requires that the Cyt diffuse away from the pairing Gua to a distance such that the C-G hydrogen bond cannot reform easily. More recently, Colizzi and Bussi (75) performed molecular dynamics (MD) pulling simulations for an RNA duplex and construct free energy landscape from the pulling simulation. The simulation showed that the base pair opening reaction starts with the unbinding of the 5′-base, followed by the unbinding of the 3′-base (i.e., the 5′-base is less stable than the 3′-base). These previous unfolding simulations offered significant insights into the pathways and transition states. However, as shown below, several important issues remain.One intriguing problem is the rate model for base pairing. There are currently three main types of models. In the first type of model, the barrier ΔG+‡ for closing a base pair is dominated by the entropic cost ΔS for positioning the nucleotides to the base-paired configuration and the barrier ΔG−‡ for opening a base pair is the enthalpic cost ΔH for disrupting the hydrogen bonds and base stacking interactions (18, 59, 60). In the second type of model, ΔG+‡ is the net free energy change for base pairing ΔG = ΔH ? TΔS and ΔG−‡ is zero (76, 77). In the third type of model, ΔG±‡=±ΔG/2 is used (78). In addition to the above three main types, other models, such as more sophisticated hybrid rate models, have been proposed (29).In this paper, we report a hybrid method (see Fig. 1) to investigate the single base pairing process. In contrast to the previous simulations for temperature- or force-induced unfolding reactions, we directly model the folding process here (i.e., the base pair closing process). Specifically, we use MD simulations to identify the conformational clusters. Based on the network of the conformational clusters as a reduced conformational ensemble, we apply kinetic Monte Carlo (KMC) and master equation (ME) methods to elucidate the detailed roles of base pairing and stacking interactions, as well as the roles of water and ions (79–82). The study reveals previously unidentified kinetics pathways, misfolded states, and rate-limiting steps. A clear understanding of the microscopic details of the elementary kinetic move is a prerequisite for further rigorous study of large-scale RNA kinetic studies. The method described here may provide a feasible way to develop a rate model for the base pair/stack-based kinetic move set. Furthermore, the mechanism of RNA single base folding may provide useful insights into many biologically significant processes, such as nucleotide flipping (83) in helicases and base pair fraying (84) (as the possible first step for nucleic duplex melting in nucleic acid enzymatic processes).Open in a separate windowFig. 1.(A) Folding of a single nucleotide (G1, red) from the unfolded (Left) to the native folded (Right) state. (B) Exhaustive sampling for the (discrete) conformations of the G1 nucleotide (Right) through enumeration of the torsion angles (formed by the blue bonds). (C) Schematic plot shows the trajectories on the energy landscape (depicted with two reaction coordinates for clarity) explored by the MD simulations. The lines, open circles, and hexagons denote the trajectories; the initial states; and the (centroid structures of the) clusters, respectively. (D) Conformational network based on six clusters. (E) The rmsds to the different clusters provide information about the structural changes in a MD trajectory.