PubMed Health⌕ Search

Biomedical subjects

Federico Fogolari

Publications and source records attributed to Federico Fogolari.

18 recordsLinked to original sources
PubMed2024

Assessing nanobody interaction with SARS-CoV-2 Nsp9.

The interaction between SARS-CoV-2 non-structural protein Nsp9 and the nanobody 2NSP90 was investigated by NMR spectroscopy using the paramagnetic perturbation methodology PENELOP (Paramagnetic Equilibrium vs Nonequilibrium magnetization Enhancement or LOss Perturbation). The Nsp9 monomer is an essential component of the replication and transcription complex (RTC) that reproduces the viral gRNA for subsequent propagation. Therefore preventing Nsp9 recruitment in RTC would represent an efficient antiviral strategy that could be applied to different coronaviruses, given the Nsp9 relative invariance. The NMR results were consistent with a previous characterization suggesting a 4:4 Nsp9-to-nanobody stoichiometry with the occurrence of two epitope pairs on each of the Nsp9 units that establish the inter-dimer contacts of Nsp9 tetramer. The oligomerization state of Nsp9 was also analyzed by molecular dynamics simulations and both dimers and tetramers resulted plausible. A different distribution of the mapped epitopes on the tetramer surface with respect to the former 4:4 complex could also be possible, as well as different stoichiometries of the Nsp9-nanobody assemblies such as the 2:2 stoichiometry suggested by the recent crystal structure of the Nsp9 complex with 2NSP23 (PDB ID: 8dqu), a nanobody exhibiting essentially the same affinity as 2NSP90. The experimental NMR evidence, however, ruled out the occurrence in liquid state of the relevant Nsp9 conformational change observed in the same crystal structure.

Viral Nonstructural Proteins↗
PubMed2006

Collagen plays an active role in the aggregation of beta2-microglobulin under physiopathological conditions of dialysis-related amyloidosis.

Dialysis-related amyloidosis is characterized by the deposition of insoluble fibrils of beta(2)-microglobulin (beta(2)-m) in the musculoskeletal system. Atomic force microscopy inspection of ex vivo amyloid material reveals the presence of bundles of fibrils often associated to collagen fibrils. Aggregation experiments were undertaken in vitro with the aim of reproducing the physiopathological fibrillation process. To this purpose, atomic force microscopy, fluorescence techniques, and NMR were employed. We found that in temperature and pH conditions similar to those occurring in periarticular tissues in the presence of flogistic processes, beta(2)-m fibrillogenesis takes place in the presence of fibrillar collagen, whereas no fibrils are obtained without collagen. Moreover, the morphology of beta(2)-m fibrils obtained in vitro in the presence of collagen is extremely similar to that observed in the ex vivo sample. This result indicates that collagen plays a crucial role in beta(2)-m amyloid deposition under physiopathological conditions and suggests an explanation for the strict specificity of dialysis-related amyloidosis for the tissues of the skeletal system. We hypothesize that positively charged regions along the collagen fiber could play a direct role in beta(2)-m fibrillogenesis. This hypothesis is sustained by aggregation experiments performed by replacing collagen with a poly-L-lysine-coated mica surface. As shown by NMR measurements, no similar process occurs when poly-L-lysine is dissolved in solution with beta(2)-m. Overall, the findings are consistent with the estimates resulting from a simplified collagen model whereby electrostatic effects can lead to high local concentrations of oppositely charged species, such as beta(2)-m, that decay on moving away from the fiber surface.

Amyloid↗
PubMed2006

NMR dynamic studies suggest that allosteric activation regulates ligand binding in chicken liver bile acid-binding protein.

Apo chicken liver bile acid-binding protein has been structurally characterized by NMR. The dynamic behavior of the protein in its apo- and holo-forms, complexed with chenodeoxycholate, has been determined via (15)N relaxation and steady state heteronuclear (15)N((1)H) nuclear Overhauser effect measurements. The dynamic parameters were obtained at two pH values (5.6 and 7.0) for the apoprotein and at pH 7.0 for the holoprotein, using the model free approach. Relaxation studies, performed at three different magnetic fields, revealed a substantial conformational flexibility on the microsecond to millisecond time scales, mainly localized in the C-terminal face of the beta-barrel. The observed dynamics are primarily caused by the protonation/deprotonation of a buried histidine residue, His(98), located on this flexible face. A network of polar buried side chains, defining a spine going from the E to J strand, is likely to provide the long range connectivity needed to communicate motion from His(98) to the EF loop region. NMR data are accompanied by molecular dynamics simulations, suggesting that His(98) protonation equilibrium is the triggering event for the modulation of a functionally important motion, i.e. the opening/closing at the protein open end, whereas ligand binding stabilizes one of the preexisting conformations (the open form). The results presented here, complemented with an analysis of proteins belonging to the intracellular lipid-binding protein family, are consistent with a model of allosteric activation governing the binding mechanism. The functional role of this mechanism is thoroughly discussed within the framework of the mechanism for the enterohepatic circulation of bile acids.

Allosteric Regulation↗
PubMed2006

Molecular analysis of a human PAX6 homeobox mutant.

Pax6 controls eye, pancreas and brain morphogenesis. In humans, heterozygous PAX6 mutations cause aniridia and various other congenital eye abnormalities. Most frequent PAX6 missense mutations are located in the paired domain (PD), while very few missense mutations have been identified in the homeodomain (HD). In the present report, we describe a molecular analysis of the human PAX6 R242T missense mutation, which is located in the second helix of the HD. It was identified in a male child with partial aniridia in the left eye, presenting as a pseudo-coloboma. Gel-retardation assays revealed that the mutant HD binds DNA as well as the wild-type HD. In addition, the mutation does not modify the DNA-binding properties of the PD. Cell transfection assays indicated that the steady-state levels of the full length mutant protein are higher than those of the wild-type one. In cotransfection assays a PAX6 responsive promoter is activated to a higher extent by the mutant protein than by the wild-type protein. In vitro limited proteolysis assays indicated that the presence of the mutation reduces the sensitivity to trypsin digestion. Thus, we suggest that the R242T human phenotype could be due to abnormal increase of PAX6 protein, in keeping with the reported sensitivity of the eye phenotype to increased PAX6 dosage.

Animals↗
PubMed2005

A decoy set for the thermostable subdomain from chicken villin headpiece, comparison of different free energy estimators.

BACKGROUND: Estimators of free energies are routinely used to judge the quality of protein structural models. As these estimators still present inaccuracies, they are frequently evaluated by discriminating native or native-like conformations from large ensembles of so-called decoy structures. RESULTS: A decoy set is obtained from snapshots taken from 5 long (100 ns) molecular dynamics (MD) simulations of the thermostable subdomain from chicken villin headpiece. An evaluation of the energy of the decoys is given using: i) a residue based contact potential supplemented by a term for the quality of dihedral angles; ii) a recently introduced combination of four statistical scoring functions for model quality estimation (FRST); iii) molecular mechanics with solvation energy estimated either according to the generalized Born surface area (GBSA) or iv) the Poisson-Boltzmann surface area (PBSA) method. CONCLUSION: The decoy set presented here has the following features which make it attractive for testing energy scoring functions:1) it covers a broad range of RMSD values (from less than 2.0 A to more than 12 A);2) it has been obtained from molecular dynamics trajectories, starting from different non-native-like conformations which have diverse behaviour, with secondary structure elements correctly or incorrectly formed, and in one case folding to a native-like structure. This allows not only for scoring of static structures, but also for studying, using free energy estimators, the kinetics of folding;3) all structures have been obtained from accurate MD simulations in explicit solvent and after molecular mechanics (MM) energy minimization using an implicit solvent method. The quality of the covalent structure therefore does not suffer from steric or covalent problems. The statistical and physical effective energy functions tested on the set behave differently when native simulation snapshots are included or not in the set and when averaging over the trajectory is performed.

Animals↗
PubMed2005

Functional analysis of a novel RUNX2 missense mutation found in a family with cleidocranial dysplasia.

Mutations of the RUNX2 gene result in dominantly inherited cleidocranial dysplasia (CCD). RUNX2 encodes for an osteoblast-specific transcription factor, which recognizes specific DNA sequences by the runt domain. DNA binding is stabilized by the interaction with the protein CBFbeta, which induces structural modifications of the runt domain. A novel 574G > A RUNX2 missense mutation has been found in members of a family clinically diagnosed with CCD. This mutation causes the glycine at position 192 to change to arginine (G192R), in loop 9 of the runt domain. Unlike other residues of loop 9, G192 does not establish DNA contacts. Accordingly, the G192R mutant showed a 50% reduction in binding activity compared to the wild-type runt domain. However, the mutation completely abolished the activating properties of the protein on osteocalcin promoter. Moreover, the G192R mutant exerts a dominant-negative effect when overexpressed. Computer modeling indicated that the G192R mutation perturbs not only loop 9, but also other parts of the runt domain, suggesting impairment of the interaction with CBFbeta.

Adult↗
PubMed2005

Solution structure of beta(2)-microglobulin and insights into fibrillogenesis.

The solution structure of human beta(2)-microglobulin (beta(2)-m) was determined by (1)H NMR spectroscopy and restrained modeling calculations. Compared to the crystal structure of type I major histocompatibility complex (MHC-I), where the protein is associated to the heavy-chain component, several differences are observed, i.e., increased separation between strands A and B, displacements of strand C' and loop DE, shortening of strands D and E. These modifications can be considered as the prodromes of the amyloid transition. Even minor charge changes in response to pH, as is the case with H31 imidazole protonation, trigger the transition that starts with unpairing of strand A. The same mechanism accounts for the partial unfolding and fiber formation subsequent to Cu(2+) binding which is shown to occur primarily at H31. Solvation of the protected regions in MHC-I decreases the tertiary packing by breaking the contiguity of the surface hydrophobic patches via surface charge cluster. Mutants or truncated forms of beta(2)-m can be designed to remove the instability from H31 titration or to enhance the instability through surface charge suppression. By monitoring the conformational evolution of wild-type protein and variants thereof, either in response or absence of external perturbation, valuable insights into intermediate structure and fibrillogenesis mechanisms are gained.

Amyloid↗
PubMed2005

Comparative molecular dynamics simulations of HIV-1 integrase and the T66I/M154I mutant: binding modes and drug resistance to a diketo acid inhibitor.

HIV-1 IN is an essential enzyme for viral replication and an interesting target for the design of new pharmaceuticals for use in multidrug therapy of AIDS. L-731,988 is one of the most active molecules of the class of beta-diketo acids. Individual and combined mutations of HIV-1 IN at residues T66, S153, and M154 confer important degrees of resistance to one or more inhibitors belonging to this class. In an effort to understand the molecular mechanism of the resistance of T66I/M154I IN to the inhibitor L-731,988 and its specific binding modes, we have carried out docking studies, explicit solvent MD simulations, and binding free energy calculations. The inhibitor was docked against different protein conformations chosen from prior MD trajectories, resulting in 2 major orientations within the active site. MD simulations have been carried out for the T66I/M154I DM IN, DM IN in complex with L-731,988 in 2 different orientations, and 1QS4 IN in complex with L-731,988. The results of these simulations show a similar dynamical behavior between T66I/M154I IN alone and in complex with L-731,988, while significant differences are observed in the mobility of the IN catalytic loop (residues 138-149). Water molecules bridging the inhibitor to residues from the active site have been identified, and residue Gln62 has been found to play an important role in the interactions between the inhibitor and the protein. This work provides information about the binding modes of L-731,988, as well as insight into the mechanism of inhibitor-resistance in HIV-1 integrase.

Amino Acid Substitution↗
PubMed2005

MM/PBSA analysis of molecular dynamics simulations of bovine beta-lactoglobulin: free energy gradients in conformational transitions?

The pH-driven opening and closure of beta-lactoglobulin EF loop, acting as a lid and closing the internal cavity of the protein, has been studied by molecular dynamics (MD) simulations and free energy calculations based on molecular mechanics/Poisson-Boltzmann (PB) solvent-accessible surface area (MM/PBSA) methodology. The forms above and below the transition pH differ presumably only in the protonation state of residue Glu89. MM/PBSA calculations are able to reproduce qualitatively the thermodynamics of the transition. The analysis of MD simulations using a combination of MM/PBSA methodology and the colony energy approach is able to highlight the driving forces implied in the transition. The analysis suggests that global rearrangements take place before the equilibrium local conformation is reached. This conclusion may bear general relevance to conformational transitions in all lipocalins and proteins in general.

Animals↗
PubMed2005

Application of MM/PBSA colony free energy to loop decoy discrimination: toward correlation between energy and root mean square deviation.

Accurate free energy estimation is needed in many predictive tasks. The molecular mechanics/Poisson-Boltzmann solvent accessible surface area (MM/PBSA) approach has proven to be accurate. However, the correlation between the estimated free energy and the distance (e.g., root mean square deviation [RMSD]) from the most stable conformation is hindered by the strong free energy dependence on minor conformational variations. In this paper, a protocol for MM/PBSA free energy estimation is designed and tested on several loop decoy sets. We show that further integration of MM/PBSA free energy estimator with the colony energy approach makes the correlation between the free energy and RMSD from the native structure apparent, for the test sets on which it could be applied. Our results suggest that (1) the MM/PBSA free energy estimator is able to detect native-like structures for most decoy sets, and (2) application of the colony energy approach greatly hampers the MM/energy strong dependence on minor conformational changes.

Amino Acids↗
PubMed2004

Constraint Logic Programming approach to protein structure prediction.

BACKGROUND: The protein structure prediction problem is one of the most challenging problems in biological sciences. Many approaches have been proposed using database information and/or simplified protein models. The protein structure prediction problem can be cast in the form of an optimization problem. Notwithstanding its importance, the problem has very seldom been tackled by Constraint Logic Programming, a declarative programming paradigm suitable for solving combinatorial optimization problems. RESULTS: Constraint Logic Programming techniques have been applied to the protein structure prediction problem on the face-centered cube lattice model. Molecular dynamics techniques, endowed with the notion of constraint, have been also exploited. Even using a very simplified model, Constraint Logic Programming on the face-centered cube lattice model allowed us to obtain acceptable results for a few small proteins. As a test implementation their (known) secondary structure and the presence of disulfide bridges are used as constraints. Simplified structures obtained in this way have been converted to all atom models with plausible structure. Results have been compared with a similar approach using a well-established technique as molecular dynamics. CONCLUSIONS: The results obtained on small proteins show that Constraint Logic Programming techniques can be employed for studying protein simplified models, which can be converted into realistic all atom models. The advantage of Constraint Logic Programming over other, much more explored, methodologies, resides in the rapid software prototyping, in the easy way of encoding heuristics, and in exploiting all the advances made in this research area, e.g. in constraint propagation and its use for pruning the huge search space.

Computational Biology↗
PubMed2003

NMR solution structure of viscotoxin C1 from Viscum album species Coloratum ohwi: toward a structure-function analysis of viscotoxins.

The high resolution three-dimensional structure of the newly discovered plant viscotoxin C1, from the Asiatic Viscum album ssp. Coloratum ohwi, has been determined in solution by (1)H NMR spectroscopy at pH 3.6 and 285 K. The viscotoxin C1-fold, consisting of a helix-turn-helix motif and a short stretch of an antiparralel beta-sheet is very similar to that found for the highly similar viscotoxins A2 and A3 and for other related thionins. Different functional properties of members of the thionin family are discussed here in light of the structural and electrostatic properties. Among the very homologous family of alpha- and beta-thionins, known for their antimicrobial activity, the viscotoxin subfamily differs from the other members because of its high toxicity against tumoral cells. Key residues for the modulation of viscotoxin cytotoxicity have been identified on the basis of sequence and structural alignment.

Amino Acid Sequence↗
PubMed2003

EF loop conformational change triggers ligand binding in beta-lactoglobulins.

Beta-lactoglobulins, belonging to the lipocalin family, are a widely studied group of proteins, characterized by the ability to solubilize and transport hydrophobic ligands, especially fatty acids. Despite many reports, the mechanism of ligand binding and the functional role of these proteins is still unclear, and many contradicting concepts are often encountered in the literature. In the present paper the comparative analysis of the binding properties of beta-lactoglobulins has been performed using sequence-derived information, structure-based electrostatic calculations, docking simulations, and NMR experiments. Our results reveal for the first time the mechanism of beta-lactoglobulin ligand binding, which is completely determined by the opening-closing of EF loop, triggered by Glu89 protonation. The alkaline shift observed for Glu89 pKa in porcine beta-lactoglobulin (pKa 9.7) with respect to the bovine species (pKa 5.5) depends upon the interplay of electrostatic effects of few nearby key residues. Porcine protein is therefore able to bind fatty acids provided that the appropriate pH solution conditions are met (pH > 8.6), where the EF loop conformational change can take place. The unusually high pH of binding detected for porcine beta-lactoglobulin seems to be functional to lipases activity. Theoretical pKa calculations extended to representative beta-lactoglobulins allowed the identification of key residues involved in structurally and functionally important electrostatic interactions. The results presented here provide a strong indication that the described conformational change is a common feature of all beta-lactoglobulins.

Amino Acid Sequence↗
PubMed2003

Amino acid empirical contact energy definitions for fold recognition in the space of contact maps.

BACKGROUND: Contradicting evidence has been presented in the literature concerning the effectiveness of empirical contact energies for fold recognition. Empirical contact energies are calculated on the basis of information available from selected protein structures, with respect to a defined reference state, according to the quasi-chemical approximation. Protein-solvent interactions are estimated from residue solvent accessibility. RESULTS: In the approach presented here, contact energies are derived from the potential of mean force theory, several definitions of contact are examined and their performance in fold recognition is evaluated on sets of decoy structures. The best definition of contact is tested, on a more realistic scenario, on all predictions including sidechains accepted in the CASP4 experiment. In 30 out of 35 cases the native structure is correctly recognized and best predictions are usually found among the 10 lowest energy predictions. CONCLUSION: The definition of contact based on van der Waals radii of alpha carbon and side chain heavy atoms is seen to perform better than other definitions involving only alpha carbons, only beta carbons, all heavy atoms or only backbone atoms. An important prerequisite for the applicability of the approach is that the protein structure under study should not exhibit anomalous solvent accessibility, compared to soluble proteins whose structure is deposited in the Protein Data Bank. The combined evaluation of a solvent accessibility parameter and contact energy allows for an effective gross screening of predictive models.

Amino Acids↗
PubMed2003

Investigations of Sso7d catalytic residues by NMR titration shifts and electrostatic calculations.

Sso7d is a small basic protein consisting of 62 amino acids isolated from the thermoacidophilic archeobacterium Sulfolobus solfataricus. The protein is endowed with DNA binding properties, RNase activity, and the capability of rescuing aggregated proteins in the presence of ATP. In this study, the electrostatic properties of Sso7d are investigated by using the Poisson-Boltzmann calculation of the surface potential distribution and following by NMR spectroscopy the proton chemical shift pH titration of acidic residues. Although the details of the catalytic mechanism still have to be defined, the results from NMR experiments confirm the possible involvement of Glu35 as the proton acceptor in the catalytic reaction, as seen by its abnormally high pK(a) value. Poisson-Boltzmann calculations and NMR titration shifts suggest the presence of a possible hydrogen bond between Glu35 and Tyr33, with a consequent rather rigid arrangement at these positions. Comparison with RNase T1 suggests that Tyr7 may be a good candidate for acting as a proton donor in the active site of Sso7d as shown by its low phenolic pK(a) of approximately 9.3. Titration experiments performed with the UpA, a RNA dinucleotide model, showed that the protein residues affected by the interaction are mainly located in a different region with respect to the surface affected by DNA recognition, in good agreement with the surface potential distribution found with electrostatic calculations.

Archaeal Proteins↗
PubMed2003

Protocol for MM/PBSA molecular dynamics simulations of proteins.

Continuum solvent models have been employed in past years for understanding processes such as protein folding or biomolecular association. In the last decade, several attempts have been made to merge atomic detail molecular dynamics simulations with solvent continuum models. Among continuum models, the Poisson-Boltzmann solvent accessible surface area model is one of the oldest and most fundamental. Notwithstanding its wide usage for simulation of biomolecular electrostatic potential, the Poisson-Boltzmann equation has been very seldom used to obtain solvation forces for molecular dynamics simulation. We propose here a fast and reliable methodology to implement continuum forces in standard molecular mechanics and dynamics algorithms. Results for a totally unrestrained 1 ns molecular dynamics simulation of a small protein are quantitatively similar to results obtained by explicit solvent molecular dynamics simulations.

Computer Simulation↗
PubMed2002

Peptide models of folding initiation sites of bovine beta-lactoglobulin: identification of nativelike hydrophobic interactions involving G and H strands.

In an attempt to characterize the early folding events in bovine beta-lactoglobulin (BLG), a set of peptides, covering the flexible N-terminal region and the stable C-terminus beta-core, was synthesized and analyzed by circular dichroism and by nuclear magnetic resonance in water, trifluoroethanol (TFE), and sodium dodecyl sulfate (SDS) below and above the critical micellar concentration. The role of local and long-range hydrophobic interactions in guiding the folding has been investigated. For the peptide fragment covering the more flexible N-terminal region of BLG (beta-strands A, B), where both theoretical predictions and kinetic refolding experiments suggested the formation of non-native alpha-helix, no native long-range contacts were identified, and a helical secondary structure was stabilized only in the presence of 25 mM SDS. At variance, in 50% (v/v) TFE, native, long-range hydrophobic interactions were observed in the peptide covering the core region comprising G and H beta-strands. The side chains involved in these interactions form a nativelike hydrophobic cluster, thus suggesting that the GH region may act as the folding initiation site for BLG. This result is reinforced by the identification, in the urea denaturated BLG, of residual structure located at the level of the GH interface, as evidenced by NMR analysis. These results, in excellent agreement with kinetic, thermodynamic, and cold denaturation folding data, once more underline the utmost importance of the GH region for the stability and folding of BLG. Severe aggregation effects prevented the structural analysis of the peptide covering the EFGH region, indicating that this larger segment does not represent an independent folding domain and that the terminal alpha-helix is necessary for stabilizing the BLG folding core.

Amino Acid Sequence↗