The C-terminus of CFAP410 forms a tetrameric helical bundle that is essential for its localization to the basal body.

Cilia are antenna-like organelles protruding from the surface of many cell types in the human body. Defects in ciliary structure or function often lead to diseases that are collectively called ciliopathies. Cilia and flagella-associated protein 410 (CFAP410) localizes at the basal body of cilia/flagella and plays essential roles in ciliogenesis, neuronal development and DNA damage repair. It remains unknown how its specific basal body location is achieved. Multiple single amino acid mutations in CFAP410 have been identified in patients with various ciliopathies. One of the mutations, L224P, is located in the C-terminal domain (CTD) of human CFAP410 and causes severe spondylometaphyseal dysplasia, axial (SMDAX). However, the molecular mechanism for how the mutation causes the disorder remains unclear. Here, we report our structural studies on the CTD of CFAP410 from three distantly related organisms, Homo sapiens, Trypanosoma brucei and Chlamydomonas reinhardtii. The crystal structures reveal that the three proteins all adopt the same conformation as a tetrameric helical bundle. Our work further demonstrates that the tetrameric assembly of the CTD is essential for the correct localization of CFAP410 in T. brucei, as the L224P mutation that disassembles the tetramer disrupts its basal body localization. Taken together, our studies reveal that the basal body localization of CFAP410 is controlled by the CTD and provide a mechanistic explanation for how the mutation L224P in CFAP410 causes ciliopathies in humans.

An Atomistic View on the Mechanism of Diatom Peptide-Guided Biomimetic Silica Formation.

Deciphering nature's remarkable way of encoding functions in its biominerals holds the potential to enable the rational development of nature-inspired materials with tailored properties. However, the complex processes that convert solution-state precursors into solid biomaterials remain largely unknown. In this study, an unconventional approach is presented to characterize these precursors for the diatom-derived peptides R5 and synthetic Silaffin-1A1 (synSil-1A1). These molecules can form defined supramolecular assemblies in solution, which act as templates for solid silica structures. Using a tailored structural biology toolbox, the structure-function relationships of these self-assemblies are unveiled. NMR-derived constraints are employed to enable a recently developed fractal-cluster formalism and then reveal the architecture of the peptide assemblies in atomistic detail. Finally, by monitoring the self-assembly activities during silica formation at simultaneous high temporal and residue resolution using real-time spectroscopy, the mechanism is elucidated underlying template-driven silica formation. Thus, it is demonstrated how to exercise morphology control over bioinorganic solids by manipulating the template architectures. It is found that the morphology of the templates is translated into the shape of bioinorganic particles via a mechanism that includes silica nucleation on the solution-state complexes' surfaces followed by complete surface coating and particle precipitation.

N6-methyladenosine in 5' UTR does not promote translation initiation.

The most abundant N6-methyladenosine (m6A) modification on mRNAs is installed non-stoichiometrically across transcripts, with 5' untranslated regions (5' UTRs) being the least conductive. 5' UTRs are essential for translation initiation, yet the molecular mechanisms orchestrated by m6A remain poorly understood. Here, we combined structural, biochemical, and single-molecule approaches and show that at the most common position, a single m6A does not affect translation yields, the kinetics of translation initiation complex assembly, or start codon recognition both under permissive growth and following exposure to oxidative stress. Cryoelectron microscopy (cryo-EM) structures of the late preinitiation complex reveal that m6A purine ring established stacking interactions with an arginine side chain of the initiation factor eIF2α, although with only a marginal energy contribution, as estimated computationally. These findings provide molecular insights into m6A interactions with the initiation complex and suggest that the subtle stabilization is unlikely to affect the translation dynamics under homeostatic conditions or stress.

Atomistic mechanism of the constitutive activation of PDGFRA via its transmembrane domain.

Single-point mutations in the transmembrane (TM) region of receptor tyrosine kinases (RTKs) can lead to abnormal ligand-independent activation. We use a combination of computational modeling, NMR spectroscopy and cell experiments to analyze in detail the mechanism of how TM domains contribute to the activation of wild-type (WT) PDGFRA and its oncogenic V536E mutant. Using a computational framework, we scan all positions in PDGFRA TM helix for identification of potential functional mutations for the WT and the mutant and reveal the relationship between the receptor activity and TM dimerization via different interfaces. This strategy also allows us design a novel activating mutation in the WT (I537D) and a compensatory mutation in the V536E background eliminating its constitutive activity (S541G). We show both computationally and experimentally that single-point mutations in the TM region reshape the TM dimer ensemble and delineate the structural and dynamic determinants of spontaneous activation of PDGFRA via its TM domain. Our atomistic picture of the coupling between TM dimerization and PDGFRA activation corroborates the data obtained for other RTKs and provides a foundation for developing novel modulators of the pathological activity of PDGFRA.

Effect of Oxidative Damage on the Stability and Dimerization of Superoxide Dismutase 1.

During their life cycle, proteins are subject to different modifications involving reactive oxygen species. Such oxidative damage to proteins may lead to the formation of insoluble aggregates and cytotoxicity and is associated with age-related disorders including neurodegenerative diseases, cancer, and diabetes. Superoxide dismutase 1 (SOD1), a key antioxidant enzyme in human cells, is particularly susceptible to such modifications. Moreover, this homodimeric metalloenzyme has been directly linked to both familial and sporadic amyotrophic lateral sclerosis (ALS), a devastating, late-onset motor neuronal disease, with more than 150 ALS-related mutations in the SOD1 gene. Importantly, oxidatively damaged SOD1 aggregates have been observed in both familial and sporadic forms of the disease. However, the molecular mechanisms as well as potential implications of oxidative stress in SOD1-induced cytotoxicity remain elusive. In this study, we examine the effects of oxidative modification on SOD1 monomer and homodimer stability, the key molecular properties related to SOD1 aggregation. We use molecular dynamics simulations in combination with thermodynamic integration to study microscopic-level site-specific effects of oxidative "mutations" at the dimer interface, including lysine, arginine, proline and threonine carbonylation, and cysteine oxidation. Our results show that oxidative damage of even single residues at the interface may drastically destabilize the SOD1 homodimer, with several modifications exhibiting a comparable effect to that of the most drastic ALS-causing mutations known. Additionally, we show that the SOD1 monomer stability decreases upon oxidative stress, which may lead to partial local unfolding and consequently to increased aggregation propensity. Importantly, these results suggest that oxidative stress may play a key role in development of ALS, with the mutations in the SOD1 gene being an additional factor.

Estimation of conformational entropy in protein-ligand interactions: a computational perspective.

Conformational entropy is an important component of the change in free energy upon binding of a ligand to its target protein. As a consequence, development of computational techniques for reliable estimation of conformational entropies is currently receiving an increased level of attention in the context of computational drug design. Here, we review the most commonly used techniques for conformational entropy estimation from classical molecular dynamics simulations. Although by-and-large still not directly used in practical drug design, these techniques provide a golden standard for developing other, computationally less-demanding methods for such applications, in addition to furthering our understanding of protein-ligand interactions in general. In particular, we focus on the quasi-harmonic approximation and discuss different approaches that can be used to go beyond it, most notably, when it comes to treating anharmonic and/or correlated motions. In addition to reviewing basic theoretical formalisms, we provide a concrete set of steps required to successfully calculate conformational entropy from molecular dynamics simulations, as well as discuss a number of practical issues that may arise in such calculations.

Microscopic analysis of protein oxidative damage: effect of carbonylation on structure, dynamics, and aggregability of villin headpiece.

One of the most important irreversible oxidative modifications of proteins is carbonylation, the process of introducing a carbonyl group in reaction with reactive oxygen species. Notably, carbonylation increases with the age of cells and is associated with the formation of intracellular protein aggregates and the pathogenesis of age-related disorders such as neurodegenerative diseases and cancer. However, it is still largely unclear how carbonylation affects protein structure, dynamics, and aggregability at the atomic level. Here, we use classical molecular dynamics simulations to study structure and dynamics of the carbonylated headpiece domain of villin, a key actin-organizing protein. We perform an exhaustive set of molecular dynamics simulations of a native villin headpiece together with every possible combination of carbonylated versions of its seven lysine, arginine, and proline residues, quantitatively the most important carbonylable amino acids. Surprisingly, our results suggest that high levels of carbonylation, far above those associated with cell death in vivo, may be required to destabilize and unfold protein structure through the disruption of specific stabilizing elements, such as salt bridges or proline kinks, or tampering with the hydrophobic effect. On the other hand, by using thermodynamic integration and molecular hydrophobicity potential approaches, we quantitatively show that carbonylation of hydrophilic lysine and arginine residues is equivalent to introducing hydrophobic, charge-neutral mutations in their place, and, by comparison with experimental results, we demonstrate that this by itself significantly increases the intrinsic aggregation propensity of both structured, native proteins and their unfolded states. Finally, our results provide a foundation for a novel experimental strategy to study the effects of carbonylation on protein structure, dynamics, and aggregability using site-directed mutagenesis.

Structure and dynamics of two beta-peptides in solution from molecular dynamics simulations validated against experiment.

We have studied two different beta-peptides in methanol using explicit solvent molecular dynamics simulations and the GROMOS 53A6 force field: a heptapeptide (peptide 1) expected to form a left-handed 3(14)-helix, and a hexapeptide (peptide 2) expected to form a beta-hairpin in solution. Our analysis has focused on identifying and analyzing the stability of the dominant secondary structure conformations adopted by the peptides, as well as on comparing the experimental NOE distance upper bounds and 3J-coupling values with their counterparts calculated on the basis of the simulated ensembles. Moreover, we have critically compared the present results with the analogous results obtained with the GROMOS 45A3 (peptide 1) and 43A1 (peptide 2) force fields. We conclude that within the limits of conformational sampling employed here, the GROMOS 53A6 force field satisfactorily reproduces experimental findings regarding the behavior of short beta-peptides, with accuracy that is comparable to but not exceeding that of the previous versions of the force field. GCE legend Conformational clustering analysis of the simulated ensemble of a ss-hexapeptide with two different simulation setups (a and b). The central members of all of the clusters populating more than 5% of all of the structures are shown, together with the most dominant hydrogen bonds and the corresponding percentages of cluster members containing them.

Mechanism and thermodynamics of binding of the polypyrimidine tract binding protein to RNA.

The polypyrimidine tract binding protein (PTB) is involved in many physiological processes, including alternative splicing, internal ribosomal entry side (IRES)-mediated initiation of translation, and polyadenylation, as well as in ensuring mRNA stability. However, the role of PTB in these processes is not fully understood, and this has motivated us to undertake a computational study of the protein. PTB RNA binding domains (RBDs) 3 and 4 and their complexes with oligopyrimidine RNAs were simulated using the GROMOS simulation software using the GROMOS 45A4 force field. First, the stability and fluctuations of the tertiary fold and of the secondary structural elements in individual domains, the combined RBD34 domain, and their complexes with RNA were studied. Second, the simulation results were validated against the experimental NMR NOE data. The analysis of hydrogen bonding patterns, salt bridge networks, and stacking interactions of the RNA to the binding pockets of the protein domains showed that binding is not sequence-specific and that many RNA fragments can bind to them successfully. Further calculations of the relative free energy of binding for different polypyrimidine sequences were carried out using the thermodynamic integration (TI) and single-step perturbation (SSP) methods. It is was not possible to calculate the relative free energies with high accuracy, but the obtained results do give qualitative insights into PTB's affinity for different RNA sequences. Furthermore, the low-energy conformations of the complexes that were found provided additional information about the mechanism of binding.

Simulated unfolded-state ensemble and the experimental NMR structures of villin headpiece yield similar wide-angle solution X-ray scattering profiles.

With the advent of powerful synchrotron sources, solution X-ray scattering is being increasingly used to get basic information about the structure of polypeptides. The solution scattering technique essentially provides one-dimensional data, which are then interpreted in terms of a three-dimensional structure through model building. Here we calculate wide-angle solution scattering patterns for an ensemble of simulated unfolded structures of villin headpiece, which differ from the native structure by rmsd = 8.8 +/- 1.0 A and have only negligible amounts of native secondary structure. We show that the wide-angle solution scattering pattern of such an ensemble shares significant similarity with the one based on the experimental NMR structures of the molecule. Our results suggest that solution scattering in the wide-angle limit, by itself, provides very little information about the secondary structure content of a polypeptide or its side-chain packing.

Comparing atomistic simulation data with the NMR experiment: how much can NOEs actually tell us?

Simulated molecular dynamics trajectories of proteins and nucleic acids are often compared with nuclear magnetic resonance (NMR) data for the purposes of assessing the quality of the force field used or, equally important, trying to interpret ambiguous experimental data. In particular, nuclear Overhauser enhancement (NOE) intensities or atom-atom distances derived from them are frequently calculated from the simulated ensembles because the distance restraints derived from NOEs are the key ingredient in NMR-based protein structure determination. In this study, we ask how diverse and nonnative-like an ensemble of structures can be and still match the experimental NOE distance upper bounds well. We present two examples in which simulated ensembles of highly nonnative polypeptide structures (an unfolded state ensemble of the villin headpiece and a high-temperature denatured ensemble of lysozyme) are shown to match fairly well the experimental NOE distance upper bounds from which the corresponding native structures were derived. For example, the unfolded ensemble of villin headpiece, which is on average 0.90 +/- 0.13 nm root-mean-square deviation away from the native NMR structure, deviates from the experimental restraints by only 0.027 nm on average. However, this artificially good agreement is largely a consequence of 1) the highly nonlinear effects of r(-6) (or r(-3)) averaging and 2) focusing only on the experimentally observed set of NOE bounds. Namely, in addition to the experimentally observed NOEs, both simulated ensembles (especially the villin ensemble) also predict a large number of NOEs, which are not seen in the experiment. If these are taken into account, the agreement between simulation and experiment gets markedly worse, as it should, given the nonnative nature of the underlying simulated ensembles. In light of the examples given, we conclude that comparing experimental NOE distance restraints with large simulated ensembles provides just by itself only limited information about the quality of simulation.

How large is an alpha-helix? Studies of the radii of gyration of helical peptides by small-angle X-ray scattering and molecular dynamics.

Using synchrotron radiation and the small-angle X-ray scattering technique we have measured the radii of gyration of a series of alanine-based alpha-helix-forming peptides of the composition Ace-(AAKAA)(n)-GY-NH(2), n=2-7, in aqueous solvent at 10(+/-1) degrees C. In contrast to other techniques typically used to study alpha-helices in isolation (such as nuclear magnetic resonance and circular dichroism), small-angle X-ray scattering reports on the global structure of a molecule and, as such, provides complementary information to these other, more sequence-local measuring techniques. The radii of gyration that we measure are, except for the 12-mer, lower than the radii of gyration of ideal alpha-helices or helices with frayed ends of the equivalent sequence-length. For example, the measured radius of gyration of the 37-mer is 14.2(+/-0.6)A, which is to be compared with the radius of gyration of an ideal 37-mer alpha-helix of 17.6A. Attempts are made to analyze the origin of this discrepancy in terms of the analytical Zimm-Bragg-Nagai (ZBN) theory, as well as distributed computing explicit solvent molecular dynamics simulations using two variants of the AMBER force-field. The ZBN theory, which treats helices as cylinders connected by random walk segments, predicts markedly larger radii of gyration than those measured. This is true even when the persistence length of the random walk parts is taken to be extremely short (about one residue). Similarly, the molecular dynamics simulations, at the level of sampling available to us, give inaccurate values of the radii of gyration of the molecules (by overestimating them by around 25% for longer peptides) and/or their helical content. We conclude that even at the short sequences examined here (< or =37 amino acid residues), these alpha-helical peptides behave as fluctuating semi-broken rods rather than straight cylinders with frayed ends.

Unusual compactness of a polyproline type II structure.

Polyproline type II (PPII) helix has emerged recently as the dominant paradigm for describing the conformation of unfolded polypeptides. However, most experimental observables used to characterize unfolded proteins typically provide only short-range, sequence-local structural information that is both time- and ensemble-averaged, giving limited detail about the long-range structure of the chain. Here, we report a study of a long-range property: the radius of gyration of an alanine-based peptide, Ace-(diaminobutyric acid)2-(Ala)7-(ornithine)2-NH2. This molecule has previously been studied as a model for the unfolded state of proteins under folding conditions and is believed to adopt a PPII fold based on short-range techniques such as NMR and CD. By using synchrotron radiation and small-angle x-ray scattering, we have determined the radius of gyration of this peptide to be 7.4 +/- 0.5 angstroms, which is significantly less than the value expected from an ideal PPII helix in solution (13.1 angstroms). To further study this contradiction, we have used molecular dynamics simulations using six variants of the AMBER force field and the GROMOS 53A6 force field. However, in all cases, the simulated ensembles underestimate the PPII content while overestimating the experimental radius of gyration. The conformational model that we propose, based on our small angle x-ray scattering results and what is known about this molecule from before, is that of a very flexible, fluctuating structure that on the level of individual residues explores a wide basin around the ideal PPII geometry but is never, or only rarely, in the ideal extended PPII helical conformation.

Random-coil behavior and the dimensions of chemically unfolded proteins.

Spectroscopic studies have identified a number of proteins that appear to retain significant residual structure under even strongly denaturing conditions. Intrinsic viscosity, hydrodynamic radii, and small-angle x-ray scattering studies, in contrast, indicate that the dimensions of most chemically denatured proteins scale with polypeptide length by means of the power-law relationship expected for random-coil behavior. Here we further explore this discrepancy by expanding the length range of characterized denatured-state radii of gyration (R(G)) and by reexamining proteins that reportedly do not fit the expected dimensional scaling. We find that only 2 of 28 crosslink-free, prosthetic-group-free, chemically denatured polypeptides deviate significantly from a power-law relationship with polymer length. The R(G) of the remaining 26 polypeptides, which range from 16 to 549 residues, are well fitted (r(2) = 0.988) by a power-law relationship with a best-fit exponent, 0.598 +/- 0.028, coinciding closely with the 0.588 predicted for an excluded volume random coil. Therefore, it appears that the mean dimensions of the large majority of chemically denatured proteins are effectively indistinguishable from the mean dimensions of a random-coil ensemble.

Structural correspondence between the alpha-helix and the random-flight chain resolves how unfolded proteins can have native-like properties.

Recently, we have proposed that, on average, the structure of the unfolded state of small, mostly alpha-helical proteins may be similar to the native structure (the 'mean-structure' hypothesis). After examining thousands of simulations of both the folded and the unfolded states of five polypeptides in atomistic detail at room temperature, we report here a result that seems at odds with the mean-structure hypothesis. Specifically, the average inter-residue distances in the collapsed unfolded structures agree well with the statistics of the ideal random-flight chain with link length of 3.8 A (the length of one amino acid). A possible resolution of this apparent contradiction is offered by the observation that the inter-residue distances in a typical alpha-helix over short stretches are close to the average distances in an ideal random-flight chain.

Atomistic protein folding simulations on the submillisecond time scale using worldwide distributed computing.

Atomistic simulations of protein folding have the potential to be a great complement to experimental studies, but have been severely limited by the time scales accessible with current computer hardware and algorithms. By employing a worldwide distributed computing network of tens of thousands of PCs and algorithms designed to efficiently utilize this new many-processor, highly heterogeneous, loosely coupled distributed computing paradigm, we have been able to simulate hundreds of microseconds of atomistic molecular dynamics. This has allowed us to directly simulate the folding mechanism and to accurately predict the folding rate of several fast-folding proteins and polymers, including a nonbiological helix, polypeptide alpha-helices, a beta-hairpin, and a three-helix bundle protein from the villin headpiece. Our results demonstrate that one can reach the time scales needed to simulate fast folding using distributed computing, and that potential sets used to describe interatomic interactions are sufficiently accurate to reach the folded state with experimentally validated rates, at least for small proteins.

The Trp cage: folding kinetics and unfolded state topology via molecular dynamics simulations.

Using over 75 mus of molecular dynamics simulation, we have generated several thousand folding simulations of the 20-residue Trp cage at experimental temperature and solvent viscosity. A total of 116 independent folding simulations reach RMSDcalpha values below 3 A RMSDcalpha, some as close as 1.4 A RMSDcalpha. We estimate a folding time of 5.5+/-3.5 mus, a rate that is in reasonable agreement with experimental kinetics. Finally, we characterize both the folded and unfolded ensemble under native conditions and note that the average topology of the unfolded ensemble is very similar to the topology of the native state.

Native-like mean structure in the unfolded ensemble of small proteins.

The nature of the unfolded state plays a great role in our understanding of proteins. However, accurately studying the unfolded state with computer simulation is difficult, due to its complexity and the great deal of sampling required. Using a supercluster of over 10,000 processors we have performed close to 800 micros of molecular dynamics simulation in atomistic detail of the folded and unfolded states of three polypeptides from a range of structural classes: the all-alpha villin headpiece molecule, the beta hairpin tryptophan zipper, and a designed alpha-beta zinc finger mimic. A comparison between the folded and the unfolded ensembles reveals that, even though virtually none of the individual members of the unfolded ensemble exhibits native-like features, the mean unfolded structure (averaged over the entire unfolded ensemble) has a native-like geometry. This suggests several novel implications for protein folding and structure prediction as well as new interpretations for experiments which find structure in ensemble-averaged measurements.