Formation and characterization of BMP2/GDF5 and BMP4/GDF5 heterodimers.

BACKGROUND: Proteins of the TGFß family, which are largely studied as homodimers, are also known to form heterodimers with biological activity distinct from their component homodimers. For instance, heterodimers of bone morphogenetic proteins, including BMP2/BMP7, BMP2/BMP6, and BMP9/BMP10, among others, have illustrated the importance of these heterodimeric proteins within the context of TGFß signaling. RESULTS: In this study, we have determined that mature GDF5 can be combined with mature BMP2 or BMP4 to form BMP2/GDF5 and BMP4/GDF5 heterodimer. Intriguingly, this combination of a BMP2 or BMP4 monomer, which exhibit high affinity to heparan sulfate characteristic to the BMP class, with a GDF5 monomer with low heparan sulfate affinity produces a heterodimer with an intermediate affinity. Using heparin affinity chromatography to purify the heterodimeric proteins, we then determined that both the BMP2/GDF5 and BMP4/GDF5 heterodimers consistently signaled potently across an array of cellular and in vivo systems, while the activities of their homodimeric counterparts were more context dependent. These differences were likely driven by an increase in the combined affinities for the type 1 receptors, Alk3 and Alk6. Furthermore, the X-ray crystal structure of BMP2/GDF5 heterodimer was determined, highlighting the formation of two asymmetric type 1 receptor binding sites that are both unique relative to the homodimers. CONCLUSIONS: Ultimately, this method of heterodimer production yielded a signaling molecule with unique properties relative to the homodimeric ligands, including high affinity to multiple type 1 and moderate heparan binding affinity.

Networks of necessity: Simulating COVID-19 mitigation strategies for disabled people and their caregivers.

A major strategy to prevent the spread of COVID-19 is the limiting of in-person contacts. However, limiting contacts is impractical or impossible for the many disabled people who do not live in care facilities but still require caregivers to assist them with activities of daily living. We seek to determine which interventions can best prevent infections of disabled people and their caregivers. To accomplish this, we simulate COVID-19 transmission with a compartmental model that includes susceptible, exposed, asymptomatic, symptomatically ill, hospitalized, and removed/recovered individuals. The networks on which we simulate disease spread incorporate heterogeneity in the risk levels of different types of interactions, time-dependent lockdown and reopening measures, and interaction distributions for four different groups (caregivers, disabled people, essential workers, and the general population). Of these groups, we find that the probability of becoming infected is largest for caregivers and second largest for disabled people. Consistent with this finding, our analysis of network structure illustrates that caregivers have the largest modal eigenvector centrality of the four groups. We find that two interventions-contact-limiting by all groups and mask-wearing by disabled people and caregivers-most reduce the number of infections in disabled and caregiver populations. We also test which group of people spreads COVID-19 most readily by seeding infections in a subset of each group and comparing the total number of infections as the disease spreads. We find that caregivers are the most potent spreaders of COVID-19, particularly to other caregivers and to disabled people. We test where to use limited infection-blocking vaccine doses most effectively and find that (1) vaccinating caregivers better protects disabled people from infection than vaccinating the general population or essential workers and that (2) vaccinating caregivers protects disabled people from infection about as effectively as vaccinating disabled people themselves. Our results highlight the potential effectiveness of mask-wearing, contact-limiting throughout society, and strategic vaccination for limiting the exposure of disabled people and their caregivers to COVID-19.

Semi-automated curation of metabolic models via flux balance analysis: a case study with Mycoplasma gallisepticum.

Primarily used for metabolic engineering and synthetic biology, genome-scale metabolic modeling shows tremendous potential as a tool for fundamental research and curation of metabolism. Through a novel integration of flux balance analysis and genetic algorithms, a strategy to curate metabolic networks and facilitate identification of metabolic pathways that may not be directly inferable solely from genome annotation was developed. Specifically, metabolites involved in unknown reactions can be determined, and potentially erroneous pathways can be identified. The procedure developed allows for new fundamental insight into metabolism, as well as acting as a semi-automated curation methodology for genome-scale metabolic modeling. To validate the methodology, a genome-scale metabolic model for the bacterium Mycoplasma gallisepticum was created. Several reactions not predicted by the genome annotation were postulated and validated via the literature. The model predicted an average growth rate of 0.358±0.12[Formula: see text], closely matching the experimentally determined growth rate of M. gallisepticum of 0.244±0.03[Formula: see text]. This work presents a powerful algorithm for facilitating the identification and curation of previously known and new metabolic pathways, as well as presenting the first genome-scale reconstruction of M. gallisepticum.

A fine kinetic balance of interactions directs transcription factor hubs to genes.

Eukaryotic gene regulation relies on the binding of sequence-specific transcription factors (TFs). TFs bind chromatin transiently yet occupy their target sites by forming high-local concentration microenvironments (hubs and condensates) that increase the frequency of binding events. Despite their ubiquity, such microenvironments have been difficult to study in endogenous contexts due to technical limitations. Here, we overcome these limitations and investigate how hubs drive TF occupancy at their targets. Using a DNA binding perturbation to a hub-forming TF, Zelda, in Drosophila embryos, we find that hub properties, including the stability and frequencies of associations to targets, are key determinants of TF occupancy. Our data suggest that the targeting of these hubs is driven not just by specific DNA motif recognition, but also by a fine-tuned kinetic balance of interactions between TFs and their co-binding partners.

EpicTope: narrating protein sequence features to identify non-disruptive epitope tagging sites.

Epitope tagging is an invaluable technique enabling the identification, tracking, and purification of proteins in vivo. We developed a tool, EpicTope, to facilitate this method by identifying amino acid positions suitable for epitope insertion. Our method uses a scoring function that considers multiple protein sequence and structural features to determine locations least disruptive to the protein's function. We validated our approach on the zebrafish Smad5 protein, showing that multiple predicted internally tagged Smad5 proteins rescue zebrafish smad5 mutant embryos, while the N- and C-terminal tagged variants do not, also as predicted. We further show that the internally tagged Smad5 proteins are accessible to antibodies in wholemount zebrafish embryo immunohistochemistry and by western blot. Our work demonstrates that EpicTope is an accessible and effective tool for designing epitope tag insertion sites. EpicTope is available under a GPL-3 license from: https://github.com/FriedbergLab/Epictope.

A Primary Care Provider's Guide to Bone Health in Spinal Cord-Related Paralysis.

Individuals with spinal cord injury/disorder (SCI/D) are at high risk for developing secondary osteoporosis. Bone loss after neurologic injury is multifactorial and is dependent on the time from and extent of neurologic injury. Most bone loss occurs in the first year after complete motor paralysis, and fractures occur most commonly in the distal femur and proximal tibia (paraplegic fracture). The 2019 International Society for Clinical Densitometry Position Statement in SCI establishes that dual-energy X-ray absorptiometry (DXA) can be used to both diagnose osteoporosis and predict lower extremity fracture risk in individuals with SCI/D. Pharmacologic treatments used in primary osteoporosis have mixed results when used for SCI/D-related osteoporosis. Ambulation, standing, and electrical stimulation may be helpful at increasing bone mineral density (BMD) in individuals with SCI/D but do not necessarily correlate with fracture risk reduction. Clinicians caring for individuals with spinal cord-related paralysis must maintain a high index of suspicion for fragility fractures and consider referral for surgical evaluation and management.

Imaging and Quantification of P-Smad1/5 in Zebrafish Blastula and Gastrula Embryos.

Spatiotemporal patterns of morphogen activity drive differential gene expression with a high degree of precision within a developing embryo and reproducibly between embryos. Understanding the formation and function of a morphogen gradient during development requires quantitative measurement of morphogen activity throughout an individual embryo and also between embryos within a population. Quantification of morphogen gradients in to presents unique challenges in imaging and image processing to minimize error and maximize the quality of the data so it may be used in computational models of development and in statistically testing hypotheses. Here we present methods for the preparation, immunostaining, imaging, and quantification of a bone morphogenetic protein (BMP) activity gradient in individual zebrafish embryos as well as methods for acquiring population-level statistics after embryo grouping and alignment. This quantitative approach can be extended to other morphogen systems, and the computational codes can be adapted to other imaging contexts in zebrafish and other organisms.

TGF-ß Family Signaling in Early Vertebrate Development.

TGF-ß family ligands function in inducing and patterning many tissues of the early vertebrate embryonic body plan. Nodal signaling is essential for the specification of mesendodermal tissues and the concurrent cellular movements of gastrulation. Bone morphogenetic protein (BMP) signaling patterns tissues along the dorsal-ventral axis and simultaneously directs the cell movements of convergence and extension. After gastrulation, a second wave of Nodal signaling breaks the symmetry between the left and right sides of the embryo. During these processes, elaborate regulatory feedback between TGF-ß ligands and their antagonists direct the proper specification and patterning of embryonic tissues. In this review, we summarize the current knowledge of the function and regulation of TGF-ß family signaling in these processes. Although we cover principles that are involved in the development of all vertebrate embryos, we focus specifically on three popular model organisms: the mouse Mus musculus, the African clawed frog of the genus Xenopus, and the zebrafish Danio rerio, highlighting the similarities and differences between these species.

Variant BMP receptor mutations causing fibrodysplasia ossificans progressiva (FOP) in humans show BMP ligand-independent receptor activation in zebrafish.

The large majority of cases of the autosomal dominant human disease fibrodysplasia ossificans progressiva (FOP) are caused by gain-of-function Arg206His mutations in the BMP type I receptor ACVR1 (ALK2). The Arg206His mutation is located in the GS domain of the type I receptor. This region is normally phosphorylated by the BMP type II receptor, which activates the type I receptor to phosphorylate its substrate, the signal transducer Smad1/5/8. A small subset of patients with FOP carry variant mutations in ACVR1 altering Gly328 to Trp, Glu or Arg. Since these mutations lie outside the GS domain, the mechanism through which ACVR1 Gly328 mutations cause disease remains unclear. We used a zebrafish embryonic development assay to test the signaling of human ACVR1 Gly328 mutant receptors comparing them to the Arg206His mutant. In this assay increased or decreased BMP pathway activation alters dorsal-ventral axial patterning, providing a sensitive assay for altered BMP signaling levels. We expressed the human ACVR1 Gly328 mutant receptors in zebrafish embryos to investigate their signaling activities. We found that all ACVR1 Gly328 human mutations ventralized wild-type embryos and could partially rescue Bmp7-deficient embryos, indicating that these mutant receptors can activate BMP signaling in a BMP ligand-independent manner. The degree of ventralization or rescue was similar among all three Gly328 mutants. Smad1/5 phosphorylation, a readout of BMP receptor signaling, was mildly increased by ACVR1 Gly328 mutations. Gene expression analyses demonstrate expanded ventral and reciprocal loss of dorsal cell fate markers. This study demonstrates that Gly328 mutants increase receptor activation and BMP ligand-independent signaling through Smad phosphorylation.

Systems biology derived source-sink mechanism of BMP gradient formation.

A morphogen gradient of Bone Morphogenetic Protein (BMP) signaling patterns the dorsoventral embryonic axis of vertebrates and invertebrates. The prevailing view in vertebrates for BMP gradient formation is through a counter-gradient of BMP antagonists, often along with ligand shuttling to generate peak signaling levels. To delineate the mechanism in zebrafish, we precisely quantified the BMP activity gradient in wild-type and mutant embryos and combined these data with a mathematical model-based computational screen to test hypotheses for gradient formation. Our analysis ruled out a BMP shuttling mechanism and a bmp transcriptionally-informed gradient mechanism. Surprisingly, rather than supporting a counter-gradient mechanism, our analyses support a fourth model, a source-sink mechanism, which relies on a restricted BMP antagonist distribution acting as a sink that drives BMP flux dorsally and gradient formation. We measured Bmp2 diffusion and found that it supports the source-sink model, suggesting a new mechanism to shape BMP gradients during development.