Deprotonation of formic, acetic acids and bicarbonate ion in slit silica nanopores at infinite dilution and in the presence of electrolytes.

Dielectric effects and the coupled electrostatics between the nanoconfined and the internal/external aqueous media contribute to the observed deviations of chemistry within the nanoconfined environment when compared with unconfined systems. A systematic understanding has remained elusive, especially with respect to background salt concentration and boundary condition effects like the nanopore surface chemistry and the reference state used to calculate free energies. We utilize molecular dynamics simulations along with thermodynamic integration to determine the free energy difference associated with acid-base chemistry in 2 nm and 4 nm slit pores open to a bulk-like reservoir. pKa increases are predicted when confining acetic acid, formic acid, and bicarbonate in the slits at infinite dilution conditions. We find that confinement weakens the acids, and the modulation of outer pore surface dipole magnitudes can tune the pKa shift values, suggesting that purely "intrinsic" electrostatic effect on confinement may not exist. At sufficiently high salt concentrations, the dielectric/electrostatic effects on pKa values diminish due to charge screening effects. These discoveries enable future modifications of nanopore chemistries to achieve desirable properties for industrial applications.

Formic and acetic acid pKa values increase under nanoconfinement.

Organic acids are prevalent in the environment and their acidity and the corresponding dissociation constants can change under varying environmental conditions. The impact of nanoconfinement (when acids are confined within nanometer-scale domains) on physicochemical properties of chemical species is poorly understood and is an emerging field of study. By combining infrared and Raman spectroscopies with molecular dynamics (MD) simulations, we quantified the effect of nanoconfinement in silica nanopores on one of the fundamental chemical reactions-the dissociation of organic acids. The pKa of formic and acetic acids confined within cylindrical silica nanopores with 4 nm diameters were measured. MD models were constructed to calculate the shifts in the pKa values of acetic acid nanoconfined within 1, 2, 3, and 4 nm silica slit pores. Both experiments and MD models indicate a decrease in the apparent acid dissociation constants (i.e., increase in the pKa values) when organic acids are nanoconfined. Therefore, nanoconfinement stabilizes the protonated species. We attribute this observation to (1) a decrease in the average dielectric response of nanoconfined aqueous solutions where charge screening may be decreased; or (2) an increase in proton concentration inside nanopores, which would shift the equilibrium towards the protonated form. Overall, the results of this study provide the first quantification of the pKa values for nanoconfined formic and acetic acids and pave the way for a unifying theory predicting the impact of nanoconfinement on acid-base chemistry.

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.

Computational and biophysical determination of pathogenicity of variants of unknown significance in cardiac thin filament.

Point mutations within sarcomeric proteins have been associated with altered function and cardiomyopathy development. Difficulties remain, however, in establishing the pathogenic potential of individual mutations, often limiting the use of genotype in management of affected families. To directly address this challenge, we utilized our all-atom computational model of the human full cardiac thin filament (CTF) to predict how sequence substitutions in CTF proteins might affect structure and dynamics on an atomistic level. Utilizing molecular dynamics calculations, we simulated 21 well-defined genetic pathogenic cardiac troponin T and tropomyosin variants to establish a baseline of pathogenic changes induced in computational observables. Computational results were verified via differential scanning calorimetry on a subset of variants to develop an experimental correlation. Calculations were performed on 9 independent variants of unknown significance (VUS), and results were compared with pathogenic variants to identify high-resolution pathogenic signatures. Results for VUS were compared with the baseline set to determine induced structural and dynamic changes, and potential variant reclassifications were proposed. This unbiased, high-resolution computational methodology can provide unique structural and dynamic information that can be incorporated into existing analyses to facilitate classification both for de novo variants and those where established approaches have provided conflicting information.

Investigation of the Recovery Stroke and ATP Hydrolysis and Changes Caused Due to the Cardiomyopathic Point Mutations in Human Cardiac ß Myosin.

Human cardiac ß myosin undergoes the cross-bridge cycle as part of the force-generating mechanism of cardiac muscle. The recovery stroke is considered one of the key steps of the kinetic cycle as it is the conformational rearrangement required to position the active site residues for hydrolysis of ATP and interaction with actin. We explored the free-energy surface of the transition and investigated the effect of the genetic cardiomyopathy causing mutations R453C, I457T, and I467T on this step using metadynamics. This work extends previous studies on Dictyostelium myosin II with engineered mutations. Here, like previously, we generated an unbiased thermodynamic ensemble of reactive trajectories for the chemical step using transition path sampling. Our methodologies were able to predict the changes to the dynamics of the recovery stroke as well as predict the pathway of breakdown of ATP to ADP and HPO42- with the stabilization of the metaphosphate intermediate. We also observed clear differences between the Dictyostelium myosin II and human cardiac ß myosin for ATP hydrolysis as well as predict the effect of the mutation I467T on the chemical step.

A Proposed Mechanism for the Initial Myosin Binding Event on the Cardiac Thin Filament: A Metadynamics Study.

The movement of tropomyosin over filamentous actin regulates the cross-bridge cycle of the thick with thin filament of cardiac muscle by blocking and revealing myosin binding sites. Tropomyosin exists in three, distinct equilibrium states with one state blocking myosin-actin interactions (blocked position) and the remaining two allowing for weak (closed position) and strong myosin binding (open position). However, experimental information illuminating how myosin binds to the thin filament and influences tropomyosin's transition across the actin surface is lacking. Using metadynamics, we mimic the effect of a single myosin head binding by determining the work required to pull small segments of tropomyosin toward the open position in several distinct regions of the thin filament. We find differences in required work due to the influence of cardiac troponin T lead to preferential binding sites and determine the mechanism of further myosin head recruitment.

Mechanochemical Function of Myosin II: Investigation into the Recovery Stroke and ATP Hydrolysis.

Myosin regulates muscle function through a complex cycle of conformational rearrangements coupled with the hydrolysis of adenosine triphosphate (ATP). The recovery stroke reorganizes the myosin active site to hydrolyze ATP and cross bridge with the thin filament to produce muscle contraction. Engineered mutations K84M and R704E in Dictyostelium myosin have been designed to specifically inhibit the recovery stroke and have been shown to indirectly affect the ATPase activity of myosin. We investigated these mutagenic perturbations to the recovery stroke and generated thermodynamically correct and unbiased trajectories for native ATP hydrolysis with computationally enhanced sampling methods. Our methodology was able to resolve experimentally observed changes to kinetic and equilibrium dynamics for the recovery stroke with the correct prediction in the severity of these changes. For ATP hydrolysis, the sequential nature along with the stabilization of a metaphosphate intermediate was observed in agreement with previous studies. However, we observed glutamate 459 being utilized as a proton abstractor to prime the attacking water instead of a lytic water, a phenomenon not well categorized in myosin but has in other ATPases. Both rare event methodologies can be extended to human myosin to investigate isoformic differences from Dictyostelium and scan cardiomyopathic mutations to see differential perturbations to kinetics of other conformational changes in myosin such as the power stroke.

FRET-based analysis of the cardiac troponin T linker region reveals the structural basis of the hypertrophic cardiomyopathy-causing Δ160E mutation.

Mutations in the cardiac thin filament (TF) have highly variable effects on the regulatory function of the cardiac sarcomere. Understanding the molecular-level dysfunction elicited by TF mutations is crucial to elucidate cardiac disease mechanisms. The hypertrophic cardiomyopathy-causing cardiac troponin T (cTnT) mutation Δ160Glu (Δ160E) is located in a putative "hinge" adjacent to an unstructured linker connecting domains TNT1 and TNT2. Currently, no high-resolution structure exists for this region, limiting significantly our ability to understand its role in myofilament activation and the molecular mechanism of mutation-induced dysfunction. Previous regulated in vitro motility data have indicated mutation-induced impairment of weak actomyosin interactions. We hypothesized that cTnT-Δ160E repositions the flexible linker, altering weak actomyosin electrostatic binding and acting as a biophysical trigger for impaired contractility and the observed remodeling. Using time-resolved FRET and an all-atom TF model, here we first defined the WT structure of the cTnT-linker region and then identified Δ160E mutation-induced positional changes. Our results suggest that the WT linker runs alongside the C terminus of tropomyosin. The Δ160E-induced structural changes moved the linker closer to the tropomyosin C terminus, an effect that was more pronounced in the presence of myosin subfragment (S1) heads, supporting previous findings. Our in silico model fully supported this result, indicating a mutation-induced decrease in linker flexibility. Our findings provide a framework for understanding basic pathogenic mechanisms that drive severe clinical hypertrophic cardiomyopathy phenotypes and for identifying structural targets for intervention that can be tested in silico and in vitro.

Proof of Principle that Molecular Modeling Followed by a Biophysical Experiment Can Develop Small Molecules that Restore Function to the Cardiac Thin Filament in the Presence of Cardiomyopathic Mutations.

This article reports a coupled computational experimental approach to design small molecules aimed at targeting genetic cardiomyopathies. We begin with a fully atomistic model of the cardiac thin filament. To this we dock molecules using accepted computational drug binding methodologies. The candidates are screened for their ability to repair alterations in biophysical properties caused by mutation. Hypertrophic and dilated cardiomyopathies caused by mutation are initially biophysical in nature, and the approach we take is to correct the biophysical insult prior to irreversible cardiac damage. Candidate molecules are then tested experimentally for both binding and biophysical properties. This is a proof of concept study-eventually candidate molecules will be tested in transgenic animal models of genetic (sarcomeric) cardiomyopathies.