Human cardiac ß-myosin powerstroke energetics: Thin filament, Pi displacement, and mutation effects.

The powerstroke of human cardiac ß-myosin is an important stage of the cross-bridge cycle that generates force for muscle contraction. However, the starting structure of this process has never been resolved, and the relative timing of the powerstroke and inorganic phosphate (Pi) release is still controversial. In this study, we generated an atomistic model of myosin on the thin filament and utilized metadynamics simulations to predict the absent starting structure of the powerstroke. We demonstrated that the displacement of Pi from the active site during the powerstroke is likely necessary, reducing the energy barrier of the conformation change. The effects of the presence of the thin filament, the hypertrophic cardiomyopathy mutation R712L, and the binding of mavacamten on the powerstroke process were also investigated.

Myosin-Catalyzed ATP Hydrolysis in the Presence of Disease-Causing Mutations: Mavacamten as a Way to Repair Mechanism.

Hypertrophic cardiomyopathy is one of the most common forms of genetic cardiomyopathy. Mavacamten is a first-in-class myosin modulator that was identified via activity screening on the wild type, and it is FDA-approved for the treatment of obstructive hypertrophic cardiomyopathy (HCM). The drug selectively binds to the cardiac ß-myosin, inhibiting myosin function to decrease cardiac contractility. Though the drug is thought to affect multiple steps of the myosin cross-bridge cycle, its detailed mechanism of action is still under investigation. Individual steps in the overall cross-bridge cycle must be queried to elucidate the full mechanism of action. In this study, we utilize the rare-event method of transition path sampling to generate reactive trajectories to gain insights into the action of the drug on the dynamics and rate of the ATP hydrolysis step for human cardiac ß-myosin. We study three known HCM causative myosin mutations: R453C, P710R, and R712L to observe the effect of the drug on the alterations caused by these mutations in the chemical step. Since the crystal structure of the drug-bound myosin was not available at the time of this work, we created a model of the drug-bound system utilizing a molecular docking approach. We find a significant effect of the drug in one case, where the actual mechanism of the reaction is altered from the wild type by mutation. The drug restores both the rate of hydrolysis to the wildtype level and the mechanism of the reaction. This is a way to check the effect of the drug on untested mutations.

Atomistic description of the relationship between protein dynamics and catalysis with transition path sampling.

Despite initial resistance, it has been increasingly accepted that protein dynamics plays a role in enzymatic catalysis. There have been two lines of research. Some works study slow conformational motions that are not coupled to the reaction coordinate, but guide the system towards catalytically competent conformations. Understanding at the atomistic level how this is accomplished has remained elusive except for a few systems. In this review we focus on fast sub-picosecond motions that are coupled to the reaction coordinate. The use of Transition Path Sampling has allowed us an atomistic description of how these rate-promoting vibrational motions are incorporated in the reaction mechanism. We will also show how we used insights from rate-promoting motions in protein design.

Protein Dynamics and Enzymatic Catalysis.

This Perspective presents a review of our work and that of others in the highly controversial topic of the coupling of protein dynamics to reaction in enzymes. We have been involved in studying this topic for many years. Thus, this perspective will naturally present our own views, but it also is designed to present an overview of the variety of viewpoints of this topic, both experimental and theoretical. This is obviously a large and contentious topic.

Frosted Branch Angiitis in the Setting of Active COVID-19 Infection and Underlying Mixed Connective Tissue Disease.

Frosted branch angiitis (FBA) is an uncommon form of retinal vasculitis and is typically associated with vision loss. We report a unique case of FBA that manifested in the setting of an active COVID-19 infection in a patient with Mixed Connective Tissue Disease (MCTD). A 34-year-old female with a history of MCTD, including overlapping findings of dermatomyositis, systemic lupus erythematosus, and rheumatoid arthritis, on immunosuppressive medications, presented for left-sided vision loss. She was also found to have an active COVID-19 infection with symptoms including sore throat and dry cough. The patient's visual acuity was counting fingers in her affected eye with a fundus exam revealing diffuse retinal hemorrhages, retinal whitening, cystoid macular edema, and perivascular sheathing of tertiary arterioles and venules, characteristic of FBA. Labs showed mildly elevated inflammatory markers. She exhibited no other signs or symptoms concerning systemic rheumatologic flare. There was no evidence of COVID-19 on viral PCR testing of intraocular fluid but given her positive nasopharyngeal PCR, COVID-induced retinal vasculitis with FBA remained high on the differential. The patient's retinal vasculitis later improved with heightened immunosuppressive therapy including high-dose intravenous corticosteroids. Clinicians should be aware of the possibility of COVID-related FBA, particularly in patients with an underlying predisposition to autoimmune inflammation. Our experience with this patient highlights the utility of high-dose systemic immunosuppressive therapy in treating such inflammatory occlusive retinal vasculitis. Further studies are needed to characterize retinal manifestations of COVID-19 in the setting of autoimmune disease.

Classical Molecular Dynamics Simulation of Glyonic Liquids: Structural Insights and Relation to Conductive Properties.

Rhamnolipids are biosurfactants that have obtained wide industrial and environmental interests with their biodegradability and great surface activity. Besides their important roles as surfactants, they are found to function as a new type of glycolipid-based protic ionic liquids (ILs)─glyonic liquids (GLs). GLs are reported to have impressive physicochemical properties, especially superionic conductivity, and it was reported in experiments that specific ion selections and the fraction of water content have a strong effect on the conductivity. Also, the shape of the conductivity curve as a function of water fraction in GLs is interesting with a sharp increase first and a long plateau. We related the conductivities to the three-dimensional (3D) networks composed of -OH inside the GLs utilizing classical molecular dynamics (MD) simulations. The amount and size of these networks vary with both ion species and water fractions. Before reaching the first hydration layer, the -OH networks with higher projection/box length ratios indicate better conductivity; after reaching the first hydration layer and forming continuous structures, the conductivity retains with more water molecules participating in the continuous networks. Therefore, networks are found to be a qualitative predictor of actual conductivity. This is explained by the analysis of the atomic structures, including radial distribution function, fraction free volume, anion conformations, and hydrogen bond occupancies, of GLs and their water mixtures under different chemical conditions.

Connecting Conformational Motions to Rapid Dynamics in Human Purine Nucleoside Phosphorylase.

The influence of protein motions on enzyme catalysis remains a topic of active discussion. Protein motions occur across a variety of time scales, from vibrational fluctuations in femtoseconds, to collective motions in milliseconds. There have been numerous studies that show conformational motions may assist in catalysis, protein folding, and substrate specificity. It is also known through transition path sampling studies that rapid promoting vibrations contribute to enzyme catalysis. Human purine nucleoside phosphorylase (PNP) is one enzyme that contains both an important conformational motion and a rapid promoting vibration. The slower motion in this enzyme is associated with a loop motion, that when open allows substrate entry and product release but closes over the active site during catalysis. We examine the differences between an unconstrained PNP structure and a PNP structure with constraints on the loop motion. To investigate possible coupling between the slow and fast protein dynamics, we employed transition path sampling, reaction coordinate identification, electric field calculations, and free energy calculations reported here.

Insights into the Mechanism of the Cardiac Drug Omecamtiv Mecarbil─A Computational Study.

Omecamtiv mecarbil (OM) is a positive inotrope that is thought to bind directly to an allosteric site of the ß-cardiac myosin. The drug is under investigation for the treatment of systolic heart failure. The drug is classified as a cardiac myosin modulator and has been observed to affect multiple vital steps of the cross-bridge cycle including the recovery stroke and the chemical step. We explored the free-energy surface of the recovery stroke of the human cardiac ß-myosin in the presence of OM to determine its influence on this process. We also investigated the effects of OM on the recovery stroke in the presence of genetic cardiomyopathic mutations R712L, F764L, and P710R using metadynamics. We also utilized the method of transition path sampling to generate an unbiased ensemble of reactive trajectories for the ATP hydrolysis step in the presence of OM that were able to provide insight into the differences observed due to OM in the dynamics and mechanism of the decomposition of ATP to ADP and HPO42-, a central part of the power generation in cardiac muscle. We studied chemistry in the presence of the same three mutations to further elucidate the effect of OM, and its use in the treatment of cardiac disease.

Perspective: Path Sampling Methods Applied to Enzymatic Catalysis.

This Perspective reviews the use of Transition Path Sampling methods to study enzymatically catalyzed chemical reactions. First applied by our group to an enzymatic reaction over 15 years ago, the method has uncovered basic principles in enzymatic catalysis such as the protein promoting vibration, and it has also helped harmonize such ideas as electrostatic preorganization with dynamic views of enzyme function. It is now being used to help uncover principles of protein design necessary to artificial enzyme creation.

Transition Path Sampling Based Calculations of Free Energies for Enzymatic Reactions: The Case of Human Methionine Adenosyl Transferase and Plasmodium vivax Adenosine Deaminase.

Transition path sampling (TPS) is widely used for the calculations of reaction rates, transition state structures, and reaction coordinates of condensed phase systems. Here we discuss a scheme for the calculation of free energies using the ensemble of TPS reactive trajectories in combination with a window-based sampling technique for enzyme-catalyzed reactions. We calculate the free energy profiles of the reactions catalyzed by the human methionine S-adenosyltransferase (MAT2A) enzyme and the Plasmodium vivax adenosine deaminase (pvADA) enzyme to assess the accuracy of this method. MAT2A catalyzes the formation of S-adenosine-l-methionine following a SN2 mechanism, and using our method, we estimate the free energy barrier for this reaction to be 16 kcal mol-1, which is in excellent agreement with the experimentally measured activation energy of 17.27 kcal mol-1. Furthermore, for the pvADA enzyme-catalyzed reaction we estimate a free energy barrier of 21 kcal mol-1, and the calculated free energy profile is similar to that predicted from experimental observations. Calculating free energies by employing our simple method within TPS provides significant advantages over methods such as umbrella sampling because it is free from any applied external bias, is accurate compared to experimental measurements, and has a reasonable computational cost.

Structure and Dynamics of the Flexible Cardiac Troponin T Linker Domain in a Fully Reconstituted Thin Filament.

The structural analysis of large protein complexes has been greatly enhanced through the application of electron microscopy techniques. One such multiprotein complex, the cardiac thin filament (cTF), has cyclic interactions with thick filament proteins to drive contraction of the heart that has recently been the subject of such studies. As important as these studies are, they provide limited or no information on highly flexible regions that in isolation would be characterized as inherently disordered. One such region is the extended cardiac troponin T (cTnT) linker between the regions of cTnT which have been labeled TNT1 and TNT2. It comprises a hinge region (residues 158-166) and a highly flexible region (residues 167-203). Critically, this region modulates the troponin/tropomyosin complex's position across the actin filament. Thus, the cTnT linker structure and dynamics are central to the regulation of the function of cardiac muscles, but up to now, it was ill-understood. To establish the cTnT linker structure, we coupled an atomistic computational cTF model with time-resolved fluorescence resonance energy transfer measurements in both ±Ca2+ conditions utilizing fully reconstituted cTFs. We mapped the cTnT linker's positioning across the actin filament, and by coupling the experimental results to computation, we found mean structures and ranges of motion of this part of the complex. With this new insight, we can now address cTnT linker structural dynamics in both myofilament activation and disease.

Engineered Tryptophan Synthase Balances Equilibrium Effects and Fast Dynamic Effects.

Creating efficient and stable enzymes for catalysis in pharmaceutical and industrial laboratories is an important research goal. Arnold et al. used directed evolution to engineer a natural tryptophan synthase to create a mutant that is operable under laboratory conditions without the need for a natural allosteric effector. The use of directed evolution allows researchers to improve enzymes without understanding the structure-activity relationship. Here, we present a transition path sampling study of a key chemical transformation in the tryptophan synthase catalytic cycle. We observed that while directed evolution does mimic the natural allosteric effect from a stability perspective, fast protein dynamics associated with chemistry has been dramatically altered. This work provides further evidence of the role of protein dynamics in catalysis and clearly demonstrates the multifaceted complexity of mutations associated with protein engineering. This study also demonstrates a fascinating contrast between allosteric and stand-alone functions at the femtosecond time scale.

Free-Energy Surfaces of Two Cardiac Thin Filament Conformational Changes during Muscle Contraction.

The troponin core is an important regulatory complex in cardiac sarcomeres. Contraction is initiated by a calcium ion binding to cardiac troponin C (cTnC), initiating a conformational shift within the protein, altering its interactions with cardiac troponin I (cTnI). The change in cTnC-cTnI interactions prompts the C-terminal domain of cTnI to dissociate from actin, allowing tropomyosin to reveal myosin-binding sites on actin. Each of the concerted movements in the cardiac thin filament (CTF) is crucial for allowing the contraction of cardiomyocytes, yet little is known about the free energy associated with each transition, which is vital for understanding contraction on a molecular level. Using metadynamics, we calculated the free-energy surface of two transitions in the CTF: cTnC opening in the presence and absence of Ca2+ and cTnI dissociating from actin with both open and closed cTnC. These results not only provide the free-energy surface of the transitions but will also be shown to determine if the order of transitions in the contraction cycle is important. From our calculations, we found that the calcium ion helps stabilize the open conformation of cTnC and that the C-terminus of cTnI is stabilized by cTnC in the open conformation when dissociating from the actin surface.

Method for Identifying Common Features in Reactive Trajectories of a Transition Path Sampling Ensemble.

Simulation methods like transition path sampling (TPS) generate an abundance of information buried in the collection of reactive trajectories that they generate. However, only limited use has been made of this information, mainly for the identification of the reaction coordinate. The standard TPS tools have been designed for monitoring the progress of the system from reactants to products. However, the reaction coordinate does not contain all the information regarding the mechanism. In our earlier work, we have used TPS on enzymatic systems and have identified important motions in the reactant well that prepares the system for the reaction. Since these events take place in the reactant well, they are beyond the reach of standard TPS postprocessing methods. We present a simple scheme for identifying the common trends in enzymatic trajectories. This scheme was designed for a specific class of enzymatic reactions: it can be used for identifying motions that guide the system to reaction-ready conformations. We have applied it to two enzymatic systems that we have studied in the past, formate dehydrogenase and purine nucleoside phosphorylase, and we were able to identify interactions, far from the transition state, that are important for preparing the system for the reaction but that had been overlooked in earlier work.

RAVIR: A Dataset and Methodology for the Semantic Segmentation and Quantitative Analysis of Retinal Arteries and Veins in Infrared Reflectance Imaging.

The retinal vasculature provides important clues in the diagnosis and monitoring of systemic diseases including hypertension and diabetes. The microvascular system is of primary involvement in such conditions, and the retina is the only anatomical site where the microvasculature can be directly observed. The objective assessment of retinal vessels has long been considered a surrogate biomarker for systemic vascular diseases, and with recent advancements in retinal imaging and computer vision technologies, this topic has become the subject of renewed attention. In this paper, we present a novel dataset, dubbed RAVIR, for the semantic segmentation of Retinal Arteries and Veins in Infrared Reflectance (IR) imaging. It enables the creation of deep learning-based models that distinguish extracted vessel type without extensive post-processing. We propose a novel deep learning-based methodology, denoted as SegRAVIR, for the semantic segmentation of retinal arteries and veins and the quantitative measurement of the widths of segmented vessels. Our extensive experiments validate the effectiveness of SegRAVIR and demonstrate its superior performance in comparison to state-of-the-art models. Additionally, we propose a knowledge distillation framework for the domain adaptation of RAVIR pretrained networks on color images. We demonstrate that our pretraining procedure yields new state-of-the-art benchmarks on the DRIVE, STARE, and CHASE_DB1 datasets. Dataset link: https://ravirdataset.github.io/data.

Ocular Toxoplasmosis: No Stranger to the Masquerade Ball.

Purpose: This article illustrates multiple atypical manifestations of ocular toxoplasmosis masquerading as acute retinal necrosis and vitreoretinal lymphoma. Methods: Two case presentations are discussed, and the body of pertinent literature is reviewed and discussed. Results: In these cases, an extensive workup and attention to history lead to the correct diagnosis and management. Conclusions: Aggressive cases of ocular toxoplasmosis may present in a variety of phenotypes that may mimic other vision- and potentially life-threatening conditions, particularly in a milieu of inadequate endogenous and exogenous antimicrobial defenses.

Epiretinal proliferation after rhegmatogenous retinal detachment.

PURPOSE: To determine the characteristics and appearance rate of epiretinal proliferation (ERP) on SD-OCT after surgery for rhegmatogenous retinal detachment (RRD) repair. METHODS: One hundred eight eyes of 108 patients who underwent one or more surgeries for RRD were enrolled. The eyes with other maculopathies that were directly related to RRD were excluded. Image acquisition was performed with SD-OCT (Heidelberg Engineering, Germany). Clinical charts were reviewed to assess clinical and surgical findings. Statistical analyses were performed using XLSTAT (Assinsoft, Paris, France). RESULTS: ERP was found in 9.3% eyes (n = 10). The mean initial visual acuity (logMAR) was 1.34 ± 0.82 in the ERP group compared to 0.49 ± 0.70 in the non-ERP group. PVR was present in 70.0% and chronic macular edema was found in 80.0% of eyes which developed ERP. The mean number of vitreoretinal surgeries in eyes with ERP was 3.3 ± 1.19 and only 1.44 ± 1.02 in eyes without. Silicone oil was used in 60.0% of eyes which developed ERP compared to 13.9% in the non-ERP group. CONCLUSION: ERP is a late-onset postoperative finding in eyes with RRD and can occur in absence of macular holes. Overall, ERP is more frequent in eyes with complicated courses of RRD including multiple operations, PVR, usage of silicone oil, and chronic macular edema.