Evaluating sheep hemoglobins with MD simulations as an animal model for sickle cell disease

Sickle cell disease (SCD) affects millions worldwide, yet there are few therapeutic options. To develop effective treatments, preclinical models that recapitulate human physiology and SCD pathophysiology are needed. SCD arises from a single Glu-to-Val substitution at position 6 in the β subunit of hemoglobin (Hb), promoting Hb polymerization and subsequent disease. Sheep share important physiological and developmental characteristics with humans, including the same developmental pattern of fetal to adult Hb switching. Herein, we investigated whether introducing the SCD mutation into the sheep β-globin locus would recapitulate SCD's complex pathophysiology by generating high quality SWISS-MODEL sheep Hb structures and performing MD simulations of normal/sickle human (huHbA/huHbS) and sheep (shHbB/shHbS) Hb, establishing how accurately shHbS mimics huHbS behavior. shHbS, like huHbS, remained stable with low RMSD, while huHbA and shHbB had higher and fluctuating RMSD. shHbB and shHbS also behaved identically to huHbA and huHbS with respect to β2-Glu6 and β1-Asp73 (β1-Asn72 in sheep) solvent interactions. These data demonstrate that introducing the single SCD-causing Glu-to-Val substitution into sheep β-globin causes alterations consistent with the Hb polymerization that drives RBC sickling, supporting the development of a SCD sheep model to pave the way for alternative cures for this debilitating, globally impactful disease.

How soluble misfolded proteins bypass chaperones at the molecular level

Subpopulations of soluble, misfolded proteins can bypass chaperones within cells. The extent of this phenomenon and how it happens at the molecular level are unknown. Through a meta-analysis of the experimental literature we find that in all quantitative protein refolding studies there is always a subpopulation of soluble but misfolded protein that does not fold in the presence of one or more chaperones, and can take days or longer to do so. Thus, some misfolded subpopulations commonly bypass chaperones. Using multi-scale simulation models we observe that the misfolded structures that bypass various chaperones can do so because their structures are highly native like, leading to a situation where chaperones do not distinguish between the folded and near-native-misfolded states. More broadly, these results provide a mechanism by which long-time scale changes in protein structure and function can persist in cells because some misfolded states can bypass components of the proteostasis machinery.

Computational and Experimental Characterization of rDNA and rRNA G-Quadruplexes

DNA G-quadruplexes in human telomeres and gene promoters are being extensively studied for their role in controlling the growth of cancer cells. G-quadruplexes have been unambiguously shown to exist both in vitro and in vivo, including in the guanine (G)-rich DNA genes encoding pre-ribosomal RNA (pre-rRNA), which is transcribed in the cell's nucleolus. Recent studies strongly suggest that these DNA sequences ("rDNA"), and the transcribed rRNA, are a potential anticancer target through the inhibition of RNA polymerase I (Pol I) in ribosome biogenesis, but the structures of ribosomal G-quadruplexes at atomic resolution are unknown and very little biophysical characterization has been performed on them to date. In the present study, circular dichroism (CD) spectroscopy is used to show that two putative rDNA G-quadruplex sequences, NUC 19P and NUC 23P and their counterpart rRNAs, predominantly adopt parallel topologies, reminiscent of the analogous telomeric quadruplex structures. Based on this information, we modeled parallel topology atomistic structures of the putative ribosomal G-quadruplexes. We then validated and refined the modeled ribosomal G-quadruplex structures using all-atom molecular dynamics (MD) simulations with the CHARMM36 force field in the presence and absence of stabilizing K⁺. Motivated by preliminary MD simulations of the telomeric parallel G-quadruplex (TEL 24P) in which the K⁺ ion is expelled, we used updated CHARMM36 force field K⁺ parameters that were optimized, targeting the data from quantum mechanical calculations and the polarizable Drude model force field. In subsequent MD simulations with optimized CHARMM36 parameters, the K⁺ ions are predominantly in the G-quadruplex channel and the rDNA G-quadruplexes have more well-defined, predominantly parallel-topology structures as compared to rRNA. In addition, NUC 19P is more structured than NUC 23P, which contains extended loops. Results from this study set the structural foundation for understanding G-quadruplex functions and the design of novel chemotherapeutics against these nucleolar targets and can be readily extended to other DNA and RNA G-quadruplexes.

Martini Coarse-Grained Force Field: Extension to RNA

RNA has an important role not only as the messenger of genetic information but also as a regulator of gene expression. Given its central role in cell biology, there is significant interest in studying the structural and dynamic behavior of RNA in relation to other biomolecules. Coarse-grain molecular dynamics simulations are a key tool to that end. Here, we have extended the coarse-grain Martini force field to include RNA after our recent extension to DNA. In the same way DNA was modeled, the tertiary structure of RNA is constrained using an elastic network. This model, therefore, is not designed for applications involving RNA folding but rather offers a stable RNA structure for studying RNA interactions with other (bio)molecules. The RNA model is compatible with all other Martini models and opens the way to large-scale explicit-solvent molecular dynamics simulations of complex systems involving RNA.

Nucleic acid polymeric properties and electrostatics: Directly comparing theory and simulation with experiment

Nucleic acids are biopolymers that carry genetic information and are also involved in various gene regulation functions such as gene silencing and protein translation. Because of their negatively charged backbones, nucleic acids are polyelectrolytes. To adequately understand nucleic acid folding and function, we need to properly describe its i) polymer/polyelectrolyte properties and ii) associating ion atmosphere. While various theories and simulation models have been developed to describe nucleic acids and the ions around them, many of these theories/simulations have not been well evaluated due to complexities in comparison with experiment. In this review, I discuss some recent experiments that have been strategically designed for straightforward comparison with theories and simulation models. Such data serve as excellent benchmarks to identify limitations in prevailing theories and simulation parameters.

MD Simulations of tRNA and Aminoacyl-tRNA Synthetases: Dynamics, Folding, Binding, and Allostery

While tRNA and aminoacyl-tRNA synthetases are classes of biomolecules that have been extensively studied for decades, the finer details of how they carry out their fundamental biological functions in protein synthesis remain a challenge. Recent molecular dynamics (MD) simulations are verifying experimental observations and providing new insight that cannot be addressed from experiments alone. Throughout the review, we briefly discuss important historical events to provide a context for how far the field has progressed over the past few decades. We then review the background of tRNA molecules, aminoacyl-tRNA synthetases, and current state of the art MD simulation techniques for those who may be unfamiliar with any of those fields. Recent MD simulations of tRNA dynamics and folding and of aminoacyl-tRNA synthetase dynamics and mechanistic characterizations are discussed. We highlight the recent successes and discuss how important questions can be addressed using current MD simulations techniques. We also outline several natural next steps for computational studies of AARS:tRNA complexes.