Integrating cellular electron microscopy with multimodal data to explore biology across space and time.

Biology spans a continuum of length and time scales. Individual experimental methods only glimpse discrete pieces of this spectrum but can be combined to construct a more holistic view. In this Review, we detail the latest advancements in volume electron microscopy (vEM) and cryo-electron tomography (cryo-ET), which together can visualize biological complexity across scales from the organization of cells in large tissues to the molecular details inside native cellular environments. In addition, we discuss emerging methodologies for integrating three-dimensional electron microscopy (3DEM) imaging with multimodal data, including fluorescence microscopy, mass spectrometry, single-particle analysis, and AI-based structure prediction. This multifaceted approach fills gaps in the biological continuum, providing functional context, spatial organization, molecular identity, and native interactions. We conclude with a perspective on incorporating diverse data into computational simulations that further bridge and extend length scales while integrating the dimension of time.

Architecture and self-assembly of the jumbo bacteriophage nuclear shell.

Bacteria encode myriad defences that target the genomes of infecting bacteriophage, including restriction-modification and CRISPR-Cas systems1. In response, one family of large bacteriophages uses a nucleus-like compartment to protect its replicating genomes by excluding host defence factors2-4. However, the principal composition and structure of this compartment remain unknown. Here we find that the bacteriophage nuclear shell assembles primarily from one protein, which we name chimallin (ChmA). Combining cryo-electron tomography of nuclear shells in bacteriophage-infected cells and cryo-electron microscopy of a minimal chimallin compartment in vitro, we show that chimallin self-assembles as a flexible sheet into closed micrometre-scale compartments. The architecture and assembly dynamics of the chimallin shell suggest mechanisms for its nucleation and growth, and its role as a scaffold for phage-encoded factors mediating macromolecular transport, cytoskeletal interactions, and viral maturation.

Targeted protein S-nitrosylation of ACE2 inhibits SARS-CoV-2 infection.

Prevention of infection and propagation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a high priority in the Coronavirus Disease 2019 (COVID-19) pandemic. Here we describe S-nitrosylation of multiple proteins involved in SARS-CoV-2 infection, including angiotensin-converting enzyme 2 (ACE2), the receptor for viral entry. This reaction prevents binding of ACE2 to the SARS-CoV-2 spike protein, thereby inhibiting viral entry, infectivity and cytotoxicity. Aminoadamantane compounds also inhibit coronavirus ion channels formed by envelope (E) protein. Accordingly, we developed dual-mechanism aminoadamantane nitrate compounds that inhibit viral entry and, thus, the spread of infection by S-nitrosylating ACE2 via targeted delivery of the drug after E protein channel blockade. These non-toxic compounds are active in vitro and in vivo in the Syrian hamster COVID-19 model and, thus, provide a novel avenue to pursue therapy.

Structure and dynamics of SARS-CoV-2 proofreading exoribonuclease ExoN.

High-fidelity replication of the large RNA genome of coronaviruses (CoVs) is mediated by a 3'-to-5' exoribonuclease (ExoN) in nonstructural protein 14 (nsp14), which excises nucleotides including antiviral drugs misincorporated by the low-fidelity viral RNA-dependent RNA polymerase (RdRp) and has also been implicated in viral RNA recombination and resistance to innate immunity. Here, we determined a 1.6-Å resolution crystal structure of severe acute respiratory syndrome CoV 2 (SARS-CoV-2) ExoN in complex with its essential cofactor, nsp10. The structure shows a highly basic and concave surface flanking the active site, comprising several Lys residues of nsp14 and the N-terminal amino group of nsp10. Modeling suggests that this basic patch binds to the template strand of double-stranded RNA substrates to position the 3' end of the nascent strand in the ExoN active site, which is corroborated by mutational and computational analyses. We also show that the ExoN activity can rescue a stalled RNA primer poisoned with sofosbuvir and allow RdRp to continue its extension in the presence of the chain-terminating drug, biochemically recapitulating proofreading in SARS-CoV-2 replication. Molecular dynamics simulations further show remarkable flexibility of multidomain nsp14 and suggest that nsp10 stabilizes ExoN for substrate RNA binding to support its exonuclease activity. Our high-resolution structure of the SARS-CoV-2 ExoN-nsp10 complex serves as a platform for future development of anticoronaviral drugs or strategies to attenuate the viral virulence.

Electrophysical cardiac remodeling at the molecular level: Insights into ryanodine receptor activation and calcium-induced calcium release from a stochastic explicit-particle model.

We present the first-ever, fully discrete, stochastic model of triggered cardiac calcium dynamics. Using anatomically accurate subcellular cardiac myocyte geometries, we simulate the molecular players involved in calcium handling using high-resolution stochastic and explicit-particle methods at the level of an individual cardiac dyadic junction. Integrating data from multiple experimental sources, the model not only replicates the findings of traditional in silico studies and complements in vitro experimental data but also reveals new insights into the molecular mechanisms driving cardiac dysfunction under stress and disease conditions. We improve upon older, nondiscrete models using the same realistic geometry by incorporating molecular mechanisms for spontaneous, as well as triggered calcium-induced calcium release (CICR). Action potentials are used to activate L-type calcium channels, triggering CICR through ryanodine receptors (RyRs) on the surface of the sarcoplasmic reticulum. These improvements allow for the specific focus on the couplon: the structure-function relationship between L-type calcium channel and RyR. We investigate the electrophysical effects of normal and diseased action potentials on CICR and interrogate the effects of dyadic junction deformation through detubulation and orphaning of RyR. Our work demonstrates the importance of the electrophysical integrity of the calcium release unit on CICR fidelity, giving insights into the molecular basis of heart disease. Finally, we provide a unique, detailed, molecular view of the CICR process using advanced rendering techniques. This easy-to-use model comes complete with tutorials and all necessary software for use and analysis so as to maximize usability and reproducibility. Our work focuses on quantifying, qualifying, and visualizing the behavior of the molecular species that underlie the function and dysfunction of subcellular cardiomyocyte systems.

Binding and Activation of LRP1-Dependent Cell Signaling in Schwann Cells Using a Peptide Derived from the Hemopexin Domain of MMP-9.

Schwann cells (SCs) undergo phenotypic transformation and then orchestrate nerve repair following a peripheral nervous system injury. The low-density lipoprotein receptor-related protein-1 (LRP1) is significantly upregulated in SCs in response to acute injury, activating cJun and promoting SC survival. Matrix-metalloproteinase-9 (MMP-9) is an LRP1 ligand that binds LRP1 through its hemopexin domain (PEX) and activates SC survival signaling and migration. To identify novel peptide mimetics within the hemopexin domain of MMP-9, we examined the crystal structure of PEX, synthesized four peptides, and examined their potential to bind and activate LRP1. We demonstrate that a 22 amino acid peptide, peptide 2, was the only peptide that activated Akt and ERK1/2 signaling in SCs, similar to a glutathione s-transferase (GST)-fused holoprotein, GST-PEX. Intraneural injection of peptide 2, but not vehicle, into crush-injured sciatic nerves activated cJun greater than 2.5-fold in wild-type mice, supporting that peptide 2 can activate the SC repair signaling in vivo. Peptide 2 also bound to Fc-fusion proteins containing the ligand-binding motifs of LRP1, clusters of complement-like repeats (CCRII and CCRIV). Pulldown and computational studies of alanine mutants of peptide 2 showed that positively charged lysine and arginine amino acids within the peptide are critical for stability and binding to CCRII. Collectively, these studies demonstrate that a novel peptide derived from PEX can serve as an LRP1 agonist and possesses qualities previously associated with LRP1 binding and SC signaling in vitro and in vivo.

Breaking Down the Bottlebrush: Atomically Detailed Structural Dynamics of Mucins.

Mucins, the biomolecular components of mucus, are glycoproteins that form a thick physical barrier at all tissue-air interfaces, forming a first line of defense against pathogens. Structural features of mucins and their interactions with other biomolecules remain largely unexplored due to the challenges associated with their high-resolution characterization. Combining limited mass spectrometry glycomics and protein sequencing data, we present all-atom, explicitly solvated molecular dynamics simulations of a major respiratory mucin, MUC5B. We detail key forces and degrees of freedom imposed by the extensive O-glycosylation, which imbue the canonically observed bottlebrush-like structures to these otherwise intrinsically disordered protein backbones. We compare our simulation results to static structures observed in recent scanning tunneling microscopy experiments as well as other published experimental efforts. Our work represents the demonstration of a workflow applied to a mucin example, which we hope will be employed by other groups to investigate the dynamics and interactions of other mucins, which can inform on structural details currently inaccessible to experimental techniques.

NetSci: A Library for High Performance Biomolecular Simulation Network Analysis Computation.

We present the NetSci program-an open-source scientific software package designed for estimating mutual information (MI) between data sets using GPU acceleration and a k-nearest-neighbor algorithm. This approach significantly enhances calculation speed, achieving improvements of several orders of magnitude over traditional CPU-based methods, with data set size limits dictated only by available hardware. To validate NetSci, we accurately compute MI for an analytically verifiable two-dimensional Gaussian distribution and replicate the generalized correlation (GC) analysis previously conducted on the B1 domain of protein G. We also apply NetSci to molecular dynamics simulations of the Sarcoendoplasmic Reticulum Calcium-ATPase (SERCA) pump, exploring the allosteric mechanisms and pathways influenced by ATP and 2'-deoxy-ATP (dATP) binding. Our analysis reveals distinct allosteric effects induced by ATP compared to dATP, with predicted information pathways from the bound nucleotide to the calcium-binding domain differing based on the nucleotide involved. NetSci proves to be a valuable tool for estimating MI and GC in various data sets and is particularly effective for analyzing intraprotein communication and information transfer.

Mining for Potent Inhibitors through Artificial Intelligence and Physics: A Unified Methodology for Ligand Based and Structure Based Drug Design.

Determining the viability of a new drug molecule is a time- and resource-intensive task that makes computer-aided assessments a vital approach to rapid drug discovery. Here we develop a machine learning algorithm, iMiner, that generates novel inhibitor molecules for target proteins by combining deep reinforcement learning with real-time 3D molecular docking using AutoDock Vina, thereby simultaneously creating chemical novelty while constraining molecules for shape and molecular compatibility with target active sites. Moreover, through the use of various types of reward functions, we have introduced novelty in generative tasks for new molecules such as chemical similarity to a target ligand, molecules grown from known protein bound fragments, and creation of molecules that enforce interactions with target residues in the protein active site. The iMiner algorithm is embedded in a composite workflow that filters out Pan-assay interference compounds, Lipinski rule violations, uncommon structures in medicinal chemistry, and poor synthetic accessibility with options for cross-validation against other docking scoring functions and automation of a molecular dynamics simulation to measure pose stability. We also allow users to define a set of rules for the structures they would like to exclude during the training process and postfiltering steps. Because our approach relies only on the structure of the target protein, iMiner can be easily adapted for the future development of other inhibitors or small molecule therapeutics of any target protein.

Interfacial Enrichment of Lauric Acid Assisted by Long-Chain Fatty Acids, Acidity and Salinity at Sea Spray Aerosol Surfaces.

Surfactant monolayers at sea spray aerosol (SSA) surfaces regulate various atmospheric processes including gas transfer, cloud interactions, and radiative properties. Most experimental studies of SSA employ a simplified surfactant mixture of long-chain fatty acids (LCFAs) as a proxy for the sea surface microlayer or SSA surface. However, medium-chain fatty acids (MCFAs) make up nearly 30% of the FA fraction in nascent SSA. Given that LCFA monolayers are easily disrupted upon the introduction of chemical heterogeneity (such as mixed chain lengths), simple FA proxies are unlikely to represent realistic SSA interfaces. Integrating experimental and computational techniques, we characterize the impact that partially soluble MCFAs have on the properties of atmospherically relevant LCFA mixtures. We explore the extent to which the MCFA lauric acid (LA) is surface stabilized by varying acidity, salinity, and monolayer composition. We also discuss the impacts of pH on LCFA-assisted LA retention, where the presence of LCFAs may shift the surface-adsorption equilibria of laurate─the conjugate base─toward higher surface activities. Molecular dynamic simulations suggest a mechanism for the enhanced surface retention of laurate. We conclude that increased FA heterogeneity at SSA surfaces promotes surface activity of soluble FA species, altering monolayer phase behavior and impacting climate-relevant atmospheric processes.

Independent Markov decomposition: Toward modeling kinetics of biomolecular complexes.

To advance the mission of in silico cell biology, modeling the interactions of large and complex biological systems becomes increasingly relevant. The combination of molecular dynamics (MD) simulations and Markov state models (MSMs) has enabled the construction of simplified models of molecular kinetics on long timescales. Despite its success, this approach is inherently limited by the size of the molecular system. With increasing size of macromolecular complexes, the number of independent or weakly coupled subsystems increases, and the number of global system states increases exponentially, making the sampling of all distinct global states unfeasible. In this work, we present a technique called independent Markov decomposition (IMD) that leverages weak coupling between subsystems to compute a global kinetic model without requiring the sampling of all combinatorial states of subsystems. We give a theoretical basis for IMD and propose an approach for finding and validating such a decomposition. Using empirical few-state MSMs of ion channel models that are well established in electrophysiology, we demonstrate that IMD models can reproduce experimental conductance measurements with a major reduction in sampling compared with a standard MSM approach. We further show how to find the optimal partition of all-atom protein simulations into weakly coupled subunits.

SARS-CoV-2 escape from a highly neutralizing COVID-19 convalescent plasma.

To investigate the evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the immune population, we coincupi bated the authentic virus with a highly neutralizing plasma from a COVID-19 convalescent patient. The plasma fully neutralized the virus for seven passages, but, after 45 d, the deletion of F140 in the spike N-terminal domain (NTD) N3 loop led to partial breakthrough. At day 73, an E484K substitution in the receptor-binding domain (RBD) occurred, followed, at day 80, by an insertion in the NTD N5 loop containing a new glycan sequon, which generated a variant completely resistant to plasma neutralization. Computational modeling predicts that the deletion and insertion in loops N3 and N5 prevent binding of neutralizing antibodies. The recent emergence in the United Kingdom, South Africa, Brazil, and Japan of natural variants with similar changes suggests that SARS-CoV-2 has the potential to escape an effective immune response and that vaccines and antibodies able to control emerging variants should be developed.

Selectivity and Ranking of Tight-Binding JAK-STAT Inhibitors Using Markovian Milestoning with Voronoi Tessellations.

Janus kinases (JAK), a group of proteins in the nonreceptor tyrosine kinase (NRTKs) family, play a crucial role in growth, survival, and angiogenesis. They are activated by cytokines through the Janus kinase-signal transducer and activator of a transcription (JAK-STAT) signaling pathway. JAK-STAT signaling pathways have significant roles in the regulation of cell division, apoptosis, and immunity. Identification of the V617F mutation in the Janus homology 2 (JH2) domain of JAK2 leading to myeloproliferative disorders has stimulated great interest in the drug discovery community to develop JAK2-specific inhibitors. However, such inhibitors should be selective toward JAK2 over other JAKs and display an extended residence time. Recently, novel JAK2/STAT5 axis inhibitors (N-(1H-pyrazol-3-yl)pyrimidin-2-amino derivatives) have displayed extended residence times (hours or longer) on target and adequate selectivity excluding JAK3. To facilitate a deeper understanding of the kinase-inhibitor interactions and advance the development of such inhibitors, we utilize a multiscale Markovian milestoning with Voronoi tessellations (MMVT) approach within the Simulation-Enabled Estimation of Kinetic Rates v.2 (SEEKR2) program to rank order these inhibitors based on their kinetic properties and further explain the selectivity of JAK2 inhibitors over JAK3. Our approach investigates the kinetic and thermodynamic properties of JAK-inhibitor complexes in a user-friendly, fast, efficient, and accurate manner compared to other brute force and hybrid-enhanced sampling approaches.

SEEKR2: Versatile Multiscale Milestoning Utilizing the OpenMM Molecular Dynamics Engine.

We present SEEKR2 (simulation-enabled estimation of kinetic rates version 2)─the latest iteration in the family of SEEKR programs for using multiscale simulation methods to computationally estimate the kinetics and thermodynamics of molecular processes, in particular, ligand-receptor binding. SEEKR2 generates equivalent, or improved, results compared to the earlier versions of SEEKR but with significant increases in speed and capabilities. SEEKR2 has also been built with greater ease of usability and with extensible features to enable future expansions of the method. Now, in addition to supporting simulations using NAMD, calculations may be run with the fast and extensible OpenMM simulation engine. The Brownian dynamics portion of the calculation has also been upgraded to Browndye 2. Furthermore, this version of SEEKR supports hydrogen mass repartitioning, which significantly reduces computational cost, while showing little, if any, loss of accuracy in the predicted kinetics.

Benchmarking ensemble docking methods in D3R Grand Challenge 4.

The discovery of new drugs is a time consuming and expensive process. Methods such as virtual screening, which can filter out ineffective compounds from drug libraries prior to expensive experimental study, have become popular research topics. As the computational drug discovery community has grown, in order to benchmark the various advances in methodology, organizations such as the Drug Design Data Resource have begun hosting blinded grand challenges seeking to identify the best methods for ligand pose-prediction, ligand affinity ranking, and free energy calculations. Such open challenges offer a unique opportunity for researchers to partner with junior students (e.g., high school and undergraduate) to validate basic yet fundamental hypotheses considered to be uninteresting to domain experts. Here, we, a group of high school-aged students and their mentors, present the results of our participation in Grand Challenge 4 where we predicted ligand affinity rankings for the Cathepsin S protease, an important protein target for autoimmune diseases. To investigate the effect of incorporating receptor dynamics on ligand affinity rankings, we employed the Relaxed Complex Scheme, a molecular docking method paired with molecular dynamics-generated receptor conformations. We found that Cathepsin S is a difficult target for molecular docking and we explore some advanced methods such as distance-restrained docking to try to improve the correlation with experiments. This project has exemplified the capabilities of high school students when supported with a rigorous curriculum, and demonstrates the value of community-driven competitions for beginners in computational drug discovery.

Structural basis for ligand modulation of the CCR2 conformational landscape.

CC chemokine receptor 2 (CCR2) is a part of the chemokine receptor family, an important class of therapeutic targets. These class A G-protein coupled receptors (GPCRs) are involved in mammalian signaling pathways and control cell migration toward endogenous CC chemokine ligands, named for the adjacent cysteine motif on their N terminus. Chemokine receptors and their associated ligands are involved in a wide range of diseases and thus have become important drug targets. CCR2, in particular, promotes the metastasis of cancer cells and is also implicated in autoimmunity-driven type-1 diabetes, diabetic nephropathy, multiple sclerosis, asthma, atherosclerosis, neuropathic pain, and rheumatoid arthritis. Although promising, CCR2 antagonists have been largely unsuccessful to date. Here, we investigate the effect of an orthosteric and an allosteric antagonist on CCR2 dynamics by coupling long-timescale molecular dynamics simulations with Markov-state model theory. We find that the antagonists shift CCR2 into several stable inactive conformations that are distinct from the crystal structure conformation and disrupt a continuous internal water and sodium ion pathway, preventing transitions to an active-like state. Several metastable conformations present a cryptic drug-binding pocket near the allosteric site that may be amenable to targeting with small molecules. Without antagonists, the apo dynamics reveal intermediate conformations along the activation pathway that provide insight into the basal dynamics of CCR2 and may also be useful for future drug design.

The flexibility of ACE2 in the context of SARS-CoV-2 infection.

The coronavirus disease 2019 (COVID-19) pandemic has swept over the world in the past months, causing significant loss of life and consequences to human health. Although numerous drug and vaccine development efforts are underway, there are many outstanding questions on the mechanism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral association to angiotensin-converting enzyme 2 (ACE2), its main host receptor, and host cell entry. Structural and biophysical studies indicate some degree of flexibility in the viral extracellular spike glycoprotein and at the receptor-binding domain (RBD)-receptor interface, suggesting a role in infection. Here, we perform explicitly solvated, all-atom, molecular dynamics simulations of the glycosylated, full-length, membrane-bound ACE2 receptor in both an apo and spike RBD-bound state to probe the intrinsic dynamics of the ACE2 receptor in the context of the cell surface. A large degree of fluctuation in the full-length structure is observed, indicating hinge bending motions at the linker region connecting the head to the transmembrane helix while still not disrupting the ACE2 homodimer or ACE2-RBD interfaces. This flexibility translates into an ensemble of ACE2 homodimer conformations that could sterically accommodate binding of the spike trimer to more than one ACE2 homodimer and suggests a mechanical contribution of the host receptor toward the large spike conformational changes required for cell fusion. This work presents further structural and functional insights into the role of ACE2 in viral infection that can potentially be exploited for the rational design of effective SARS-CoV-2 therapeutics.

A multiscale coarse-grained model of the SARS-CoV-2 virion.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and ongoing development of a largely "bottom-up" coarse-grained (CG) model of the SARS-CoV-2 virion. Data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data become publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion.

A potential interaction between the SARS-CoV-2 spike protein and nicotinic acetylcholine receptors.

Changeux et al. (Changeux et al. C. R. Biol. 343:33-39.) recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs) and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs and of the spike to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4ß2 and α7 subtypes and the muscle-like αßγδ receptor from Tetronarce californica. Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behavior of the bound Y674-R685 is highly dependent on the receptor subtype; it adopts extended conformations in the α4ß2 and α7 complexes but is more compact when bound to the muscle-like receptor. In the α4ß2 and αßγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket in which it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1, and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a nonparallel arrangement to one another.

OCRE Domains of Splicing Factors RBM5 and RBM10: Tyrosine Ring-Flip Frequencies Determined by Integrated Use of 1 H NMR Spectroscopy and Molecular Dynamics Simulations.

The 55-residue OCRE domains of the splicing factors RBM5 and RBM10 contain 15 tyrosines in compact, globular folds. At 25 °C, all 15 tyrosines show symmetric 1 H NMR spectra, with averaged signals for the pairs of δ- and ϵ-ring hydrogens. At 4 °C, two tyrosines were identified as showing 1 H NMR line-broadening due to lowered frequency of the ring-flipping. For the other 13 tyrosine rings, it was not evident, from the 1 H NMR data alone, whether they were either all flipping at high frequencies, or whether slowed flipping went undetected due to small chemical-shift differences between pairs of exchanging ring hydrogen atoms. Here, we integrate 1 H NMR spectroscopy and molecular dynamics (MD) simulations to determine the tyrosine ring-flip frequencies. In the RBM10-OCRE domain, we found that, for 11 of the 15 tyrosines, these frequencies are in the range 2.0×106 to 1.3×108  s-1 , and we established an upper limit of <1.0×106  s-1 for the remaining four residues. The experimental data and the MD simulation are mutually supportive, and their combined use extends the analysis of aromatic ring-flip events beyond the limitations of routine 1 H NMR line-shape analysis into the nanosecond frequency range.