Examinations of tRNA Range of Motion Using Simulations of Cryo-EM Microscopy and X-Ray Data

In Silico Investigation of Parkin-Activating Mutations Using Simulations and Network Modeling

Complete loss-of-function mutations in the PRKN gene are a major cause of early-onset Parkinson's disease (PD). PRKN encodes the Parkin protein, an E3 ubiquitin ligase that works in conjunction with the ubiquitin kinase PINK1 in a distinct quality control pathway to tag damaged mitochondria for autophagic clearance, i.e., mitophagy. According to previous structural investigations, Parkin protein is typically kept in an inactive conformation via several intramolecular, auto-inhibitory interactions. Here, we performed molecular dynamics simulations (MDS) to provide insights into conformational changes occurring during the de-repression of Parkin and the gain of catalytic activity. We analyzed four different Parkin-activating mutations that are predicted to disrupt certain aspects of its auto-inhibition. All four variants showed greater conformational motions compared to wild-type protein, as well as differences in distances between domain interfaces and solvent-accessible surface area, which are thought to play critical roles as Parkin gains catalytic activity. Our findings reveal that the studied variants exert a notable influence on Parkin activation as they alter the opening of its closed inactive structure, a finding that is supported by recent structure- and cell-based studies. These findings not only helped further characterize the hyperactive variants but overall improved our understanding of Parkin's catalytic activity and nominated targets within Parkin's structure for potential therapeutic designs.

Protein structure-based in-silico approaches to drug discovery: Guide to COVID-19 therapeutics

With more than 5 million fatalities and close to 300 million reported cases, COVID-19 is the first documented pandemic due to a coronavirus that continues to be a major health challenge. Despite being rapid, uncontrollable, and highly infectious in its spread, it also created incentives for technology development and redefined public health needs and research agendas to fast-track innovations to be translated. Breakthroughs in computational biology peaked during the pandemic with renewed attention to making all cutting-edge technology deliver agents to combat the disease. The demand to develop effective treatments yielded surprising collaborations from previously segregated fields of science and technology. The long-standing pharmaceutical industry's aversion to repurposing existing drugs due to a lack of exponential financial gain was overrun by the health crisis and pressures created by front-line researchers and providers. Effective vaccine development even at an unprecedented pace took more than a year to develop and commence trials. Now the emergence of variants and waning protections during the booster shots is resulting in breakthrough infections that continue to strain health care systems. As of now, every protein of SARS-CoV-2 has been structurally characterized and related host pathways have been extensively mapped out. The research community has addressed the druggability of a multitude of possible targets. This has been made possible due to existing technology for virtual computer-assisted drug development as well as new tools and technologies such as artificial intelligence to deliver new leads. Here in this article, we are discussing advances in the drug discovery field related to target-based drug discovery and exploring the implications of known target-specific agents on COVID-19 therapeutic management. The current scenario calls for more personalized medicine efforts and stratifying patient populations early on for their need for different combinations of prognosis-specific therapeutics. We intend to highlight target hotspots and their potential agents, with the ultimate goal of using rational design of new therapeutics to not only end this pandemic but also uncover a generalizable platform for use in future pandemics.

Substitution of PINK1 Gly411 modulates substrate receptivity and turnover

The ubiquitin (Ub) kinase-ligase pair PINK1-PRKN mediates the degradation of damaged mitochondria by macroautophagy/autophagy (mitophagy). PINK1 surveils mitochondria and upon stress accumulates on the mitochondrial surface where it phosphorylates serine 65 of Ub to activate PRKN and to drive mitochondrial turnover. While loss of either PINK1 or PRKN is genetically linked to Parkinson disease (PD) and activating the pathway seems to have great therapeutic potential, there is no formal proof that stimulation of mitophagy is always beneficial. Here we used biochemical and cell biological methods to study single nucleotide variants in the activation loop of PINK1 to modulate the enzymatic function of this kinase. Structural modeling and in vitro kinase assays were used to investigate the molecular mechanism of the PINK1 variants. In contrast to the PD-linked PINK1^G411S mutation that diminishes Ub kinase activity, we found that the PINK1^G411A variant significantly boosted Ub phosphorylation beyond levels of PINK1 wild type. This resulted in augmented PRKN activation, mitophagy rates and increased viability after mitochondrial stress in midbrain-derived, gene-edited neurons. Mechanistically, the G411A variant stabilizes the kinase fold of PINK1 and transforms Ub to adopt the preferred, C-terminally retracted conformation for improved substrate turnover. In summary, we identify a critical role of residue 411 for substrate receptivity that may now be exploited for drug discovery to increase the enzymatic function of PINK1. The genetic substitution of Gly411 to Ala increases mitophagy and may be useful to confirm neuroprotection in vivo and might serve as a critical positive control during therapeutic development.Abbreviations: ATP: adenosine triphosphate; CCCP: carbonyl cyanide m-chlorophenyl hydrazone; Ub-CR: ubiquitin with C-terminally retracted tail; CTD: C-terminal domain (of PINK1); ELISA: enzyme-linked immunosorbent assay; HCI: high-content imaging; IB: immunoblot; IF: immunofluorescence; NPC: neuronal precursor cells; MDS: molecular dynamics simulation; PD: Parkinson disease; p-S65-Ub: ubiquitin phosphorylated at Ser65; RMSF: root mean scare fluctuation; TOMM: translocase of outer mitochondrial membrane; TVLN: ubiquitin with T66V and L67N mutation, mimics Ub-CR; Ub: ubiquitin; WT: wild-type.

Molecular Modeling and Phenotypic Description of a Patient with a Novel Exonic Deletion of GALNS with Resultant Morquio Syndrome with Two Successful Pregnancies

In this report, we describe phenotypic features of a patient with mucopolysaccharidosis type IVA (Morquio syndrome) harboring a novel exon 1 deletion in GALNS with enzymatic confirmation consistent with Morquio syndrome. To our knowledge, this is the first reported case of this variant. Additionally, we protein modelled wild-type GALNS and the pathogenic variant with an exon 1 deletion for comparative analysis using statistical mechanics methods described herein. We demonstrate that, even when the protein is translated, the mutation would affect protein stability and function via homodimer interaction modifications. Lastly, given the patient's 2 successful pregnancies, data about the management of pregnancies in mucopolysaccharidoses are reviewed, and we discuss the management of pregnancy in patients with Morquio syndrome.

Attacking COVID-19 Progression Using Multi-Drug Therapy for Synergetic Target Engagement

COVID-19 is a devastating respiratory and inflammatory illness caused by a new coronavirus that is rapidly spreading throughout the human population. Over the past 12 months, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus responsible for COVID-19, has already infected over 160 million (>20% located in United States) and killed more than 3.3 million people around the world (>20% deaths in USA). As we face one of the most challenging times in our recent history, there is an urgent need to identify drug candidates that can attack SARS-CoV-2 on multiple fronts. We have therefore initiated a computational dynamics drug pipeline using molecular modeling, structure simulation, docking and machine learning models to predict the inhibitory activity of several million compounds against two essential SARS-CoV-2 viral proteins and their host protein interactors-S/Ace2, Tmprss2, Cathepsins L and K, and Mpro-to prevent binding, membrane fusion and replication of the virus, respectively. All together, we generated an ensemble of structural conformations that increase high-quality docking outcomes to screen over >6 million compounds including all FDA-approved drugs, drugs under clinical trial (>3000) and an additional >30 million selected chemotypes from fragment libraries. Our results yielded an initial set of 350 high-value compounds from both new and FDA-approved compounds that can now be tested experimentally in appropriate biological model systems. We anticipate that our results will initiate screening campaigns and accelerate the discovery of COVID-19 treatments.

Clinical description & molecular modeling of novel MAX pathogenic variant causing pheochromocytoma in family, supports paternal parent-of-origin effect

The titular member of the MAX network of proteins, MYC-associated factor X (MAX), serves an important regulatory function in transcription of E-box genes associated with cell proliferation, differentiation, and apoptosis. Wild type MAX dimerizes with both MYC and MAD, both of which are members of the MAX network, and can promote or repress cell functions as needed. However, pathogenic variants in MAX are known to upset this balance, leading to uncontrolled oncogenic activity and disease phenotypes such as paragangliomas and pheochromocytomas. We report a 58-year-old male and his 32-year-old daughter, both of which have a history of pheochromocytoma and the unique nonsense MAX variant c.271C>T (p.Q91X). These individuals were diagnosed with pheochromocytomas in their early twenties that were later removed through corrective surgery. The father now presents with recurrent symptoms of hypertension, hyperhidrosis, and headaches, which accompany new pheochromocytomas of his remaining adrenal gland. Pathogenicity of this MAX variant is proven through molecular modeling. The case of this father-daughter pair supports both heritability of pheochromocytoma and the paternal parent-of-origin effect for MAX pathogenic variants.

A Virtual Screening Platform Identifies Chloroethylagelastatin A as a Potential Ribosomal Inhibitor

Chloroethylagelastatin A (CEAA) is an analogue of agelastatin A (AA), a natural alkaloid derived from a marine sponge. It is under development for therapeutic use against brain tumors as it has excellent central nervous system (CNS) penetration and pre-clinical therapeutic activity against brain tumors. Recently, AA was shown to inhibit protein synthesis by binding to the ribosomal A-site. In this study, we developed a novel virtual screening platform to perform a comprehensive screening of various AA analogues showing that AA analogues with proven therapeutic activity including CEAA have significant ribosomal binding capacity whereas therapeutically inactive analogues show poor ribosomal binding and revealing structural fingerprint features essential for drug-ribosome interactions. In particular, CEAA was found to have greater ribosomal binding capacity than AA. Biological tests showed that CEAA binds the ribosome and contributes to protein synthesis inhibition. Our findings suggest that CEAA may possess ribosomal inhibitor activity and that our virtual screening platform may be a useful tool in discovery and development of novel ribosomal inhibitors.

Genomics combined with a protein informatics platform to assess a novel pathogenic variant c.1024 A>G (p.K342E) in OPA1 in a patient with autosomal dominant optic atrophy

BACKGROUND

Autosomal Dominant Optic Atrophy (ADOA) is caused by mutations in the Optic Atrophy 1 Gene which disrupts the OPA1 protein. This disruption affects the normal function of the protein; impairs fusion of the mitochondrial inner membrane; and prevents normal OPA1 protein degradation. These events cause damage in retinal ganglion cells that could affect the patients with symptoms ranging from none to legally blind.

MATERIALS AND METHODS

Our study identifies a missense variant mutation, c.1024 A > G (p.K342E), in OPA1 gene causing ADOA. Diagnosed clinically in three family members and the presence of this mutation was confirmed in two members by genetic testing. Pathogenic variants in OPA1 impact the secondary protein structure and function by causing non-conservative amino acid substitutions. We also modeled this mutation and compared it to the wild type using statistical mechanics.

RESULTS AND CONCLUSIONS

The proband's pathogenic variant, c.1024 A > G (p.K342E), is located in the GTPase domain of OPA1 and causes changes in the protein structure by affecting the oligomerization pattern thus resulting in ADOA. Identifying the pathogenic potential of the missense mutations in the OPA1 gene using neoteric protein modeling techniques would help in the early detection of ADOA in patients who have family history of blindness. This action would help in providing early follow up, possible treatment in the future, and genetic counseling. Abbreviations: ADOA: Autosomal Dominant Optic Atrophy; CYCS: Caspase Activator Cytochrome C; OPA1: Optic Atrophy Gene 1; RGC: Retinal Ganglion Cells; VUS: Variant of Uncertain Significance.

Structural Models for the Dynamic Effects of Loss-of-Function Variants in the Human SIM1 Protein Transcriptional Activation Domain

Single-minded homologue 1 (SIM1) is a transcription factor with numerous different physiological and developmental functions. SIM1 is a member of the class I basic helix-loop-helix-PER-ARNT-SIM (bHLH-PAS) transcription factor family, that includes several other conserved proteins, including the hypoxia-inducible factors, aryl hydrocarbon receptor, neuronal PAS proteins, and the CLOCK circadian regulator. Recent studies of HIF-a-ARNT and CLOCK-BMAL1 protein complexes have revealed the organization of their bHLH, PASA, and PASB domains and provided insight into how these heterodimeric protein complexes form; however, experimental structures for SIM1 have been lacking. Here, we describe the first full-length atomic structural model for human SIM1 with its binding partner ARNT in a heterodimeric complex and analyze several pathogenic variants utilizing state-of-the-art simulations and algorithms. Using local and global positional deviation metrics, deductions to the structural basis for the individual mutants are addressed in terms of the deleterious structural reorganizations that could alter protein function. We propose new experiments to probe these hypotheses and examine an interesting SIM1 dynamic behavior. The conformational dynamics demonstrates conformational changes on local and global regions that represent a mechanism for dysfunction in variants presented. In addition, we used our ab initio hybrid model for further prediction of variant hotspots that can be engineered to test for counter variant (restoration of wild-type function) or basic research probe.

Examination of Molecular Effects of MYLK Deletion in a Patient with Extensive Aortic, Carotid, and Abdominal Dissections That Underlie the Genetic Dysfunction

We describe the phenotype of a patient with extensive aortic, carotid, and abdominal dissections. The proband was found to have a heterozygous deletion of exons 21-34 in MYLK, which is a rare finding, as deletions in this gene have been infrequently reported. We describe this finding following detection in a proband with an extensive history of aortic, carotid, and abdominal dissections. Neoteric molecular modeling techniques to help determine the impact of this deletion on protein function indicated loss of function due to lack of any kinase domain. We also provide the electrostatics calculations from the wild type and mutant variant. Through a combined multiomic approach of clinical, functional, and protein informatics, we arrive at a data fusion for determination of pathogenicity embedded within the genetic code for this particular genetic variant, which, as a platform, continues to broaden its scope across the field of variants of uncertain significance classification.

Structural And Computational Perspectives of Selectively Targeting Mutant Proteins

Diseases are often caused by mutant proteins. Many drugs have limited effectiveness and/or toxic side effects because of a failure to selectively target the disease-causing mutant variant, rather than the functional wild type protein. Otherwise, the drugs may even target different proteins with similar structural features. Designing drugs that successfully target mutant proteins selectively represents a major challenge. Decades of cancer research have led to an abundance of potential therapeutic targets, often touted to be "master regulators". For many of these proteins, there are no FDA-approved drugs available; for others, off-target effects result in dose-limiting toxicity. Cancer-related proteins are an excellent medium to carry the story of mutant-specific targeting, as the disease is both initiated and sustained by mutant proteins; furthermore, current chemotherapies generally fail at adequate selective distinction. This review discusses some of the challenges associated with selective targeting from a structural biology perspective, as well as some of the developments in algorithm approach and computational workflow that can be applied to address those issues. One of the most widely researched proteins in cancer biology is p53, a tumor suppressor. Here, p53 is discussed as a specific example of a challenging target, with contemporary drugs and methodologies used as examples of burgeoning successes. The oncogene KRAS, which has been described as "undruggable", is another extensively investigated protein in cancer biology. This review also examines KRAS to exemplify progress made towards selective targeting of diseasecausing mutant proteins. Finally, possible future directions relevant to the topic are discussed.

Characterization of a Pathogenic Variant in the ABCD1 Gene Through Protein Molecular Modeling

Background

The ATP-binding cassette, subfamily D, member 1 (ABCD1) protein is a peroxisomal half-transporter that allows for very long chain fatty acid (VLCFA) degradation. Pathogenic variants of ABCD1 cause VLCFAs to build up in various tissues and bodily fluids, resulting in a disorder called X-linked adrenoleukodystrophy (X-ALD). This disorder is most commonly marked by adrenocortical insufficiency and high VLCFA concentration, and has varying levels of neurological involvement depending on phenotype. For example, the Addison-only form of X-ALD has no neurological impact, while the cerebral form of X-ALD often causes severe sensory loss, motor function impairment, cognitive decline, and death.

Methods

A newly characterized and suspected pathogenic variant in ABCD1 cause VLCFAs to build up in various tissues and bodily fluids, resulting in a disorder called X-linked adrenoleukodystrophy (X-ALD). This disorder is most commonly marked by adrenocortical insufficiency and high VLCFA concentration, and has varying levels of neurological involvement depending on phenotype. For example, the Addison-only form of X-ALD has no neurological impact, while the cerebral form of X-ALD often causes severe sensory loss, motor function impairment, cognitive decline, and death.

Results

A case of adult onset adrenomyeloneuropathy (AMN) and a novel ABCD1 cause VLCFAs to build up in various tissues and bodily fluids, resulting in a disorder called X-linked adrenoleukodystrophy (X-ALD). This disorder is most commonly marked by adrenocortical insufficiency and high VLCFA concentration, and has varying levels of neurological involvement depending on phenotype. For example, the Addison-only form of X-ALD has no neurological impact, while the cerebral form of X-ALD often causes severe sensory loss, motor function impairment, cognitive decline, and death.

Conclusions

Data fusion from multiple sources was combined in a comprehensive approach yielding an enriched assessment of the patient's disease and prognosis. Molecular modeling was performed on the variant to better characterize its clinical significance and confirm pathogenicity.

Molecular Dynamics Simulations Suggest a Non-Doublet Decoding Model of -1 Frameshifting by tRNASer3

In-frame decoding in the ribosome occurs through canonical or wobble Watson-Crick pairing of three mRNA codon bases (a triplet) with a triplet of anticodon bases in tRNA. Departures from the triplet-triplet interaction can result in frameshifting, meaning downstream mRNA codons are then read in a different register. There are many mechanisms to induce frameshifting, and most are insufficiently understood. One previously proposed mechanism is doublet decoding, in which only codon bases 1 and 2 are read by anticodon bases 34 and 35, which would lead to -1 frameshifting. In E. coli, tRNA^Ser3_GCU can induce -1 frameshifting at alanine (GCA) codons. The logic of the doublet decoding model is that the Ala codon's GC could pair with the tRNA^Ser3's GC, leaving the third anticodon residue U36 making no interactions with mRNA. Under that model, a U36C mutation would still induce -1 frameshifting, but experiments refute this. We perform all-atom simulations of wild-type tRNA^Ser3, as well as a U36C mutant. Our simulations revealed a hydrogen bond between U36 of the anticodon and G1 of the codon. The U36C mutant cannot make this interaction, as it lacks the hydrogen-bond-donating H3. The simulation thus suggests a novel, non-doublet decoding mechanism for -1 frameshifting by tRNA^Ser3 at Ala codons.

Molecular Inhibitor of QSOX1 Suppresses Tumor Growth In Vivo

Quiescin sulfhydryl oxidase 1 (QSOX1) is an enzyme overexpressed by many different tumor types. QSOX1 catalyzes the formation of disulfide bonds in proteins. Because short hairpin knockdowns (KD) of QSOX1 have been shown to suppress tumor growth and invasion in vitro and in vivo, we hypothesized that chemical compounds inhibiting QSOX1 enzymatic activity would also suppress tumor growth, invasion, and metastasis. High throughput screening using a QSOX1-based enzymatic assay revealed multiple potential QSOX1 inhibitors. One of the inhibitors, known as "SBI-183," suppresses tumor cell growth in a Matrigel-based spheroid assay and inhibits invasion in a modified Boyden chamber, but does not affect viability of nonmalignant cells. Oral administration of SBI-183 inhibits tumor growth in 2 independent human xenograft mouse models of renal cell carcinoma. We conclude that SBI-183 warrants further exploration as a useful tool for understanding QSOX1 biology and as a potential novel anticancer agent in tumors that overexpress QSOX1.

Genomic Observations of a Rare/Pathogenic SMAD3 Variant in Loeys⁻Dietz Syndrome 3 Confirmed by Protein Informatics and Structural Investigations

Background and objectives: Loeys–Dietz syndrome 3, also known as aneurysms-–osteoarthritis syndrome, is an autosomal dominant genetic connective tissue disease caused by pathogenic variants in SMAD3, a transcription factor involved in TGF-β signaling. This disorder is characterized by early-onset osteoarthritis and arterial aneurysms. Common features include scoliosis, uvula abnormalities, striae, and velvety skin. Materials and Methods: The pathogenicity of a variant of uncertain significance in the SMAD3 gene was evaluated (variant c.220C > T) through personalized protein informatics and molecular studies. Results: The case of a 44-year-old male, who was originally presumed to have Marfan syndrome, is presented. An expanded gene panel determined the probable cause to be a variant in SMAD3, c.220C > T (p.R74W). His case was complicated by a history of stroke, but his phenotype was otherwise characteristic for Loeys–Dietz syndrome 3. Conclusion: This case emphasizes the importance of comprehensive genetic testing to evaluate patients for connective tissue disorders, as well as the potential benefit of utilizing a protein informatics platform for the assessment of variant pathogenicity.

Disulfide engineering of human Kunitz-type serine protease inhibitors enhances proteolytic stability and target affinity toward mesotrypsin

Serine protease inhibitors of the Kunitz-bovine pancreatic trypsin inhibitor (BPTI) family are ubiquitous biological regulators of proteolysis. These small proteins are resistant to proteolysis, but can be slowly cleaved within the protease-binding loop by target proteases, thereby compromising their activity. For the human protease mesotrypsin, this cleavage is especially rapid. Here, we aimed to stabilize the Kunitz domain structure against proteolysis through disulfide engineering. Substitution within the Kunitz inhibitor domain of the amyloid precursor protein (APPI) that incorporated a new disulfide bond between residues 17 and 34 reduced proteolysis by mesotrypsin 74-fold. Similar disulfide engineering of tissue factor pathway inhibitor-1 Kunitz domain 1 (^KD1TFPI1) and bikunin Kunitz domain 2 (^KD2bikunin) likewise stabilized these inhibitors against mesotrypsin proteolysis 17- and 6.6-fold, respectively. Crystal structures of disulfide-engineered APPI and ^KD1TFPI1 variants in a complex with mesotrypsin at 1.5 and 2.0 Å resolution, respectively, confirmed the formation of well-ordered disulfide bonds positioned to stabilize the binding loop. Long all-atom molecular dynamics simulations of disulfide-engineered Kunitz domains and their complexes with mesotrypsin revealed conformational stabilization of the primed side of the inhibitor-binding loop by the engineered disulfide, along with global suppression of conformational dynamics in the Kunitz domain. Our findings suggest that the Cys-17-Cys-34 disulfide slows proteolysis by dampening conformational fluctuations in the binding loop and minimizing motion at the enzyme-inhibitor interface. The generalizable approach developed here for the stabilization against proteolysis of Kunitz domains, which can serve as important scaffolds for therapeutics, may thus find applications in drug development.

Protein informatics combined with multiple data sources enriches the clinical characterization of novel TRPV4 variant causing an intermediate skeletal dysplasia

Background

Transient receptor potential cation channel subfamily V member 4 (TRPV4) is an ion channel permeable to Ca²⁺ that is sensitive to physical, hormonal, and chemical stimuli. This protein is expressed in many cell types, including osteoclasts, chondrocytes, and sensory neurons. As such, pathogenic variants of this gene are associated with skeletal dysplasias and neuromuscular disorders. Pathogenesis of these phenotypes is not yet completely understood, but it is known that genotype–phenotype correlations for TRPV4 pathogenic variants often are not present.

Methods

Newly characterized, suspected pathogenic variant in TRPV4 was analyzed using protein informatics and personalized protein‐level molecular studies, genomic exome analysis, and clinical study.

Results

This statement is demonstrated in the family of our proband, a 47‐year‐old female having the novel c.2401A>G (p.K801E) variant of TRPV4. We discuss the common symptoms between the proband, her father, and her daughter, and compare her phenotype to known TRPV4‐associated skeletal dysplasias.

Conclusions

Protein informatics and molecular modeling are used to confirm the pathogenicity of the unique TRPV4 variant found in this family. Multiple data were combined in a comprehensive manner to give complete overall perspective on the patient disease and prognosis.

Personalized molecular modeling for pinpointing associations of protein dysfunction and variants associated with hereditary cancer syndromes

Background

Although the process of reclassification of a variant of uncertain significance can be complex, they are commonly detected through molecular testing. It often takes years before enough clinical data are acquired, and it can be costly and time‐consuming to perform functional analysis of a single variant. It is important that other tools are developed to aid in clarifying how a specific genetic variant impacts a protein's function, and ultimately the health of the patient.

Methods

Two more newly characterized, suspected pathogenic variants in NBN and PTEN were analyzed through personalized protein modeling. Comparisons between the wild‐type and altered protein were studied using simulations, genomic exome analysis, and clinic study.

Results

Modeling of the new NBN and PTEN protein structures suggested loss of essential domains important for normal enzymatic function for these personalized genomic examples which matched the clinical findings.

Conclusion

The defects detected through modeling were consistent with the expected clinical effect. Personalized protein modeling is another tool for determination of correct variant classification, which can become further useful through construction of deposition archive.

Protein molecular modeling techniques investigating novel TAB2 variant R347X causing cardiomyopathy and congenital heart defects in multigenerational family

BACKGROUND

Haploinsufficiency of TAB2 is known to cause congenital heart defects and cardiomyopathy due to its important roles in cardiovascular tissue, both during development and through adult life. We report a sibling pair displaying adult-onset cardiomyopathy, hypermobility, and mild myopia. Our proband, a 39-year-old male, presents only with the above symptoms, while his 36-year-old sister was also notable for a ventricular septal defect in her infancy.

METHODS

Whole-exome sequencing was utilized to identify the molecular basis of the phenotype found in two siblings. A molecular modeling technique that takes advantage of conformational sampling advances (Maxwell's demon molecular dynamics and Monte Carlo) were used to make a model of the mutant variant for comparative analytics to the wild-type.

RESULTS

Exome sequencing revealed a novel, heterogeneous pathogenic variant in TAB2, c.1039 C>T (p.R347X), that was present in both individuals. This pathogenic variant removes just over half the residues from the TAB2 protein and severely impacts its functional ability, which we describe in detail.

CONCLUSIONS

Analysis of the proband's family showed a history of cardiomyopathy, but no congenital heart defects or connective tissue disease. We highlight the heterogeneity in phenotype of TAB2 pathogenic variants and confirm the pathogenicity of this new variant through neoteric protein modeling techniques.

Protein molecular modeling shows residue T599 is critical to wild-type function of POLG and description of a novel variant associated with the SANDO phenotype

Sensory ataxic neuropathy with dysarthria and ophthalmoparesis (SANDO) is a rare phenotype resulting from pathogenic variants of mitochondrial DNA polymerase gamma (POLG). We modeled a novel POLG variant, T599P, that causes the SANDO phenotype and another variant at the same residue, p.T599E, to observe their effect on protein function and confirm the pathogenicity of T599P. Through neoteric molecular modeling techniques, we show that changes at the T599 residue position introduce extra rigidity into the surrounding helix–loop–helix, which places steric pressure on nearby nucleotides. We also provide a clinical description of the T599P variant, which was found in a 42-year-old female proband. The proband presented a 1-year history of progressive gait instability, dysarthria and foot numbness. Her neurologic examination revealed ataxic dysarthria, restricted eye movements, head and palatal tremors, reduced lower limb reflexes, distal multimodal sensory loss and a wide, unsteady ataxic gait. Electromyography studies indicated a sensory neuropathy. Whole-exome sequencing was pursued after tests for infectious, inflammatory and paraneoplastic causes were negative.

Whole Exome Sequencing and Molecular Modeling of a Missense Variant in TNFAIP3 That Segregates with Disease in a Family with Chronic Urticaria and Angioedema

Chronic urticaria is a common condition characterized by recurrent hives lasting several weeks or months and is usually idiopathic. Approximately half of the individuals with chronic urticaria will present with episodes of angioedema that can be severe and debilitating. In this report, we describe a 47-year-old Hispanic male who presented initially for an evaluation of chronic hives following hospitalization due to hive-induced anaphylaxis. The individual had a history significant for urticaria and angioedema beginning in his early 30s. Interestingly, both the individual's 41-year-old sister and 12-year-old daughter were also affected with chronic urticaria and severe angioedema. Whole exome sequencing of the proband and several family members revealed a heterozygous variant of uncertain significance in exon 2 of TNFAIP3, denoted as c.65G>A (p.R22Q), in all affected members. Variants in TNFAIP3 have been associated with multiple autoimmune diseases, susceptibility to allergy and asthma, and periodic fever syndromes, suggesting that this variant could potentially play a role in disease.

Multifaceted peptide assisted one-pot synthesis of gold nanoparticles for plectin-1 targeted gemcitabine delivery in pancreatic cancer

An astute modification of the plectin-1-targeting peptide KTLLPTP by introducing a C-terminal cysteine preceded by a tyrosine residue imparted a reducing property to the peptide. This novel property is then exploited to fabricate gold nanoparticles (GNP) via an in situ reduction of gold(iii) chloride in a one-pot, green synthesis. The modified peptide KTLLPTPYC also acts as a template to generate highly monodispersed, spherical GNPs with a narrow size distribution and improved stability. Plectin-1 is known to be aberrantly expressed in the surface of pancreatic ductal adenocarcinoma (PDAC) cells while showing cytoplasmic expression in normal cells. The synthesized GNPs are thus in situ surface modified with the peptides via the cysteine residue leaving the N-terminal KTLLPTP sequence free for targeting plectin-1. The visual molecular dynamics based simulations support the experimental observations like particle size, gemcitabine conjugation and architecture of the peptide-grafted nanoassembly. Additionally, GNPs conjugated to gemcitabine demonstrate significantly higher cytotoxicity in vitro in two established PDAC cell lines (AsPC-1 and PANC-1) and an admirable in vivo antitumor efficacy in a PANC-1 orthotopic xenograft model through selective uptake in PDAC tumor tissues. Altogether, this strategy represents a unique method for the fabrication of a GNP based targeted drug delivery platform using a multifaceted peptide that acts as reducing agent, template for GNP synthesis and targeting agent to display remarkable selectivity towards PDAC.

The PINK1 p.I368N mutation affects protein stability and ubiquitin kinase activity

Background

Mutations in PINK1 and PARKIN are the most common causes of recessive early-onset Parkinson’s disease (EOPD). Together, the mitochondrial ubiquitin (Ub) kinase PINK1 and the cytosolic E3 Ub ligase PARKIN direct a complex regulated, sequential mitochondrial quality control. Thereby, damaged mitochondria are identified and targeted to degradation in order to prevent their accumulation and eventually cell death. Homozygous or compound heterozygous loss of either gene function disrupts this protective pathway, though at different steps and by distinct mechanisms. While structure and function of PARKIN variants have been well studied, PINK1 mutations remain poorly characterized, in particular under endogenous conditions. A better understanding of the exact molecular pathogenic mechanisms underlying the pathogenicity is crucial for rational drug design in the future.

Methods

Here, we characterized the pathogenicity of the PINK1 p.I368N mutation on the clinical and genetic as well as on the structural and functional level in patients’ fibroblasts and in cell-based, biochemical assays.

Results

Under endogenous conditions, PINK1 p.I368N is expressed, imported, and N-terminally processed in healthy mitochondria similar to PINK1 wild type (WT). Upon mitochondrial damage, however, full-length PINK1 p.I368N is not sufficiently stabilized on the outer mitochondrial membrane (OMM) resulting in loss of mitochondrial quality control. We found that binding of PINK1 p.I368N to the co-chaperone complex HSP90/CDC37 is reduced and stress-induced interaction with TOM40 of the mitochondrial protein import machinery is abolished. Analysis of a structural PINK1 p.I368N model additionally suggested impairments of Ub kinase activity as the ATP-binding pocket was found deformed and the substrate Ub was slightly misaligned within the active site of the kinase. Functional assays confirmed the lack of Ub kinase activity.

Conclusions

Here we demonstrated that mutant PINK1 p.I368N can not be stabilized on the OMM upon mitochondrial stress and due to conformational changes in the active site does not exert kinase activity towards Ub. In patients’ fibroblasts, biochemical assays and by structural analyses, we unraveled two pathomechanisms that lead to loss of function upon mutation of p.I368N and highlight potential strategies for future drug development.

Electronic supplementary material

The online version of this article (doi:10.1186/s13024-017-0174-z) contains supplementary material, which is available to authorized users.

An Acrobatic Substrate Metamorphosis Reveals a Requirement for Substrate Conformational Dynamics in Trypsin Proteolysis

The molecular basis of enzyme catalytic power and specificity derives from dynamic interactions between enzyme and substrate during catalysis. Although considerable effort has been devoted to understanding how conformational dynamics within enzymes affect catalysis, the role of conformational dynamics within protein substrates has not been addressed. Here, we examine the importance of substrate dynamics in the cleavage of Kunitz-bovine pancreatic trypsin inhibitor protease inhibitors by mesotrypsin, finding that the varied conformational dynamics of structurally similar substrates can profoundly impact the rate of catalysis. A 1.4-Å crystal structure of a mesotrypsin-product complex formed with a rapidly cleaved substrate reveals a dramatic conformational change in the substrate upon proteolysis. By using long all-atom molecular dynamics simulations of acyl-enzyme intermediates with proteolysis rates spanning 3 orders of magnitude, we identify global and local dynamic features of substrates on the nanosecond-microsecond time scale that correlate with enzymatic rates and explain differential susceptibility to proteolysis. By integrating multiple enhanced sampling methods for molecular dynamics, we model a viable conformational pathway between substrate-like and product-like states, linking substrate dynamics on the nanosecond-microsecond time scale with large collective substrate motions on the much slower time scale of catalysis. Our findings implicate substrate flexibility as a critical determinant of catalysis.

Heterozygous PINK1 p.G411S increases risk of Parkinson's disease via a dominant-negative mechanism

See Gandhi and Plun-Favreau (doi:10.1093/aww320) for a scientific commentary on this article.

Heterozygous mutations in recessive Parkinson’s disease genes have been postulated to increase disease risk. Puschmann et al. report a genetic association between heterozygous PINK1 p.G411S and Parkinson’s disease. They provide structural and functional explanations for a partial dominant-negative effect of the mutant protein, which impairs wild-type PINK1 activity through hetero-dimerization.

See Gandhi and Plun-Favreau (doi:10.1093/aww320) for a scientific commentary on this article.

It has been postulated that heterozygous mutations in recessive Parkinson’s genes may increase the risk of developing the disease. In particular, the PTEN-induced putative kinase 1 (PINK1) p.G411S (c.1231G>A, rs45478900) mutation has been reported in families with dominant inheritance patterns of Parkinson’s disease, suggesting that it might confer a sizeable disease risk when present on only one allele. We examined families with PINK1 p.G411S and conducted a genetic association study with 2560 patients with Parkinson’s disease and 2145 control subjects. Heterozygous PINK1 p.G411S mutations markedly increased Parkinson’s disease risk (odds ratio = 2.92, P = 0.032); significance remained when supplementing with results from previous studies on 4437 additional subjects (odds ratio = 2.89, P = 0.027). We analysed primary human skin fibroblasts and induced neurons from heterozygous PINK1 p.G411S carriers compared to PINK1 p.Q456X heterozygotes and PINK1 wild-type controls under endogenous conditions. While cells from PINK1 p.Q456X heterozygotes showed reduced levels of PINK1 protein and decreased initial kinase activity upon mitochondrial damage, stress-response was largely unaffected over time, as expected for a recessive loss-of-function mutation. By contrast, PINK1 p.G411S heterozygotes showed no decrease of PINK1 protein levels but a sustained, significant reduction in kinase activity. Molecular modelling and dynamics simulations as well as multiple functional assays revealed that the p.G411S mutation interferes with ubiquitin phosphorylation by wild-type PINK1 in a heterodimeric complex. This impairs the protective functions of the PINK1/parkin-mediated mitochondrial quality control. Based on genetic and clinical evaluation as well as functional and structural characterization, we established p.G411S as a rare genetic risk factor with a relatively large effect size conferred by a partial dominant-negative function phenotype.

Activation of the E3 ubiquitin ligase Parkin

The PINK1 (phosphatase and tensin homologue-induced putative kinase 1)/Parkin-dependent mitochondrial quality control pathway mediates the clearance of damaged organelles, but appears to be disrupted in Parkinson's disease (PD) [Springer and Kahle (2011) Autophagy 7, 266-278]. Upon mitochondrial stress, PINK1 activates the E3 ubiquitin (Ub) ligase Parkin through phosphorylation of the Ub-like (UBL) domain of Parkin and of the small modifier Ub itself at a conserved residue [Sauvé and Gehring (2014) Cell Res. 24, 1025-1026]. Recently resolved partial crystal structures of Parkin showed a 'closed', auto-inhibited conformation, consistent with its notoriously weak enzymatic activity at steady state [Wauer and Komander (2013) EMBO J. 32, 2099-2112; Riley et al. (2013) Nat. Commun. 4, 1982; Trempe et al. (2013) Science 340, 1451-1455; Spratt et al. (2013) Nat. Commun. 4, 1983]. It has thus become clear that Parkin must undergo major structural rearrangements in order to unleash its catalytic functions. Recent published findings derived from X-ray structures and molecular modelling present a complete structural model of human Parkin at an all-atom resolution [Caulfield et al. (2014) PLoS Comput. Biol. 10, e1003935]. The results of the combined in silico simulations-based and experimental assay-based study indicates that PINK1-dependent Ser65 phosphorylation of Parkin is required for its activation and triggering of 'opening' conformations. Indeed, the obtained structures showed a sequential release of Parkin's intertwined domains and allowed docking of an Ub-charged E2 coenzyme, which could enable its enzymatic activity. In addition, using cell-based screening, select E2 enzymes that redundantly, cooperatively or antagonistically regulate Parkin's activation and/or enzymatic functions at different stages of the mitochondrial autophagy (mitophagy) process were identified [Fiesel et al. (2014) J. Cell Sci. 127, 3488-3504]. Other work that aims to pin-point the particular pathogenic dysfunctions of Parkin mis-sense mutations have been recently disseminated (Fabienne C. Fiesel, Thomas R. Caulfield, Elisabeth L. Moussaud-Lamodiere, Daniel F.A.R. Dourado, Kotaro Ogaki, Owen A. Ross, Samuel C. Flores, and Wolfdieter Springer, submitted). Such a structure-function approach provides the basis for the dissection of Parkin's regulation and a targeted drug design to identify small-molecule activators of this neuroprotective E3 Ub ligase.

(Patho-)physiological relevance of PINK1-dependent ubiquitin phosphorylation

Mutations in PINK1 and PARKIN cause recessive, early-onset Parkinson's disease (PD). Together, these two proteins orchestrate a protective mitophagic response that ensures the safe disposal of damaged mitochondria. The kinase PINK1 phosphorylates ubiquitin (Ub) at the conserved residue S65, in addition to modifying the E3 ubiquitin ligase Parkin. The structural and functional consequences of Ub phosphorylation (pS65-Ub) have already been suggested from in vitro experiments, but its (patho-)physiological significance remains unknown. We have generated novel antibodies and assessed pS65-Ub signals in vitro and in cells, including primary neurons, under endogenous conditions. pS65-Ub is dependent on PINK1 kinase activity as confirmed in patient fibroblasts and postmortem brain samples harboring pathogenic mutations. We show that pS65-Ub is reversible and barely detectable under basal conditions, but rapidly induced upon mitochondrial stress in cells and amplified in the presence of functional Parkin. pS65-Ub accumulates in human brain during aging and disease in the form of cytoplasmic granules that partially overlap with mitochondrial, lysosomal, and total Ub markers. Additional studies are now warranted to further elucidate pS65-Ub functions and fully explore its potential for biomarker or therapeutic development.

Structural and Functional Impact of Parkinson Disease-Associated Mutations in the E3 Ubiquitin Ligase Parkin

Mutations in the PARKIN/PARK2 gene that result in loss-of-function of the encoded, neuroprotective E3 ubiquitin ligase Parkin cause recessive, familial early-onset Parkinson disease. As an increasing number of rare Parkin sequence variants with unclear pathogenicity are identified, structure-function analyses will be critical to determine their disease relevance. Depending on the specific amino acids affected, several distinct pathomechanisms can result in loss of Parkin function. These include disruption of overall Parkin folding, decreased solubility, and protein aggregation. However pathogenic effects can also result from misregulation of Parkin autoinhibition and of its enzymatic functions. In addition, interference of binding to coenzymes, substrates, and adaptor proteins can affect its catalytic activity too. Herein, we have performed a comprehensive structural and functional analysis of 21 PARK2 missense mutations distributed across the individual protein domains. Using this combined approach, we were able to pinpoint some of the pathogenic mechanisms of individual sequence variants. Similar analyses will be critical in gaining a complete understanding of the complex regulations and enzymatic functions of Parkin. These studies will not only highlight the important residues, but will also help to develop novel therapeutics aimed at activating and preserving an active, neuroprotective form of Parkin.

Structural Insights into tRNA Dynamics on the Ribosome

High-resolution structures at different stages, as well as biochemical, single molecule and computational approaches have highlighted the elasticity of tRNA molecules when bound to the ribosome. It is well acknowledged that the inherent structural flexibility of the tRNA lies at the heart of the protein synthesis process. Here, we review the recent advances and describe considerations that the conformational changes of the tRNA molecules offer about the mechanisms grounded in translation.

Nimbolide targets BCL2 and induces apoptosis in preclinical models of Waldenströms macroglobulinemia

Neem leaf extract (NLE) has medicinal properties, which have been attributed to its limonoid content. We identified the NLE tetranorterpenoid, nimbolide, as being the key limonoid responsible for the cytotoxicity of NLE in various preclinical models of human B-lymphocyte cancer. Of the models tested, Waldenströms macroglobulinemia (WM) cells were most sensitive to nimbolide, undergoing significant mitochondrial mediated apoptosis. Notably, nimbolide toxicity was also observed in drug-resistant (bortezomib or ibrutinib) WM cells. To identify putative targets of nimbolide, relevant in WM, we used chemoinformatics-based approaches comprised of virtual in silico screening, molecular modeling and target–ligand reverse docking. In silico analysis revealed the antiapoptotic protein BCL2 was the preferential binding partner of nimbolide. The significance of this finding was further tested in vitro in RS4;11 (BCL2-dependent) tumor cells, in which nimbolide induced significantly more apoptosis compared with BCL2 mutated (Jurkat BCL2^Ser70-Ala) cells. Lastly, intraperitoneal administration of nimbolide in WM tumor xenografted mice, significantly reduced tumor growth and IgM secretion in vivo, while modulating the expression of several proteins as seen on immunohistochemistry. Overall, our data demonstrate that nimbolide is highly active in WM cells, as well as other B-cell cancers, and engages BCL2 to exert its cytotoxic activity.

Phosphorylation by PINK1 releases the UBL domain and initializes the conformational opening of the E3 ubiquitin ligase Parkin

Loss-of-function mutations in PINK1 or PARKIN are the most common causes of autosomal recessive Parkinson's disease. Both gene products, the Ser/Thr kinase PINK1 and the E3 Ubiquitin ligase Parkin, functionally cooperate in a mitochondrial quality control pathway. Upon stress, PINK1 activates Parkin and enables its translocation to and ubiquitination of damaged mitochondria to facilitate their clearance from the cell. Though PINK1-dependent phosphorylation of Ser65 is an important initial step, the molecular mechanisms underlying the activation of Parkin's enzymatic functions remain unclear. Using molecular modeling, we generated a complete structural model of human Parkin at all atom resolution. At steady state, the Ub ligase is maintained inactive in a closed, auto-inhibited conformation that results from intra-molecular interactions. Evidently, Parkin has to undergo major structural rearrangements in order to unleash its catalytic activity. As a spark, we have modeled PINK1-dependent Ser65 phosphorylation in silico and provide the first molecular dynamics simulation of Parkin conformations along a sequential unfolding pathway that could release its intertwined domains and enable its catalytic activity. We combined free (unbiased) molecular dynamics simulation, Monte Carlo algorithms, and minimal-biasing methods with cell-based high content imaging and biochemical assays. Phosphorylation of Ser65 results in widening of a newly defined cleft and dissociation of the regulatory N-terminal UBL domain. This motion propagates through further opening conformations that allow binding of an Ub-loaded E2 co-enzyme. Subsequent spatial reorientation of the catalytic centers of both enzymes might facilitate the transfer of the Ub moiety to charge Parkin. Our structure-function study provides the basis to elucidate regulatory mechanisms and activity of the neuroprotective Parkin. This may open up new avenues for the development of small molecule Parkin activators through targeted drug design.

Parkinson's disease (PD) is a devastating neurological condition caused by the selective and progressive degeneration of dopaminergic neurons in the brain. Loss-of-function mutations in the PINK1 or PARKIN genes are the most common causes of recessively inherited PD. Together the encoded proteins coordinate a protective cellular quality control pathway that allows elimination of impaired mitochondria in order to prevent further cellular damage and ultimately death. Although it is known that the kinase PINK1 operates upstream and activates the E3 Ubiquitin ligase Parkin, the molecular mechanisms remain elusive. Here, we combined state-of-the art computational and functional biological methods to demonstrate that Parkin is sequentially activated through PINK1-dependent phosphorylation and subsequent structural rearrangement. The induced motions result in release of Parkin's closed, auto-inhibited conformation to liberate its enzymatic functions. We provide for the first time a complete protein structure of Parkin at an all atom resolution and a comprehensive molecular dynamics simulation of its activation and opening conformations. The generated models will allow uncovering the exact mechanisms of regulation and enzymatic activity of Parkin and potentially the development of novel therapeutics through a structure-function-based drug design.

Local structural and environmental factors define the efficiency of an RNA pseudoknot involved in programmed ribosomal frameshift process

In programmed -1 ribosomal frameshift, an RNA pseudoknot stalls the ribosome at specific sequence and restarts translation in a new reading frame. A precise understanding of structural characteristics of these pseudoknots and their PRF inducing ability has not been clear to date. To investigate this phenomenon, we have studied various structural aspects of a -1 PRF inducing RNA pseudoknot from BWYV using extensive molecular dynamics simulations. A set of functional and poorly functional forms, for which previous mutational data were available, were chosen for analysis. These structures differ from each other by either single base substitutions or base-pair replacements from the native structure. We have rationalized how certain mutations in RNA pseudoknot affect its function; e.g., a specific base substitution in loop 2 stabilizes the junction geometry by forming multiple noncanonical hydrogen bonds, leading to a highly rigid structure that could effectively resist ribosome-induced unfolding, thereby increasing efficiency. While, a CG to AU pair substitution in stem 1 leads to loss of noncanonical hydrogen bonds between stems and loop, resulting in a less stable structure and reduced PRF inducing ability, inversion of a pair in stem 2 alters specific base-pair geometry that might be required in ribosomal recognition of nucleobase groups, negatively affecting pseudoknot functioning. These observations illustrate that the ability of an RNA pseudoknot to induce -1 PRF with an optimal rate depends on several independent factors that contribute to either the local conformational variability or geometry.

Structures of yeast 80S ribosome-tRNA complexes in the rotated and nonrotated conformations

The structural understanding of eukaryotic translation lags behind that of translation on bacterial ribosomes. Here, we present two subnanometer resolution structures of S. cerevisiae 80S ribosome complexes formed with either one or two tRNAs and bound in response to an mRNA fragment containing the Kozak consensus sequence. The ribosomes adopt two globally different conformations that are related to each other by the rotation of the small subunit. Comparison with bacterial ribosome complexes reveals that the global structures and modes of intersubunit rotation of the yeast ribosome differ significantly from those in the bacterial counterpart, most notably in the regions involving the tRNA, small ribosomal subunit, and conserved helix 69 of the large ribosomal subunit. The structures provide insight into ribosome dynamics implicated in tRNA translocation and help elucidate the role of the Kozak fragment in positioning an open reading frame during translation initiation in eukaryotes.

Motion of transfer RNA from the A/T state into the A-site using docking and simulations

The ribosome catalyzes peptidyl transfer reactions at the growing nascent polypeptide chain. Here, we present a structural mechanism for selecting cognate over near-cognate A/T transfer RNA (tRNA). In part, the structural basis for the fidelity of translation relies on accommodation to filter cognate from near-cognate tRNAs. To examine the assembly of tRNAs within the ribonucleic-riboprotein complex, we conducted a series of all-atom molecular dynamics (MD) simulations of the entire solvated 70S Escherichia coli ribosome, along with its associated cofactors, proteins, and messenger RNA (mRNA). We measured the motion of the A/T state of tRNA between initial binding and full accommodation. The mechanism of rejection was investigated. Using novel in-house algorithms, we determined trajectory pathways. Despite the large intersubunit cavity, the available space is limited by the presence of the tRNA, which is equally large. This article describes a "structural gate," formed between helices 71 and 92 on the ribosomal large subunit, which restricts tRNA motion. The gate and the interacting protein, L14, of the 50S ribosome act as steric filters in two consecutive substeps during accommodation, each requiring: (1) sufficient energy contained in the hybrid tRNA kink and (2) sufficient energy in the Watson-Crick base pairing of the codon-anticodon. We show that these barriers act to filter out near-cognate tRNA and promote proofreading of the codon-anticodon. Since proofreading is essential for understanding the fidelity of translation, our model for the dynamics of this process has substantial biomedical implications.

The concept of RNA-assisted protein folding: the role of tRNA

We suggest that tRNA actively participates in the transfer of 3D information from mRNA to peptides--in addition to its well-known, "classical" role of translating the 3-letter RNA codes into the one letter protein code. The tRNA molecule displays a series of thermodynamically favored configurations during translation, a movement which places the codon and coded amino acids in proximity to each other and make physical contact between some amino acids and their codons possible. This specific codon-amino acid interaction of some selected amino acids is necessary for the transfer of spatial information from mRNA to coded proteins, and is known as RNA-assisted protein folding.

Computational studies of molecular machines: the ribosome

The past decade has produced an avalanche of experimental data on the structure and dynamics of the ribosome. Groundbreaking studies in structural biology and kinetics have placed important constraints on ribosome structural dynamics. However, a gulf remains between static structures and time dependent data. In particular, X-ray crystallography and cryo-EM studies produce static models of the ribosome in various states, but lack dynamic information. Single molecule studies produce information on the rates of transitions between these states but do not have high-resolution spatial information. Computational studies have aided in bridging this gap by providing atomic resolution simulations of structural fluctuations and transitions between configurations.

Inter-ring rotation of apolipoprotein A-I protein monomers for the double-belt model using biased molecular dynamics

The double belt model for lipid-bound discoidal apolipoprotein A-I consists of two alpha-helical monomers bound about an unilamellar bilayer of lipids. Previous work, based on salt bridge calculations, has demonstrated that the L5/5 registration, Milano mutant, and Paris mutant are preferred conformations for apolipoprotein A-I. The salt bridge scoring indicated better energetic scoring in these alignments. The Paris (R151C) and Milano (R173C) mutants indicate a mode of change must be available. To find proper registration, one proposed change is a 'rotationally' independent circular motion of the two protein monomers about the lipid unilamellar bilayer core. Here, we present computational data for independent inter-ring rotation of the two alpha-helical monomers about the lipid unilamellar bilayer core. The simulations presented here support the existing double-belt model. We find the rotation of the two protein monomers is able to occur with biasing. We determine that a cysteine mutant at Glu107 as a possible target for future mutational studies. Since HDL remodeling is necessary for cholesterol transport, our model for remodeling through dynamics has substantial biomedical implications.