A Multiplex Homology-Directed DNA Repair Assay Reveals the Impact of More Than 1,000 BRCA1 Missense Substitution Variants on Protein Function.

Loss-of-function pathogenic variants in BRCA1 confer a predisposition to breast and ovarian cancer. Genetic testing for sequence changes in BRCA1 frequently reveals a missense variant for which the impact on cancer risk and on the molecular function of BRCA1 is unknown. Functional BRCA1 is required for the homology-directed repair (HDR) of double-strand DNA breaks, a critical activity for maintaining genome integrity and tumor suppression. Here, we describe a multiplex HDR reporter assay for concurrently measuring the effects of hundreds of variants of BRCA1 for their role in DNA repair. Using this assay, we characterized the effects of 1,056 amino acid substitutions in the first 192 residues of BRCA1. Benchmarking these results against variants with known effects on DNA repair function or on cancer predisposition, we demonstrate accurate discrimination of loss-of-function versus benign missense variants. We anticipate that this assay can be used to functionally characterize BRCA1 missense variants at scale, even before the variants are observed in results from genetic testing.

Abundance of megalin and Dab2 is reduced in syncytiotrophoblast during placental malaria, which may contribute to low birth weight.

Placental malaria caused by Plasmodium falciparum contributes to ~200,000 child deaths annually, mainly due to low birth weight (LBW). Parasitized erythrocyte sequestration and consequent inflammation in the placenta are common attributes of placental malaria. The precise molecular details of placental changes leading to LBW are still poorly understood. We hypothesized that placental malaria may disturb maternofetal exchange of vitamins, lipids, and hormones mediated by the multi-ligand (n ~ 50) scavenging/signaling receptor megalin, which is abundantly expressed in placenta but was not previously analyzed in pregnancy outcomes. We studied abundance of megalin and its intracellular adaptor protein Dab2 by immunofluorescence microscopy in placental biopsies from Ugandan women with (n = 8) and without (n = 20) active placental malaria. We found that: (a) abundances of both megalin (p = 0.01) and Dab2 (p = 0.006) were significantly reduced in brush border of syncytiotrophoblast of infected placentas; (b) amounts of megalin and Dab2 were strongly correlated (Spearman's r = 0.53, p = 0.003); (c) abundances of megalin and Dab2 (p = 0.046) were reduced in infected placentas from women with LBW deliveries. This study provides first evidence that placental malaria infection is associated with reduced abundance of megalin transport/signaling system and indicate that these changes may contribute to the pathology of LBW.

Massively Parallel Functional Analysis of BRCA1 RING Domain Variants.

Interpreting variants of uncertain significance (VUS) is a central challenge in medical genetics. One approach is to experimentally measure the functional consequences of VUS, but to date this approach has been post hoc and low throughput. Here we use massively parallel assays to measure the effects of nearly 2000 missense substitutions in the RING domain of BRCA1 on its E3 ubiquitin ligase activity and its binding to the BARD1 RING domain. From the resulting scores, we generate a model to predict the capacities of full-length BRCA1 variants to support homology-directed DNA repair, the essential role of BRCA1 in tumor suppression, and show that it outperforms widely used biological-effect prediction algorithms. We envision that massively parallel functional assays may facilitate the prospective interpretation of variants observed in clinical sequencing.

High-throughput screening platform identifies small molecules that prevent sequestration of Plasmodium falciparum-infected erythrocytes.

BACKGROUND: We developed a 2-step approach to screen molecules that prevent and/or reverse Plasmodium falciparum-infected erythrocyte (IE) binding to host receptors. IE adhesion and sequestration in vasculature causes severe malaria, and therefore antiadhesion therapy might be useful as adjunctive treatment. IE adhesion is mediated by the polymorphic family (approximately 60 members) of P. falciparum EMP1 (PfEMP1) multidomain proteins. METHODS: We constructed sets of PfEMP1 domains that bind ICAM-1, CSA, or CD36, receptors that commonly support IE binding. Combinations of domain-coated beads were assayed by Bio-Plex technology as a high-throughput molecular platform to screen antiadhesion molecules (antibodies and small molecules). Molecules identified as so-called hits in the screen (first step) then could be assayed individually for inhibition of binding of live IE to receptors (second step). RESULTS: In proof-of-principle studies, the antiadhesion activity of several antibodies was concordant in Bio-Plex and live IE assays. Using this 2-step approach, we identified several molecules in a small molecule library of 10 000 compounds that could inhibit and reverse binding of IEs to ICAM-1 and CSA receptors. CONCLUSION: This 2-step screening approach should be efficient for identification of antiadhesion drug candidates for falciparum malaria.

Structure-function-immunogenicity studies of PfEMP1 domain DBL2ßPF11_0521, a malaria parasite ligand for ICAM-1.

Plasmodium falciparum virulence has been ascribed to its ability to sequester in deep vascular beds, mediated by the variant surface antigen family PfEMP1 binding endothelial receptors like ICAM-1. We previously observed that naturally-acquired antibodies that block a PfEMP1 domain, DBL2ß of PF11_0521 allele, from binding to the human ICAM1 receptor, reduce the risk of malaria hospitalization in children. Here, we find that DBL2ßPF11_0521 binds ICAM-1 in the low nM range and relate the structure of this domain with its function and immunogenicity. We demonstrate that the interaction with ICAM-1 is not impaired by point mutations in the N-terminal subdomain or in the flexible Loop 4 of DBL2ßPF11_0521, although both substructures were previously implicated in binding ICAM-1. These data will help to refine the existing model of DBLß::ICAM-1 interactions. Antibodies raised against full-length DBL2ßPF11_0521, but not truncated forms lacking the N terminal fragment, block its interaction with ICAM-1. Our data suggest that full length domain is optimal for displaying functional epitopes and has a broad surface of interaction with ICAM-1 that is not disrupted by individual amino acid substitutions at putative key residues. This information might be important for the future design of anti-malarial vaccines based on PfEMP1 antigens.

Automated Event Detection and Activity Monitoring in Long Molecular Dynamics Simulations.

Events of scientific interest in molecular dynamics (MD) simulations, including conformational changes, folding transitions, and translocations of ligands and reaction products, often correspond to high-level structural rearrangements that alter contacts between molecules or among different parts of a molecule. Due to advances in computer architecture and software, MD trajectories representing such structure-changing events have become easier to generate, but the length of these trajectories poses a challenge to scientific interpretation and analysis. In this paper, we present automated methods for the detection of potentially important structure-changing events in long MD trajectories. In contrast with traditional tools for the analysis of such trajectories, our methods provide a detailed report of broken and formed contacts that aids in the identification of specific time-dependent side-chain interactions. Our approach employs a coarse-grained representation of amino acid side chains, a contact metric based on higher order generalizations of Delaunay tetrahedralization, techniques for detecting significant shifts in the resulting contact time series, and a new kernel-based measure of contact alteration activity. The analysis methods we describe are incorporated in a newly developed package, called TimeScapes, which is freely available and compatible with trajectories generated by a variety of popular MD programs. Tests based on actual microsecond time scale simulations demonstrate that the package can be used to efficiently detect and characterize important conformational changes in realistic protein systems.

Peptide insertion, positioning, and stabilization in a membrane: insight from an all-atom molecular dynamics simulation.

Peptide insertion, positioning, and stabilization in a model membrane are probed via an all-atom molecular dynamics (MD) simulation. One peptide (WL5) is simulated in each leaflet of a solvated dimyristoylglycero-3-phosphate (DMPC) membrane. Within the first 5 ns, the peptides spontaneously insert into the membrane and then stabilize during the remaining 70 ns of simulation time. In both leaflets, the peptides localize to the membrane interface, and this localization is attributed to the formation of peptide-lipid hydrogen bonds. We show that the single tryptophan residue in each peptide contributes significantly to these hydrogen bonds; specifically, the nitrogen heteroatom of the indole ring plays a critical role. The tilt angles of the indole rings relative to the membrane normal in the upper and lower leaflets are approximately 26 degrees and 54 degrees , respectively. The tilt angles of the entire peptide chain are 62 degrees and 74 degrees . The membrane induces conformations of the peptide that are characteristic of beta-sheets, and the peptide enhances the lipid ordering in the membrane. Finally, the diffusion rate of the peptides in the membrane plane is calculated (based on experimental peptide concentrations) to be approximately 6 A(2)/ns, thus suggesting a 500 ns time scale for intermolecular interactions.

SR protein kinase 1 is resilient to inactivation.

SR protein kinase 1 (SRPK1) is a constitutively active kinase, which processively phosphorylates multiple serines within its substrates, ASF/SF2. We describe crystallographic, molecular dynamics, and biochemical results that shed light on how SRPK1 preserves its constitutive active conformation. Our structure reveals that unlike other known active kinase structures, the activation loop remains in an active state without any specific intraprotein interactions. Moreover, SRPK1 remains active despite extensive mutation to the activation segment. Molecular dynamics simulations reveal that SRPK1 partially absorbs the effect of mutations by forming compensatory interactions that maintain a catalytically competent chemical environment. Furthermore, SRPK1 is similarly resistant to deletion of its spacer loop region. Based upon a model of SRPK1 bound to a segment encompassing the docking motif and active-site peptide of ASF/SF2, we suggest a mechanism for processive phosphorylation and propose that the atypical resiliency we observed is critical for SRPK1's processive activity.

Molecular dynamics study of MscL interactions with a curved lipid bilayer.

Mechanosensitivity is a ubiquitous sensory mechanism found in living organisms. The simplest known mechanotransducing mechanism is found in bacteria in the form of the mechanosensitive membrane channel of large conductance, MscL. This channel has been studied extensively using a variety of methods at a functional and structural level. The channel is gated by membrane tension in the lipid bilayer alone. It serves as a safety valve protecting bacterial cells against hypoosmotic shock. MscL of Escherichia coli embedded in bilayers composed of asymmetric amounts of single-tailed and double-tailed lipids has been shown to gate spontaneously, even in the absence of membrane tension. To gain insight into the effect of the lipid membrane composition and geometry on MscL structure, a fully solvated, all-atom model of MscL in a stress-free curved bilayer composed of double- and single-tailed lipids was studied using a 9.5-ns molecular dynamics simulation. The bilayer was modeled as a domed structure accommodating the asymmetric composition of the monolayers. During the course of the simulation a spontaneous restructuring of the periplasmic loops occurred, leading to interactions between one of the loops and phospholipid headgroups. Previous experimental studies of the role of the loops agree with the observation that opening starts with a restructuring of the periplasmic loop, suggesting an effect of the curved bilayer. Because of limited resources, only one simulation of the large system was performed. However, the results obtained suggest that through the geometry and composition of the bilayer the protein structure can be affected even on short timescales.

Potentials of mean force for acetylcholine unbinding from the alpha7 nicotinic acetylcholine receptor ligand-binding domain.

The nicotinic acetylcholine receptor is a prototype ligand-gated ion channel that mediates signal transduction in the neuromuscular junctions and other cholinergic synapses. The molecular basis for the energetics of ligand binding and unbinding is critical to our understanding of the pharmacology of this class of receptors. Here, we used steered molecular dynamics to investigate the unbinding of acetylcholine from the ligand-binding domain of human alpha7 nicotinic acetylcholine receptor along four different predetermined pathways. Pulling forces were found to correlate well with interactions between acetylcholine and residues in the binding site during the unbinding process. From multiple trajectories along these unbinding pathways, we calculated the potentials of mean force for acetylcholine unbinding. Four available methods based on Jarzynski's equality were used and compared for their efficiencies. The most probable pathway was identified to be along a direction approximately parallel to the membrane. The derived binding energy for acetylcholine was in good agreement with that derived from the experimental binding constant for acetylcholine binding protein, but significantly higher than that for the complete human alpha7 nicotinic acetylcholine receptor. In addition, it is likely that several intermediate states exist along the unbinding pathways.

Dynamic binding of PKA regulatory subunit RI alpha.

Recent crystal structures have revealed that regulatory subunit RIalpha of PKA undergoes a dramatic conformational change upon complex formation with the catalytic subunit. Molecular dynamics studies were initiated to elucidate the contributions of intrinsic conformational flexibility and interactions with the catalytic subunit in formation and stabilization of the complex. Simulations of a single RIalpha nucleotide binding domain (NBD), missing cAMP, showed that its C helix spontaneously occupies two distinct conformations: either packed against the nucleotide binding domain as in its cAMP bound structure, or extended into an intermediate form resembling that of the holoenzyme structure. C helix extension was not seen in a simulation of either RIalpha NBD. In a model complex containing both NBDs and the catalytic subunit, well-conserved residues at the interface between the NBDs in the cAMP bound form were found to stabilize the complex through contacts with the catalytic subunit. The model structure is consistent with available experimental data.

Lipid bilayer pressure profiles and mechanosensitive channel gating.

The function of membrane proteins often depends on the proteins' interaction with their lipid environment, spectacularly so in the case of mechanosensitive channels, which are gated through tension mediated by the surrounding lipids. Lipid bilayer tension is distributed quite inhomogeneously, but neither the scale at which relevant variation takes place nor the effect of varying lipid composition or tension has yet been investigated in atomic detail. We calculated lateral pressure profile distributions in lipid bilayers of various composition from all-atom molecular dynamics simulations totaling 110.5 ns in length. Reproducible pressure profile features at the 1 A length scale were determined. Lipids with phosphatidylcholine headgroups were found to shift the lateral pressure out of the hydrophobic core and into the headgroup region by an amount that is independent of area per lipid. POPE bilayers simulated at areas smaller than optimal exerted dramatically higher lateral pressure in a narrow region at the start of the aliphatic chain. Stretching of POPC bilayers increased tension predominantly in the same region. A simple geometric analysis for the gating of the mechanosensitive channel MscL suggests that pressure profiles affect its gating through the second moment of the profile in a tension-independent manner.

Gating of MscL studied by steered molecular dynamics.

Steered molecular dynamics simulations of the mechanosensitive channel of large conductance, MscL, were used to investigate how forces arising from membrane tension induce gating of the channel. A homology model of the closed form of MscL from Escherichia coli was subjected to external forces of 35-70 pN applied to residues near the membrane-water interface. The magnitude and location of these forces corresponded to those determined from the lateral pressure profile computed from a lipid bilayer simulation. A fully expanded state was obtained on the 10-ns timescale that revealed the mechanism for transducing membrane forces into channel opening. The expanded state agrees well with proposed models of MscL gating, in that it entails an irislike expansion of the pore accompanied by tilting of the transmembrane helices. The channel was most easily opened when force was applied predominantly on the cytoplasmic side of MscL. Comparison of simulations in which gating progressed to varying degrees identified residues that pose steric hindrance to channel opening.